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RNA-Seq of the Catharanthus roseus SRA database was done in order to detect putative 
universal stress proteins (USPs) and their possible controlling factors. Previous analysis 
indicated the existence and characterization of uspA-like genes. In silico analysis of RNA- 
Seq database in several plant tissues revealed the possible functions and régulations of 
some uspA-Yike transcripts whose transcription factors (TFs) that might drive their 
expression were detected. BLAST indicated the existence of TF superfamilies erf (ethylene- 
responsive TF), bHLH (basic helix-loop-helix) and WRKY that might regulate several uspA- 
like genes. This data was proven via semi-quantitative RT-PCR in four plant tissues. Several 
of these transcription factor superfamilies are known for their action in the plant defense 
against biotic and abiotic stresses. 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 
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ML, mature leaf 
FL, flower 
R, root 

MTI, MAMP-triggered immunity 
PTI, PAMP-triggered immunity 
ETI, effector-triggered immunity 
Trp, tryptophan 
JA, jasmonic acid 
SA, salicylic acid 


1. Introduction 

Vinca or Catharanthus roseus is an important médicinal 
plant species with several applications in pharmaceuticals 
and industry. This plant species is well known for the 
biosynthesis of two anticancer bisindole alkaloids, namely 
Vinblastine and vincristine, produced during the terpenoid 
indole alkaloids (TIAs) pathway [1,2]. This wild plant 
species can tolerate several biotic and abiotic stresses 
[3,4]; hence, we screened its transcriptome for usp-like 
gene family that might act in conferring tolérance to such 
adverse conditions. Although the exact function of this 
gene family is unknown, it acts not only through improving 
tolérance to biotic and abiotic stresses, but also tolérance 
to other stresses, e.g., starvation and DNA damage [5j. 

The RNA-Seq data of different tissues that were 
deposited in the CathaCyc database corne from a study 
originally conducted by Van Moerkercke et al. [6[. This 
database was previously utilized [7] for detecting and 
characterizing 24 uspA-like genes. Transcription factors 
(TFs) with known actions that might drive these uspA-like 
genes likely indicate the mode of action of these genes. TFs 
are key proteins required in the régulation of almost ail 
biosynthetic pathways in life [8[. They bind to spécifie 
target DNA sequences controlling cellular transcription by 
either promoting or blocking the recruitment of RNA 
polymerase to the target genes. TFs usually bind directiy to 
the target promoter régions; however, they sometimes 
bind to RNA polymerase. The action of TF directiy affect 
downstream genes functioning at a particular time on 
target cells (gene-in-time/gene-in-site). They are also 
important in driving the plant’s responses to biotic and 
abiotic stresses [9], 

In the présent study, an in silico analysis was conducted 
in order to detect the TFs that might regulate the 
expression of uspA-like genes in C. roseus in a trial to 
detect the possible rôle of these uspA-like genes under 
biotic and abiotic stresses. 

2. Materials and methods 

An RNA-seq analysis of C. roseus available at NCBI 
( http://www.ncbi.nlm.nih.gov/sra , experiment SRP005953) 
was conducted in order to detect TFs concordantly 
expressed with uspA-like genes that had been previously 
characterized [7], The data were de novo-assembled using 
Trinity-RNAseq v (r2013_08_14), and the recovered trans- 
cripts were calculated for abundance and differential 
expression using the Rsem package and the edgR method 
[7], Blast alignment was done as previously described [7], 


and the values of fpkm for the recovered uspA-like 
transcripts and their co-expressed transcripts were calcu- 
lated. The resulted transcripts were clustered based on their 
expression patterns across different tissues, treatments and 
hairy root génotypes. Clusters were studied in order to 
detect TF transcripts co-expressed with uspA-like trans- 
cripts in C. roseus. A structural analysis was done for 
selected C. roseus TF transcripts frequently co-expressed 
with uspA-like transcripts. 

To validate the RNA-Seq dataset of C. roseus, semi- 
quantitative RT-PCR (or sqRT-PCR) was conducted for 
seven transcripts of selected clusters. They included four 
uspA-like genes and three co-expressing TFs. RNAs were 
isolated from four pot-grown plant tissues (e.g., mature 
leaf, flower, root and stem) using the Trizol reagent method 
(Invitrogen) and treated with RNase-free DNase (Promega 
Inc.). Primers were designed using Netprimer software 
(http://www.premierbiosoft.com/netprimer/index.html) 
with the following criteria: length of 20-27 bases, GC 
content of ~50%, with minimal secondary structure, 
comparable annealing températures (48-55 °C) of the 
primer pairs, and PCR products of 281-300 bp (Table SI). 
sqRT-PCR was performed in forty cycles including déna- 
turation at 94 C for 15 s, annealing at 47-48 C for 30 s, 
and extension at 72 °C for 45 s. Amounts of 1 ul cDNA, 
1 x PCR buffer (with 1.5 mM MgCl 2 ), 200 uM dNTPs, 
200 nM of each gene-specifïc primers and 0.2 U of Taq DNA 
polymerase (Promega Inc.) were used in the reaction. 
Amplicons were run on a 1.2% agarose gels stained with 
ethidium bromide and visualized using the Gel Doc XR 
from Bio-Rad Laboratories (Hercules, CA, <). The actin gene 
(250 bp) was used as the house-keeping control. 

3. Results and discussion 

RNA-Seq analysis of the C. roseus SRA database of 
different tissues resulted in the recovery of 50,723 
transcripts that were calculated for abundance and 
differential expression. A number of 24 uspA-like trans- 
cripts were previously detected and characterized [7[. The 
in silico analysis of the présent study resulted in the 
recovery of 1245 transcripts involving 13, out of the 24, 
uspA-like transcripts (Table S2). The rest of the transcripts 
were not regulated in different tissues of C. roseus, hence 
not presented in the clusters generated from the analysis. 
The regulated transcripts were shown in 12 clusters based 
on their expression patterns across different tissues. BLAST 
to detect hits in the NCBI indicated the occurrence of a 
large number of TFs mainly in the erf (ethylene-responsive 
TF), bHLH (basic helix-loop-helix) and WRKY gene 
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Table 1 

Transcription factors in Catharanthus roseus assemblies and their NCBI protein and mRNA hits. 
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ID 

Transcript ID 

NCBI protein hit 

TF1 

comp30166_c0 

XP_004 162335 

TF2 

comp24014_c0 

XP_0041 64936 

TF3 

comp30976_c0 

XP_004 173262 

TF4 

compl6043_c0 

XP_004 169980 

TF5 

compl0162_c0 

XP_003521193 

TF6 

comp24246_c0 

XP_004250124 

TF7 

compl5693_c0 

XP_006602463 

TF8 

comp20416_c0 

XP_01 1460603 

TF9 

comp2848_c0 

XP_006351643 

TF10 

comp4341_c0 

XP_006356833 

TF11 

compl0042_c0 

AJF11723 

TF12 

compll785_c0 

XP_003 534796 

TF13 

comp3515_cl 

XP_004235231 

TF14 

comp20639_c0 

XP_0041 65456 

TF15 

comp21556_c0 

XP_004504739 

TF16 

comp26488_c0 

XP_004501445 

TF17 

comp8026_cl 

XP_004291261 


Description 

Cucumis sativus ethylene-responsive transcription factor erf034-like mRNA 
Cucumis sativus ethylene-responsive transcription factor 5-like mRNA 
Cucumis sativus ap2-like ethylene-responsive transcription factor bbm2-like mRNA 
Cucumis sativus ethylene-responsive transcription factor erf008-like mRNA 
Glycine max transcription factor bhlh48-like mRNA 
Solarium lycopersicum transcription factor bhlh96-like mRNA 
Glycine max transcription factor bhlh78-like transcript variant mRNA 
Transcription factor bHLH106-like [ Fragaria vesca subsp. vesca] 

Solanum tuberosum transcription factor bhlhl55-like transcript variant mRNA 

Solanum tuberosum transcription factor bhlh74-like transcript variant mRNA 

Citrus sinensis probable wrky transcription factor 40-like mRNA 

Glycine max probable wrky transcription factor 14-like transcript variant mRNA 

Solanum lycopersicum probable wrky transcription factor 70-like mRNA 

Cucumis sativus wrky transcription factor 72-like mRNA 

Cicer arietinum wrky transcription factor 9-like transcript variant mRNA 

Cicer arietinum wrky transcription factor 50-like mRNA 

Fragaria vesca wrky transcription factor 11-like mRNA 


TF1-4: transcription factors with erf domain, TF5-10: transcription factors with bFTLH domain, TF11-17: transcription factors with WRKY domain. 


superfamilies. A number of 17 (four of erf, six of bHLH and 
seven of WRKY) TFs (see Table 1 for details) showed co- 
expression with six, out of the 13, uspA-like transcripts in 
five clusters, e.g., 4, 7, 8, 9 and 12 (Fig. 1 & Table S2). In 
cluster 4, the uspA-like transcript KJ634217 co-expressed 
with bHLH48. In cluster 7, the uspA-like transcripts 
KJ634222 and KJ634224 were in concordance with 
WRKY40, while the uspA-like transcript KJ634231 co- 
expressed with erf034. The uspA-like transcripts KJ634232 
and KJ634221 co-expressed with a number of ten and five 
TFs of the three families, respectively. sqRT-PCR was used 
for validating the RNA-Seq dataset utilizing the selected 
co-expressing uspA-like and TF transcripts of clusters 4, 
7 and 8 with varying expression levels in the four selected 
tissues, e.g., immature leaf (iML), mature leaf (ML), flower 
(FL), and root (R). The results concerning sqRT-PCR shown 
in Fig. 2 indicated an alignment with the RNA-Seq dataset 
(Table S2 and Fig. 1). 

The AP-2 - or activating protein-2 of ethylene respon- 
sive (ER) TF - domain represents a superfamily of closely 
related TFs, which plays a critical rôle in regulating gene 
expression during the early development [10]. It was also 
shown to interfère with other signal transduction path- 
ways [11,12], Ethylene-responsive element-binding fac- 
tors (ERFs) are a superfamily of transcriptional factors [13] 
that induce many developmental processes, and are also 
involved in the adaptation to biotic or abiotic stresses such 
as pathogen attack, wounding, UV irradiation, extreme 
température, and drought [14-16]. This TF is a key element 
in the intégration of both signais for the régulation of 
defense response genes [17]. Recent studies hâve extended 
the influence of ethylene signaling with plant responses to 
both sait and water stress [18,19], The erf8 and erf34 were 
reported in Arabidopsis and lice to be involved in the 
régulation of gene expression by stress factors and 
components of stress signal transduction pathways 
[13], The erf was further reported to operate during the 
progression of leaf senescence [20], while erf5 was 
reported in Arabidopsis to regulate the chitin-induced 
innate immunity response [21]. Similar results were 


reached in tomato ( Solanum lycopersicum ), tobacco 
(Nicotiana tabacum ), and Arabidopsis [22,23], 

The information available for the basic helix-loop-helix 
(bHLH) superfamily indicates its existence in eukaryotes 
with numerous functions [24,25], It can be responsible for 
the incoming regulatory signais from various cellular 
pathways and the direction of the transcriptional activity 
of a target gene [26], The function of bHLH48 is not 
available, although this TF was characterized in a number of 
plant species including Arabidopsis (www.ncbi.nlm.nih.gov/ 
gene/8 18831 ), rice (www.ncbi.nlm.nih.gov/gene/ 
4341525), Glycine max (soybean) (www.ncbi.nlm.nih.gov/ 
gene/100814124), corn ( Zea mays ) (www.maizegdb.org/ 
gene_center/gene/bhlh48) and Jatropha (Jatropha curcas ) 
(www.ncbi.nlm.nih.gov/gene/1 05634733). In Arabidopsis, 
bHLH74 was reported to regulate the expression of a subset 
of genes involved in cell expansion, while bHLH78 was 
reported to trigger flowering in response to blue light 
[27]. bHLH96 and bHLH155 hâve no known rôle, except 
that the expression of the first one was reported to be 
constitutive in roots, leaves, stems and flowers, while the 
second is involved in root development [24], The bHLH106 
was recently reported to parti cipate in sait stress tolérance 
in Arabidopsis [28], Knockout lines of this TF were more 
sensitive to NaCl, KCl, LiCl, ABA, and to low températures 
than the wild type. 

The third type of TF, namely WRKY, is a member of a 
family of TFs functioning in plant responses to various 
physiological processes of végétative and reproductive 
growth [29], This type of TF is important in plant defense 
mechanisms, including MAMP-triggered immunity (MTI) 
or PAMP-triggered immunity (PTI), effector-triggered 
immunity (ETI), and systemic acquired résistance [30], It 
was also proven to hâve a regulatory rôle in sait and 
mannitol stresses in Arabidopsis [31] and in cold and 
drought responses in barley ( Hordeum vulgare ) [32]. In 
addition, WRKY TFs regulate trichome and seed coat 
development in Arabidopsis [33], sesquiterpene biosyn- 
thesis in cotton ( Gossypium hirsutum ) [34], seed develop- 
ment in barley [35] and Arabidopsis [36], respectively, and 
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Cluster 4 


3 



Tissue type 

-•-KJ634217 -«-bHLH48 

Cluster 7 


3 



Tissue type 

-•-KJ 634224 -«-KJ634222 -«-WRKY40 


Cluster 8 


2 



Tissue type 

-•-KJ634231 -•-erf034 


Cluster 9 



Tissue type 


-«-KJ634232-«-WRKY14 -«-bHLH96 -»-WRKY70 -«-WRKY72 
-«-bHLH78 -«-WRKY9 -^erf5 -»-erf bbm2 -«-bHLH106 

Cluster 12 


3 



Tissue type 

-•-KJ634221 -*-bHLH155 -«-WRKY -«-WRKY50 -*-bHLH74 -«-WRKY11 


Fig. 1. Co-expression of uspA-Uke transcripts and putative transcription factors (see Table 1 for details) of erf, bHLH and WRKY superfamilies in clusters 4, 7, 
8, 9 and 12 (see Table S2) across Catharanthus roseus plant tissues. iML: immature leaf, ML: mature leaf, FL: flower, R: root, ST: stem, S: seedling, SC: 
suspension culture, HR: hairy root. TF1 -4: transcription factors with erf domain, TF5-1 0: transcription factors with bHLH domain, TF1 1-17: transcription 
factors with WRKY domain. 
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actin 

iML ML FL R 
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KJ634217 

iML ML FL R 



KJ634222 

iML ML FL R 



KJ634231 

iML ML FL R 



bHLH48 

iML ML FL R 



KJ634224 WRKY40 

iML ML FL R iML ML FL R 



erft)34 

iML ML FL R 



Fig. 2. Semi-quantitative RT-PCR used for validating RNA-Seq analysis for three transcription factors (bHLH48, WRKY40 and erf034) and four uspA (IDs 
KJ634217, KJ634222, KJ634224and I<J634231)genes of clusters 4, 7 and 8. Complété description of different TFs and co-expressing uspA-like genes and their 
expression patterns in immature leaf (iML), mature leaf (ML), flower (FL) and root (R) are shown in Table S2 and Fig. l.The actin gene was used as the house- 
keeping control. Amplicons were run against 100-bp ladder DNA standard (Bioron). 


senescence in Arabidopsis [37]. The WRKYl of tobacco 
[38], similar to WRKY9 in Arabidopsis [39], was reported 
to be phosphorylated by the MAPK, namely SIPK, in order 
to médiate HR-like cell death. The WRKYl 1 and WRKY70 
were reported more recently to participate in regulating 
the signal transduction pathways involved in triggering 
induced systemic résistance in Arabidopsis [40], The 
WRKYl 4 was reported to act as a key regulatory rôle in 
methanol-induced tryptophan (Trp) and Trp-derived 
secondary métabolite biosynthesis during senescence 
in rice leaves [41 ]. The WRKY40 participâtes as tran- 
scriptional regulator in ABA signal transduction in 
Arabidopsis [42], In the presence of the signal, WRKY40 
relocalises to the nucleoplasm. Interestingly, the reloca- 
lisation was not observed in response to other abiotic or 
biotic stimuli. Recent reports indicated that WRKY50 is 
involved in low oleic acid- and SA-derived repression of 
jasmonic acid (JA)-inducible defense responses in Arabi- 
dopsis [43], The WRKY72 was reported to participate in 
basal immunity in tomato and Arabidopsis, while 
appeared to be largely non-responsive to the defense 
hormone salicylic acid (SA) [44], 

In general, we anticipate that the information in the 
présent study highlights the possible functions of uspA- 
like genes in C. roseus in terms of cell response to 
biotic and abiotic stress. This is because a number of 
C. roseus uspA-like transcripts co-expressed with many 
putative ethylene-responsive, bHLH or WRKY putative 
TFs, whose analogues of some of them are known for 
their action in the plant defense against biotic and 
abiotic stresses. 

Acknowledgments 

This project was funded by the Deanship of Scientific 
Research (DSR), King Abdulaziz University, Jeddah, under 
Grant No. (15-3-1432/HiCi). The authors, therefore, 
acknowledge with thanks technical and financial support 
from DSR. 


Appendix A. Supplementary data 


Supplementary data associated with this article can be 
found, in the online version, at http://dx.doi.Org/10.1016/j. 
crvi.201 6.1 0.004. 


References 

[1 ] S.K. Rijhwani, J.V. Shanks, Effect of elicitor dosage and exposure time on 
biosynthesis of indole alkaloids by Catharanthus roseus hairy root 
cultures, Biotech. Progress 14 (1998) 442-449. 

[2] H. Rischer, M. Oresic, T. Seppànen-Laakso, M. Katajamaa, F. Lammertyn, 
et al., Gene-to-metabolite networks for terpenoid indole alkaloid bio- 
synthesis in Catharanthus roseus cells, Proc. Natl. Acad. Sci. USA 103 
(2006) 5614-5619. 

[3] T. Nystrom, F.C. Neidhardt, Expression and rôle of the universal stress 
protein, UspA, of Escherichia coli during growth arrest, Mol. Microbiol. 
11 (1994) 537-544. 

[4] R. Van der Heijden, D.I. Jacobs, W. Snoeijer, D. Hallard, R. Verpoorte, The 
Catharanthus alkaloids: pharmacognosy and biotechnology, Curr. Med. 
Chem. 11 (2004) 1241-1253. 

[5] W.T. Li, Y.M. Wei, J. R. Wang, C.J. Liu, X.J. Lan, et al., Identification, 
localization, and characterization of putative USP genes in barley, 
Theor. Applied Genet. 121 (2010) 907-917. 

[6] A. Van Moerkercke, M. Fabris, J. Pollier, G.J.E. Baart, S. Rombauts, et al., 
CathaCyc, a metabolic pathway database built from Catharanthus roseus 
RNA-Seq data, Plant Cell Physiol. 54 (2013) 673-685. 

[7] A. Bahieldin, A. Atef, A.M. Shokry, S. Al-Karim, S.G. Al Attas, et al., 
Structural identification of putative USPs in Catharanthus roseus, C. R. 
Biol. 338 (2015) 643-649. 

[8] D.S. Latchman, Transcription factors: an overview, Int. J. Biochem. Cell 
Biol. 29 (1997) 1305-1312. 

[9] S. Lindemose, C. O’Shea, M.K. Jensen, K. Skriver, Structure, function and 
networks of transcription factors involved in abiotic stress responses, 
Int. J. Mol. Sci. 14 (2013) 5842-5878. 

[10] K. Hilger-Eversheim, M. Moser, H. Schorle, R. Buettner, Regulatory rôles 
of AP-2 transcription factors in vertebrate development, apoptosis and 
cell-cycle control, Gene 260 (2000) 1-12. 

[11] D. Eckert, S. Buhl, S. Weber, R. Jàger, H. Schorle, The AP-2 family of 
transcription factors, Genome Biol. 6 (2005) 246. 

[12] Q. Li, R.H. Dashwood, Activator protein 2alpha associâtes with adeno- 
matous polyposis coli/beta-catenin and inhibits beta-catenin/T-cell 
factor transcriptional activity in colorectal cancer cells, J. Biol. Chem. 
279 (2004) 45669-45675. 

[13] T. Nakano, K. Suzuki, T. Fujimura, H. Shinshi, Genome-wide analysis of 
the ERF gene family in Arabidopsis and rice, Plant Physiol. 140 (2006) 
411-432. 








6 


A. Bahieldin et al./C. R. Biologies 340 (2017) 1-6 


[14] D. Hao, M. Ohme-Takagi, A. Sarai, Unique mode of GCC box récogni- 
tion by the DNA-binding domain of ethylene-responsive element- 
binding factor (ERF domain) in plant, J. Biol. Chem. 273 (1998) 
26857-26861. 

[15] P.J. O’Donnel, C. Cal vert, R. Atzorn, C. Wasternack, H.M.O. Leyser, et al., 
Ethylene as a signal mediating the wound response of tomato plants, 
Science 274 (1996) 1914-1917. 

[16] I.A. Penninckx, K. Eggermont, F.R. Terras, B.P. Thomma, G.W. De Sam- 
blanx, et al., Pathogen-induced systemic activation of a plant defensin 
gene in Arabidopsis follows a salicylic acid-independent pathway, 
Plant Cell 8 (1996) 2309-2323. 

[17] O. Lorenzo, R. Piqueras, J.J. Sânchez-Serrano, R. Solano, Ethylene re- 
sponse factorl intégrâtes signais from ethylene and jasmonate path- 
ways in plant defense, Plant Cell 15 (2003) 165-178. 

[18] Y.R. Cao, S.Y. Chen, J.S. Zhang, Ethylene signaling régulâtes sait stress 
response: an overview, Plant Signal. Behav. 3 (2008) 761-763. 

[19] J. Cela, C. Chang, S. Munné-Bosch, Accumulation of 7- rather than 
a-tocopherol alters ethylene signaling gene expression in the 
vte4 mutant of Arabidopsis thaliana, Plant Cell Physiol. 52 (201 1 ) 
1389-1400. 

[20] T. Koyama, H. Nii, N. Mitsuda, M. Ohta, S. Kitajima, et al., A regulatory 
cascade involving class II ethylene response factor transcriptional 
repressors opérâtes in the progression of leaf senescence, Plant Physiol. 
162 (2013) 991-1005. 

[21] G.H. Son, J. Wan, H.J. Kim, X.C. Nguyen, W.S. Chung, et al., Ethylene- 
responsive element-binding factor 5, ERF5, is involved in chitin-in- 
duced innate immunity response, Mol. Plant Microbe Interact. 25 
(2012) 48-60. 

[22] H. Wang, Z. Huang, Q. Chen, Z. Zhang, H. Zhang, et al., Ectopic 
overexpression of tomato JERF3 in tobacco activâtes downstream 
gene expression and enhances sait tolérance, Plant Mol. Biol. 55 
(2004) 183-192. 

[23] C. Weiste, T. Iven, U. Fischer, L. Onate-Sânchez, W. Drôge-Laser, In 
planta ORFeome analysis by large-scale over-expression of GATEWAY 
compatible cDNA clones: screening of ERF transcription factors in- 
volved in abiotic stress defense, Plant J. 52 (2007) 382-390. 

[24] M.A. Heim, M. Jakoby, M. Werber, C. Martin, B. Weisshaar, et al., The 
basic helix-loop-helix transcription factor family in plants: a genome- 
wide study of protein structure and functional diversity, Mol. Biol. Evol. 
20 (2003) 735-747. 

[25] A.W. Moore, S. Barbel, L.Y. Jan, Y.N. Jan, A genomewide survey of basic 
helix-loop-helix factors in Drosophila, Proc. Natl. Acad. Sci. USA 97 
(2000) 10436-10441. 

[26] G. Toledo-Ortiz, E. Huq, P.H. Quail, The Arabidopsis basic/helix-loop- 
helix transcription factor family, Plant Cell 15 (2003) 1749-1770. 

[27] Y. Liu, X. Li, K. Li, H. Liu, C. Lin, Multiple bHLH proteins form hetero- 
dimers to médiate CRY2-dependent régulation of flowering-time in 
Arabidopsis, Plos Genet. 9 (2013) el 003861. 

[28] A. Ahmad, Y. Niwa, S. Goto, T. Ogawa, M. Shimizu, et al., bHLH106 
intégrâtes functions of multiple genes through their G-Box to confer 
sait tolérance on Arabidopsis, Plos One 10 (2015) eOl 26872. 

[29] P.J. Rushton, I.E. Somssich, P. Ringler, Q.J. Shen, WRKY transcription 
factors, Trends Plant Sci. 15 (2010) 247-258. 


[30] N. Ishihama, H. Yoshioka, Post-translational régulation of WRKY tran- 
scription factors in plant immunity, Curr. Opin. Plant Biol. 15 (2012) 
431-437. 

[31] L. Chen, Y. Song, S. Li, L. Zhang, C. Zou, et al., The rôle of WRKY 
transcription factors in plant abiotic stresses, Biochim. Biophys. Acta 
1819 (2012) 120-128. 

[32] C. Mare, E. Mazzucotelli, C. Crosatti, E. Francia, A.M. Stanca, et al., Hv- 
WRKY38: a new transcription factor involved in cold- and drought- 
response in barley, Plant Mol. Biol. 55 (2004) 399-416. 

[33] C.S. Johnson, B. Kolevski, D.R. Smyth, TRANSPARENT TESTA GLABRA2, a 
trichome and seed coat development gene of Arabidopsis, encodes a 
WRKY transcription factor, Plant Cell 14 (2002) 1359-1375. 

[34] Y.H. Xu, J.W. Wang, S. Wang, J.Y. Wang, X.Y. Chen, Characterization of 
GaWRKYl, a cotton transcription factor that régulâtes the sesquiter- 
pene synthase gene (+)-delta-cadinene synthase-A, Plant Physiol. 135 
(2004) 507-515. 

[35] M. Luo, E.S. Dennis, F. Berger, W.J. Peacock, A. Chaudhury, MIN1SEED3 
(MINI3), a WRKY family gene, and HAIKU2 (1KU2), a leucine-rich repeat 
(LRR) KINASE gene, are regulators of seed size in Arabidopsis, Proc. Natl. 
Acad. Sci. 102 (2005) 17531-17536. 

[36] C. Sun, S. Palmqvist, H. Olsson, M. Boren, S. Ahlandsberg, et al., A novel 
WRKY transcription factor, SUSIBA2, participâtes in sugar signaling in 
barley by binding to the sugar-responsive éléments of the isol pro- 
moter, Plant Cell 15 (2003) 2076-2092. 

[37] Y. Miao, A. Smykowski, U. Zentgraf, A novel upstream regulator of 
WRKY53 transcription during leaf senescence in Arabidopsis thaliana, 
Plant Biol. (Stuttgart, Germany) Suppl. 1 (2008) 110-120. 

[38] F.L. Menke, H.G. Kang, Z. Chen, J.M. Parle, D. Kumar, et al., Tobacco 
transcription factor WRKY1 is phosphorylated by the MAP kinase SIPK 
and médiates HR-like cell death in tobacco, Mol. Plant Microbe Interact. 
18 (2005) 1027-1034. 

[39] H. Adachi, T. Nakano, N. Miyagawa, N. Ishihama, M. Yoshioka, et al., 
WRKY transcription factors phosphorylated by MAPK regulate a plant 
immune NADPH oxidase in Nicotiana benthamiana, Plant Cell 27 (2015) 
2645-2663. 

[40] C.H. Jiang, Z. Y. Huang, P. Xie, C. Gu, K. Li, et al., Transcription factors 
WRKY70 and WRKY11 served as regulators in rhizobacterium Bacillus 
cereus AR156-induced systemic résistance to Pseudomonas syringae pv. 
tomato DC3000 in Arabidopsis, J. Exp. Bot. 67 (2016) 157-174. 

[41] K. Kang, S. Parle, U. Natsagdorj, Y.S. Kim, K. Bâcle, Methanol is an 
endogenous elicitor molécule for the synthesis of tryptophan and 
tryptophan-derived secondary métabolites upon senescence of de- 
tached rice leaves, Plant J. 66 (2011) 247-257. 

[42] K. Geilen, M. Bôhmer, Dynamic subnuclear relocalisation of WRKY40 in 
response to Abscisic acid in Arabidopsis thaliana, Sci. Rep. 5 (2015) 
(Article 13369). 

[43] Q.M. Gao, S. Venugopal, D. Navarre, A. Kachroo, Low oleic acid-derived 
repression of jasmonic acid-inducible defense responses requires the 
WRKY50 and WRKY51 proteins, Plant Physiol. 155 (2011) 464-476. 

[44] K.K. Bhattarai, H.S. Atamian, I. Kaloshian, T. Eulgem, WRKY72-type 
transcription factors contribute to basal immunity in tomato and 
Arabidopsis as well as gene-for-gene résistance mediated by the to- 
mato R gene Mi-1, Plant J. 63 (2010) 229-240. 


C. R. Biologies 340 (2017) 7-12 



ELSEVIER 


Contents lists available at ScienceDirect 

Comptes Rendus Biologies 

www.sciencedirect.com 



Molecular biology and genetics/Biologie et génétique moléculaires 


Expression of the Galanthus nivalis agglutinin ( GNA ) gene in 
transgenic potato plants confers résistance to aphids 



CrossMark 


Xiaoxiao Mi a ’ b ’ c , Xue Liu a ’ b , Haolu Yan a ’ b , Lina Liang a,b , Xiangyan Zhou a,b , 
Jiangwei Yang a,b , Huaijun Si a ’ b , Ning Zhang a ' b ’* 

11 Gansu Provincial Key Laboratoiy ofAridland Crop Science, Gansu Key Laboratory ofCrop Genetic and Germplasm Enhancement, Gansu 
Agricultural University, Lanzhou 730070, People's Republic of China 

b College ofLife Science and Technology, Gansu Agricultural University, Lanzhou 730070, People's Republic of China 
c The Eighth Middle School ofTianshui City, Tianshui 620500, Gansu Province, People's Republic of China 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 20 July 2016 

Accepted after révision 23 October 2016 

Available online 7 December 2016 


Keywords: 

Aphid 

GNA 

Transgenic potato 

Insect résistance 

Leaf and stem-specific promoter 


Aphids, the largest group of sap-sucking pests, cause significant yield losses in agricultural 
crops worldwide every year. The massive use of pesticides to combat this pest causes 
severe damage to the environment, putting in risk the human health. In this study, 
transgenic potato plants expressing Galanthus nivalis agglutinin (GNA) gene were 
developed using CaMV 35S and ST-LS1 promoters generating six transgenic Iines 
(35S1-35S3 and ST1-ST3 corresponding to the first and second promoter, respectively). 
Quantitative real-time polymerase chain reaction (qRT-PCR) analysis indicated that the 
GNA gene was expressed in leaves, stems and roots of transgenic plants under the control 
of the CaMV 35S promoter, while it was only expressed in leaves and stems under the 
control of the ST-LS1 promoter. The levels of aphid mortality after 5 days of the inoculation 
in the assessed transgenic iines ranged from 20 to 53.3%. The range of the aphid population 
in transgenic plants 15 days after inoculation was between 17.0 ±1.43 (ST2) and 
36.6 ± 0.99 (35S3) aphids per plant, which corresponds to 24.9-53.5% of the aphid population 
in non-transformed plants. The results of our study suggest that GNA expressed in transgenic 
potato plants confers a potential tolérance to aphid attack, which appears to be an alternative 
against the use of pesticides in the future. 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 


1. Introduction 

Aphids are considered as the largest group of sap- 
sucking pests that cause significant yield losses in 
agricultural crops worldwide [ 1 — 3 ]. They induce damage 
to their host plants, modifying plant metabolism, ingesting 
plant nutrients from the phloem, and vectoring plant- 
pathogenic viruses [4,5]. Macrosiphum euphorbiae (Thom- 


* Corresponding author. Gansu Provincial Key Laboratory of Aridland 
Crop Science, Gansu Key Laboratory of Crop Genetic and Germplasm 
Enhancement, Gansu Agricultural University, Lanzhou 730070, People's 
Republic of China. 

E-mail address : ningzh@gsau.edu.cn (N. Zhang). 


as), and Myzus persicae (Sulzer) aphid species constitute 
one of the major potato pests affecting the production of 
this crop [6], which is considered the fourth most 
important worldwide [7], 

Despite the improvement of biopesticides like toxins 
derived from entomopathogenic fungus [8] and bacteria 
[9], and spiders venom [10] that hâve specificity for target 
pest species [11] and do not affect mammals [12], these 
substances remain in the epidermis, being effective only 
for crewing insects. Notwithstanding the fact they combat 
aphids, which are able to penetrate their long acupuncture 
mouthparts into plant phloem, it is necessary to search for 
effective substances that must be able to reach these active 
sites with a long-time permanence and low digestion ratio 
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once it has been ingested by the insect. One candidate 
substance that compiles with the aforementioned require- 
ment is the snowdrop lectin Galanthus nivalis agglutinin 
(CNA), which is able to cross the rnidgut epithelium [13] 
remaining stable and active within the insect gut after 
having been ingested. It has been reported that GNA 
confers résistance to crewing and sap-sucking insects in 
rice, tobacco, cotton, râpe, and wolfberry [14-16], without 
toxicity to higher animais [17], The use of genetic 
engineering that allows the synthesis of GNA to improve 
the résistance of important crops like potato to aphids is 
yet a pending issue that needs to be implemented. 

Transgenic technology has allowed the expression of a 
broad spectrum of promoters, which potentially are 
involved in the synthesis of lectins, one of them, the 
CaMV 35S promoter, is most widely used in transgenic 
plants, since it not only affects the associated transgene, 
but also exerts influence in thousands of base pairs up- or 
downstream of the insertion site on a given chromosome 
[18,19], However, the foreign gene, driven by CaMV 35S 
promoter, is expressed in ail tissues during plant growth 
and development [20j. It causes consumption of excessive 
matter and energy within the cells for the expression of the 
target gene, affecting its temporal and spatial effectiveness 
[21 ]. Furthermore, there are other promoters that are only 
expressed in photosynthetic tissues; this is the case of ST- 
LS1 , a light-inducible promoter whose expression has been 
detected in leaves and photosynthetic stems of potato [22— 
24], In the présent study, the transgenic potato plants 
expressing GNA gene driven by CaMV 35S and ST-LS1 
promoters respectively were obtained to find the différ- 
ence of tissue-specific expression and to estimate the 
résistance of the transgenic potato plants to aphids. Our 
aims were to obtain the transgenic potato plants express- 
ing the GNA gene driven by CaMV 35S and ST-LS1 
promoters, respectively, to find the différence of the 
tissue-specific expression of the GNA gene and further to 
improve the résistance of potato to aphids. 

2. Materials and methods 

2.3. Plants and insects 

The potato cv. ‘Atlantic’ was propagated in vitro by sub- 
culturing single-node cuttings on Murashige and Skoog 
medium [25] supplemented with 3% sucrose and 0.45% 
agar. Plantlets were grown in 150-mL flasks under white 
fluorescent light during 1 6 h and in the dark during 8 h, at a 
température of 24 ± 2 °C. Microtubers were induced in the 
dark at 24 ± 2 C in an MS medium supplemented with 8% 
sucrose and 0.45% agar [26], Green peach aphids (M. persicae, 
Sulzer) were collected from Yuzhong County in Gansu 
Province, China. Aphids were reared in a light incubator at 
25 ± 2 °C under a 14 h light/10 h dark photoperiod, and fed 
with plantlets of potato cv. ‘Atlantic’. 

2.2. Construction of the plant expression vector 

The fragment of GNA gene with nucléotide sequence 
(GenBank accession No. M55556.1) [27] was digested with 
BamH I and Sal I from the clone vector and ligated into a 


binary vector pBI121 [28] and pBI121-ST-LSl [24] 
resulting in the recombinant vectors pBI121-CaMV35S- 
GNA and pBI121-ST-LSl-GNA. These vectors were trans- 
formed into Escherichia coli DH5a respectively, and further 
verified by the same restriction endonuclease digestion. 
After that, pBIl 21 -CaMV35S-GNA and pBIl 21 -ST-LS1 -GNA 
were transformed into Agrobacterium tumefaciens 
LBA4404, respectively, using the freeze-thaw method [29], 

2.3. Transformation of potato 

Potato transformation was performed according to the 
protocol of Si et al. [30], Microtuber slices of potato cv. 
‘Atlantic’ were co-cultured (media-MS + 1 mg/L IAA + 0.2 mg/ 
L GA 3 + 0.5 mg/L 6-BA + 2 mg/L ZT) for 2 days with 
A. tumefaciens LBA4404 containing the plasmid pBI121- 
CaMV35S-GNA and pBI121-ST-LSl-GNA, respectively, then 
transferred into a sélection media supplemented with 50 mg/ 
L kanamycin. When green buds sprouted from the surface of 
the slices and reached a length of 1 cm, they were transferred 
to a sélective rooting medium containing 100 mg/L of 
kanamycin and 200 mg/L of carbenicillin. Plantlets with 
well-developed roots were propagated for further molecular 
analysis. 

High-quality DNA was isolated from the leaves of 
putatively transformed and non-transformed (NT) control 
potato plants for PCR according to the method proposed by 
Edwards et al. [31 ]. The part of the coding sequence of the 
GNA gene was amplified using a PCR Screening Kit 
(GenStar, Beijing, China) with forward primer (5'- 
GCGGATCCATGGCTAAGGCAAGTCTCC-3') and reverse 
primer (5'-GTACGAGCTCTTACnTGCCGTCACAAGCT-3'). 
Amplification was performed in a thermal cycler 
(Tl 00™, BIO-RAD) programmed for one cycle of 3 min 
at 94 C followed by 30 cycles of 30 s at 94 C, 30 s at 66 °C, 
and 1 min at 72 °C. A final extension step was performed 
for 5 min at 72 °C. The amplification products with 500 bp 
in length were separated by electrophoresis on 1.0% 
agarose gels treated with GoldView II staining. 

2.4. Gene expression analysis by qRT-PCR 

Total RNAs were isolated from the transgenic fines and 
NT control using RNAsimple Total RNA Kit (lot#N2822, 
TIANGEN, Beijing, China) following the manufacturées 
instructions. Reverse transcription was performed in 20 p,L 
reaction mixture with the RevertAid First Strand cDNA 
Synthesis Kit (Cat No: 3K1622, Thermo Scientific) and qRT- 
PCR amplification was performed in 20 p,L of the reaction 
mixture with the SuperReal PreMix Plus (SYBR Green) 
(lot#N3113, TIANGEN, Beijing, China), 10p,M of each 
primer (e/3a as an internai control gene and forward and 
reverse primers: 5'-CAAGGATGACCCAGCCAAG-3' and 5'- 
TTCCTT ACCTGAACGCCTGT-3', and the gene-specific for- 
ward and reverse primers: 5'-CTCACCACTTACGCACAAGC- 
3' and 5'-CGGCAATATCCTCnTCTCG-3'). Reactions were 
conducted with an ABI3000 device (Applied Biosystems 
3000 Real-Time PCR System) using the default cycling 
conditions (30 s at 95 °C and 40 cycles of 95 °C for 5 s, 60 °C 
for 34 s, 1 5 s at 95 °C, 1 min at 60 °C and 95 °C for 15 s). Each 
experiment was repeated three finies independently. After 
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Fig. 1. Insect bioassay setup was used to identify the résistance to aphids of the transgenic potato plants and NT. Left: insect bioassay setup. Right: enlarged 
view of the being tested potato plants. Inset in A shows aphid culture on potato plants. The bar represents 5 mm. 


each reaction, dissociation curve analysis was carried out 
to verify the specificity of the amplification and 2~ AAct was 
used to calculate the relative expression levels. 

2.5. Evaluation of résistance ofthe transgenic potato plants to 
aphids 

Insect résistance bioassays were performed according to 
the method described by Cooper et al. [32] and Du et al. 
[33]. Thus, a 75-mm-long sample from the top of the stem 
was detached from 25-day-grown transgenic and NT plants. 
Stem wounds were wrapped in a wet tampon and each 
plant was placed into a glass cage (65 mm in diameter, 
90 mm in height) (Fig. IA). The plants were inoculated with 
two aphid larvae in the first-instar larva State. Five plants 
were used per transgenic individual and control plant. 
Three réplications were performed for each individual line. 


Ail the cages were placed in the illumination incubator 
(LRIT-300-G) at 24 ± 1 °C. The survival of insects within the 
cages was monitored at an interval 24 h for 15 days (Fig. IB). 
The mean survival and mortality per plant was calculated as 
the numbers of surviving and dead (in the four days following 
inoculation) aphids per day. ANOVA followed by Duncan’s 
multiple-range tests were conducted to compare différences 
in mortality and survival rate among transgenic and NT 
control plants for the insect bioassay experiment. 

3. Results 

3.1. Potato transformation and GNAgene expression analysis 
by qRT-PCR 

After three weeks of infection with A. tumefaciens 
LBA4404 containing the plasmid pBI121-CaMV35S-GNA 


T2 T3 



Fig. 2. Potato transformation. A. Shoot formation directly from transgenic microtuber dises of potato cv. ‘Atlantic’ after three weeks of culture in a sélective 
medium (MS + 1 mg/L IAA + 0.2 mg/L GA3 + 0.5 mg/L 6-BA + 2 mg/L ZT + 50 mg/L kanamycin + 200 mg/L carbenicillin) and incubated under a photoperiod 
with 16 h light/8 h dark cycles at 24 °C. B. Roots were formed in about 10 days when green shoots were transferred to the sélective rooting medium 
(MS + 100 mg/L kanamycin + 200 mg/L carbenicillin). NT: non-transformed potato plant as a négative control; T1-T3: three transgenic potato plants. 
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Fig. 3. Transgenic potato plants vérification by PCR assay. M: DL2000 
marker(TaKaRa); NT: non-transformed potato plant as négative control; 
1: plasmid pBI121-CaMV35S-GNA as positive control; 2-8: seven 
transgenic potato plants. 


and pBI121-ST-LSl-GNA respectively, green buds sprouted 
from the surface of the microtuber slices (Fig. 2A). Roots 
were formed in about 10 days (Fig. 2B) when the buds were 
transferred to a sélective rooting media. PCR assay 
demonstrated that transformed plants showed a 500 bp 
amplification product, whereas no product was found in 
NT plants (Fig. 3). Here, we selected three transgenic lines 
35S1-35S3 from the CaMV 35S promoter and three 
transgenic lines from the ST-LS1 promoter for further 
analysis. qRT-PCR analysis showed that the CNA gene was 
expressed in roots, stems and leaves of the transgenic 
potato plants transformed with pB!121-CaMV35S-GNA 
containing the constitutive promoter CaMV 35S (Fig. 4). 
Flowever, in the transformation carried out by pBI121-ST- 
LS1-GNA, the CNA gene was expressed in ail the organs 
except roots. CNA gene expression was higher in leaves 
than in stems of the transgenic plants (Fig. 4). 
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Fig. 4. Tissue-specific expression assay of CNA gene in the transgenic 
potato plants by qRT-PCR. A. Transgenic potato plants transformed with 
pB1121-CaMV35S-GNA. 35S1-35S3 are three different transgenic plant 
lines. B. Transgenic potato plants transformed with pBI121-ST-LSl-GNA. 
ST1-ST3 corresponds to different transgenic plant lines. Error bars 
indicate standard déviations obtained from three independent replicates. 


3.2. Bioassay of résistance of the transgenic potato plants to 
aphids 

3.2.1. Aphid mortality analysis 

The aphid mortality levels after 5 days of inoculation in 
the transgenic lines ranged between 20 (35S1) and 53.3% 
(STI ) (Fig. 5). Différences (at P < 0.05) in the mortality level 
against the NT control were found, except for 35S1 and 
35S2 lines. 



Non-transformed and transgenic potato plants 

Fig. 5. Mortality analysis of first-instar larvae in the first five days following the inoculation in non-transformed and transgenic potato plants. NT: non- 
transformed potato plant; 35S1-35S3: transgenic plants driven by the CaMV 35S promoter; ST1-ST3: transgenic plants driven by the ST-LS1 promoter. 
Different small letters indicated a significant différence at P < 0.05 by Duncan’s multiple-range test. Error bars indicate the standard déviations obtained 
from three independent replicates. 
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3.2.2. Aphid survival analysis 

The aphid population on NT plants increased at a steady 
rate and soared up to a maximum of 68.3 ± 0.95 individuals 
per plant (Fig. 6). On transgenic lines, the aphid population 
increased more slowly than NT and rose to a maximum 
ranged frorn 17.0 ± 1.43 (ST2) to 36.6 ± 0.99 (35S3) individ- 
uals per plant, which corresponded from 24.9 to 53.5% of NT 
aphid population (Fig. 6). 

4. Discussion 

Potato is the world’s fourth-largest food crop, following 
maize, wheat, and rice. Because it supplies dietary fibers, 




Fig. 6. Survival analysis of the aphid population during the first 15 days 
after inoculation with first-instar aphid larvae. NT: non-transformed 
potato plant; 35S1-35S3: transgenic plants driven by CaMV 35S 
promoter; ST1-ST3: transgenic plants driven by the ST-LS1 promoter. 
Different small letters indicate a significant différence at P <0.05 by 
Duncan’s multiple-range test. Error bars indicate the standard déviations 
obtained from three independent replicates. 


carbohydrates, high quantifies of proteins, vitamins and 
minerais, potato is usually regarded as a starchy food or a 
vegetable. China is one of the largest potato-producers 
worldwide. Aphids significantly impact agricultural and 
horticultural crops, either by causing direct damage to 
plants through feeding on the phloem, or indirectly by 
acting as vectors for plant-pathogenic viruses [34], There 
are many measures, such as genetic engineering, which 
hâve been applied to reduce the damage caused by aphids 
to crops [10,35-37], Aphids hinder worldwide potato 
production; therefore, to effectively control the damage of 
the aphid has become an important issue in potato 
production. 

Since the snowdrop lectin (GNA) was separated in 1987 
[38], its encoded gene was cloned and applied to plant 
genetic engineering [27], Because of lacking receptors in 
mammals, GNA is relatively safe to the human body. 
Compared with pesticide residue and ecological destruc- 
tion caused by the long-term use of Chemical pesticides, 
GNA is regarded as the safer and more efficient approach to 
reduce the damage for crops against insects. In our 
experiment, the resuit showed that GNA could express 
in the transgenic potato plants and represent high 
mortality to aphid. It was found that the quantity of 
aphids on the transgenic potato plants was far less than on 
non-transgenic plants. 

Similar to other phloem-specifïc promoters (RSs, ubil) 
with insecticidal effects in sap-sucking homopteran [14], 
our study detected insecticidal effects toward the green 
peach aphid JW. persicae (Sulzer) in transgenic potato plants 
with CaMV 35S and ST-LS1 promoters. However, the resuit 
from qRT-PCR assay showed that the expression of the GNA 
gene driven by the ST-LS1 promoter in the transgenic 
plants was higher in leaves than in stems, and no 
expression in roots (Fig. 4), and even the highest mortality 
(STI) and lowest survival (ST2) of aphids were found in 
lines generated with this promoter. 

Because tuber is the edible portion of potato, in 
accordance with our results, we hypothesize a no or less 
expression of GNA driven by ST-LS1 promoter in this organ 
of transgenic potato plants, guaranteeing a more safety to 
its use as food. However, this hypothesis will hâve to be 
tested in the future. Moreover, our study has been done in 
vitro conditions in plantlets, so it is needed to further grow 
plants in the field to assay their résistance to aphids. In 
addition, as potato is not the only plant-hosting aphids in 
China, the protection of other aphid hosts must be taken 
into considération in order to prevent and control the 
damage caused by aphids for potato. 


Acknowledgements 

This work was supported by the Agriculture Science and 
Technology Innovation Project of Gansu Province, China 
(GNCX-2011-49). We would like to thank Dr. Jiahe Wu 
(Institute of Microbiology, Chinese Academy of Sciences, 
Beijing, China) for kindly provided the vector of GNA gene 
and Dr. David A. Ramirez (International Potato Center, 
Lima, Peru) for editing the manuscript and valuable 
discussion. 


12 


X. Mi et al./C. R. Biologies 340 (2017) 7-12 


Référencés 

[1] S.M. Guo, L.G. Kamphuis, L.L. Gao, J.P. Klingler, J. Lichtenzveig, O. 
Edwards, K.B. Singh, Identification of distinct quantitative trait loci 
associated with defence against the closely related aphids Acythosiphon 
pisum and A. kondo in Medicago truncatula, J. Exp. Bot. 63 (2012) 3913- 
3922. 

[2] L.G. Kamphuis, J. Lichtenzveig, K. Peng, S.M. Guo, J.P. Klingler, K.H. 
Siddique, K.B. Singh, Characterization and genetic dissection of résis- 
tance to spotted alfalfa aphid ( Therioaphis trifolii ) in Medicago trunca- 
tula, J. Exp. Bot. 64 (2013) 5157-5172. 

[3] D. Wang, Q. Liu, H.D. Jones, T. Bruce, L. Xia, Comparative transcriptomic 
analyses revealed divergences of two agriculturally important aphid 
species, BMC Genomics 15 (2014) 1023. 

[4] P.J. Moran, G.A. Thompson, Molecular responses to aphid feeding in 
Arabidopsis in relation to plant defence pathways, Plant Physiol. 125 
(2001) 1074-1085. 

[5] C. Couldridge, H.J. Newbury, B.F. Lloyd, J.B.J. Pritchard, Exploring plant 
responses to aphid feeding using a full Arabidopsis microarray reveals a 
small number of genes with significantly altered expression, Bull. 
Entomol. Res. 97 (2007) 523-532. 

[6] Y. Pelletier, D. Michaud, Insect pest control on potato: genetically-based 
control, in: R.M. Duchesne, G. Boiteau (Eds.), Potato insect pest control: 
development of a sustainable approach, Gouvernement du Québec, 
1995, pp. 69-79. 

[7] Food and Agriculture Organization of the United Nations (FAO), World 
crop production statistics, 2016 http://www.faostat.fao.org , . Accessed 
lOth September 2016. 

[8] W. Liu, Y. Xie, J. Dong, J. Xue, Y. Zhang, Y. Lu, J. Wu, Pathogenicity of three 
entomopathogenic fungi to Matsucoccus matsumurae, PLoS One 9 
(2014) el 03350. 

[9] F. Wang, Y. Liu, F. Zhang, L. Chai, L. Ruan, D. Peng, M. Sun, Improvement 
of crystal solubility and increasing toxicity against Caenorhabditis 
elegans by asparagine substitution in block 3 of Bacillus thuringiensis 
crystal protein Cry5Ba, Appl. Environ. Microbiol. 78 (2012) 7197-7204. 

[10] J.R. Bloomquist, Mode of action of atracotoxin at central and peripheral 
synapses of insects, Invert. Neurosci. 5 (2003) 45-50. 

[11] P.R. Whitehorn, S.O. Connor, F.L. Wackers, D. Goulson, Neonicotinoid 
pesticide reduces bumble bee colony growth and queen production, 
Science 336 (2012) 351-352. 

[12] S.S. Pineda, D. Wilson, J.S. Mattick, G.F. King, The léthal toxin from 
Australian funnel-web spiders is encoded by an intronless gene, PLoS 
One 7 (2012) e43699. 

[13] A.M. Gatehouse, J.A. Gatehouse, M. Bharathi, J. Spence, K.S. Powell, 
Immunohistochemical and developmental studies to elucidate the 
mechanism of action of the snowdrop lectin on the rice brown plant 
hopper, Nilaparvata lugens (Stal), J. Insect Physiol. 44 (1998) 529-539. 

[14] K.V. Rao, K.S. Rathore, T.K. Hodges, X. Fu, E. Stoger, D. Sudhakar, 
Expression of snowdrop lectin (GNA) in transgenic rice plants confers 
résistance to rice brown plant hopper, Plant J. 15 (1998) 469-477. 

[15] V.A. Hilder, K.S. Powell, A.M.R. Gatehouse, J.A. Gatehouse, L.N. Gate- 
house, Y. Shi, W.J. Peumans, Expression of snowdrop lectin in trans- 
genic tobacco plants results in added protection against aphids, 
Transgenic Res. 4 (1995) 18-25. 

[16] X. Foissac, N.T. Loc, P. Christou, A.M. Gatehouse, J.A. Gatehouse, Résis- 
tance to green leaf hopper ( Nephotettix virescens) and brown plant 
hopper ( Nilaparvata lugens) in transgenic rice expressing snowdrop 
lectin ( Galanthus nivalis agglutinin; GNA), J. Insect Physiol. 46 (2000) 
573-583. 

[17] A.M.R. Gatehouse, K.S. Powell, E.J.M.V. Damme, W.J. Peumans, Insecti- 
cidal properties of plant lectins: their potential in plant protection, in: 
A. Pusztai, S. Bardocz (Eds.), Lectins: Biomédical Perspectives, Taylor & 
Francis, London, 1995, pp. 35-57. 

[18] A. Wilmink, B.C.D. Ven, J.J. Dons, Activity of constitutive promoters in 
various species from the Liliaceae, Plant Mol. Biol. 28 (1995) 949-955. 


[19] C. Potenza, L. Aleman, C. Sengupta-Gopalan, Targeting transgene ex- 
pression in research, agricultural, and environmental applications: 
promoters used in plant transformation, In Vitro Cell. Dev. Bio-Plant 
40 (2004) 1-22. 

[20] J.T. Odell, F. Nagy, N.H. Chua, Identification of DNA sequences required 
for activity of the cauliflower mosaic virus 35S promoter, Nature 313 
(1985) 810-812. 

[21] J.R. Gittins, T.K. Pellny, E.R. Hiles, C. Rosa, S. Biricolti, D.J. James, 
Transgene expression driven by heterologous ribulose-l,5-bisphos- 
phate carboxylase/oxygenase small-subunit gene promoters in the 
végétative tissues of apple ( Malus pumila mill.). Planta 210 (2000) 
232-240. 

[22] J. Stockhaus, J. Schell, L. Willmitzer, Corrélation of the expression of the 
nuclear photosynthetic gene ST-LS1 with the presence of chloroplasts, 
EMBO J 8 (1989) 2445-2451. 

[23] P. Eckes, S. Rosahl, J. Schell, L. Willmitzer, Isolation and characterization 
of a light-inducible, organ-specific gene from potato and analysis of its 
expression after tagging and transfer into tobacco and potato shoots, 
Mol. Gen. Genet. 205 (1986) 14-22. 

[24] Y. Qu, N. Zhang, J. Chang, X. Jin, Y.K. Wen, H.J. Si, D. Wang, Cloning and 
functional analysis of light-inducible, and stem and leaf-specific ex- 
pression promoter ST-LS1 in potato ( Solanum tuberosum L), J. Agri. 
Biotechnol. 21 (2013) 828-837. 

[25] T. Murashige, F. Skoog, A revised medium for rapid growth and bioas- 
says with tobacco tissue cultures, Physiol. Plant. 15 (1962) 473-497. 

[26] N. Zhang, H.J. Si, D. Wang, Cloning of rd29A gene promoter from 
Arabidopsis thaliana and its application in stress-resistance transgenic 
potato, Acta Agro. Sin. 31 (2005) 159-164. 

[27] E.J.V. Damme, H. Kaku, F. Perini, I.J. Goldstein, B. Peeters, F. Yagi, B. 
Decock, W.J. Peumans, Biosynthesis, primary structure and molecular 
cloning of snowdrop ( Galanthus nivalis L.) lectin, Eur. J. Biochem. 202 
(1991) 23-30. 

[28] T. Y. Li, Y.C. Tian, X.F. Qing, Studies on high-efficient insect résistance 
transgenic tobacco, China Sin. (B). 24 (1994) 268-276. 

[29] R. Hofgen, L. Willmitzer, Storage of competent cells for Agrobacterium 
transformation, Nucleic Acid Res. 16 (1988) 9877. 

[30] H.J. Si, C.H. Xie, J. Liu, An efficient protocol for Agrobacterium-mediated 
transformation of microtuber and the introduction of antisense class I 
patatin gene into potato, Acta Agro. Sin. 29 (2003) 801-805. 

[31] K. Edwards, C. Johnstone, C. Thompson, A simple and rapid method for 
the préparation of genomic DNA for PCR analysis, Nucleic Acids Res. 16 
(1991) 1349. 

[32] S.G. Cooper, D.S. Douches, E.J. Grafius, Combining genetic engineering 
and traditional breeding to provide elevated résistance in potatoes to 
Colorado potato beetle, Entomol Exp. Appl. 112 (2004) 37-46. 

[33] J.L. Du, D.G. Wu, T.W. Zhang, X.J. Qian, C.Z. Liu, Effect of UV-B for 
different radiation durations on biological characteristics of two color 
morphs of pea aphid Acyrthosiphon pisum (Harris) offspring, Acta 
Agrestia Sin. 20 (2012) 961-966. 

[34] E.Y. Nakasu, M.G. Edwards, E. Fitches, J.A. Gatehouse, A.M. Gatehouse, 
Transgenic plants expressing œ-ACTX-Hvla and snowdrop lectin 
(GNA) fusion protein show enhanced résistance to aphids, Front Plant 
Sci. 5 (2014) 673. 

[35] C.J. Burns, L.J. Mclntosh, P.J. Mink, A.M. Jurek, A.A. Li, Pesticide exposure 
and neurodevelopmental outcomes: review of the épidémiologie and 
animal studies, J. Toxicol. Environ. Health B Crit. Rev. 16 (2013) 127- 
283. 

[36] T. Blacquière, G. Smagghe, C.A.V. Gestel, V. Mommaerts, Neonicotinoids 
in bees: a review on concentrations, side-effects and risk assessment, 
Ecotoxicology 21 (2012) 973-992. 

[37] K. Derecka, M.J. Blythe, S. Malla, D.P. Genereux, A. Guffanti, P. Pavan, 
D.A. Barrett, Transient exposure to low levels of insecticide affects 
metabolic networks of honeybee larvae, PLoS One. 8 (2013) e68191. 

[38] E.J.V. Damme, A.K. Allen, W.J. Peumans, Isolation and characterization 
of a lectin with exclusive specificity towards mannose from snowdrop 
( Galanthus nivalis L.) bulbs, FEBS Lett. 215 (1987) 140-144. 


C. R. Biologies 340 (2017) 13-17 



ELSEVIER 


Contents lists available at ScienceDirect 

Comptes Rendus Biologies 

www.sciencedirect.com 



Molecular biology and genetics/Biologie et génétique moléculaires 

Age-related changes of metallothionein 1/2 and 
metallothionein 3 expression in rat brain 



CrossMark 


Rosaria Scudiero *, Luisa Cigliano, Mariailaria Verderame 

Department of Biology, University of Naples Federico II, via Mezzocannone 8, 80134 Napoli, Italy 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 13 October 2016 

Accepted after révision 15 November 2016 

Available online 8 December 2016 


Keywords: 

Brain aging 

Metallothionein isoforms 
MT1/2 and MT3 gene expression 
Rattus 

Real-time PCR 


Neurodegeneration is one of the main physiological conséquences of aging on brain. 
Metallothioneins (MTs), low molecular weight, cysteine-rich proteins that bind heavy- 
metal ions and oxygen-free radicals, are commonly expressed in various tissues of 
mammals. MTs are involved in the régulation of cell prolifération and protection, and may 
be engaged in aging. Expression of the ubiquitous MTs (1 and 2) and the brain spécifie MT3 
hâve been studied in many neurodegenerative disorders. The research results indicate that 
MTs may play important, although not yet fully known, rôles in brain diseases; in addition, 
data Iack the ability to identify the MT isoforms functionally involved. The aim of this 
study was to analyse the level of gene expression of selected MT isoforms during brain 
aging. By using real-time PCR analysis, we determined the MT1 /2 and MT3 expression 
profiles in cérébral cortex and hippocampus of adolescent (2 months), adult (4 and 
8 months), and middle-aged (16 months) rats. We show that the relative abundance of ail 
types of MT transcripts changes during aging in both hippocampus and cortex; the first 
effect is a generalized decrease in the content of MTs transcripts from 2- to 8-months-old 
rats. After passing middle âge, at 1 6 months, we observe a huge increase in MT3 transcripts 
in both cortical and hippocampal areas, while the MT1/2 mRNA content increases slightly, 
returning to the levels measured in adolescent rats. These findings demonstrate an age- 
related expression of the MT3 gene. A possible link between the increasing amount of MT3 
in brain aging and its different metal-binding behaviour is discussed. 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 


1. Introduction 

Metallothioneins (MTs) are small, intracellular heavy- 
metal-binding proteins involved in métal détoxification 
and protection against oxidative stress [1,2]. Ubiquitously 
présent in living organisms [3,4], the mammalian MT 
family falls into four subgroups: MT1, MT2, MT3, and MT4 
[5[. MT1 and MT2 encode through a sériés of gene 
duplication events for multiple isoforms présent in almost 
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ail tissues, while MT3 and MT4 evolved in a single form to 
accomplish spécifie rôles in brain and epithelium, respec- 
tively [6-8], In the mammalian brain, MT1 and MT2 are 
preferentially expressed in astrocytes and activated 
microglia [9[; MT3 is expressed predominantly in neurons 
and, to a lesser extent, in astrocytes [10], 

It has been demonstrated that MT1 and MT2 protect the 
central nervous System (CNS) in response to experimen- 
tally induced injury [11] and following viral infections [12], 
whereas MT3 neither affects inflammatory responses nor 
plays an important antioxidant rôle [11], Indeed, the MT3 
functional rôle and régulation in CNS is still a matter of 
debate [10,13], MT3 is known as a growth inhibitory factor 
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able to inhibit the survival and neurite formation in 
cultured neurons in vitro [14], MT3 mRNA and protein are 
upregulated following brain injury [15] and downregula- 
ted in Alzheimer’s disease (AD) [14,16], 

The known functions of MTs include metalloregulatory 
rôles in cell growth and différentiation [17], so it is not 
surprising that most mammalian tissues contain age- 
related basal levels of MTs, with the major amount found in 
developing cells [18], Enhanced synthesis of MTs is also 
observed in rapidly proliferating tissues, stressing their 
crucial rôle in normal and neoplastic cell growth [18-20], 
However, the MTs’ rôle in aging is controversial: in 
mammals, MT expression has been found to increase in 
aged kidney [21], whereas decreased significantly with 
aging in the skin [22], 

Aging is a major risk factor for brain neurodegenerative 
disorders, including cerebrovascular disease, AD, Parkin- 
son’s disease (PD) and cancer [23-26], Brain aging is 
accompanied with molecular, functional and genetic 
changes, leading to increased susceptibility to diseases 
and cognitive impairments. As it has been proven in the 
mouse, MTs downregulation results in a progressive 
neurodegeneration, leading to early aging, morbidity, 
and mortality [27]; the neurodegenerative alterations 
are attenuated by MTs overexpression, suggesting the 
neuroprotective rôle of MTs in aging [27], These studies, 
however, did not include the identification of the MT 
isoforms functionally involved. 

In this study, by using real-time PCR analysis, we 
examined and compared the expression of MTi/2 and MT3 
genes during aging in the rat brain. In particular, we 
evaluated the MT1 / 2 and MT3 expression profiles in the 
cérébral cortex and hippocampus of pubertal (2 months), 
adult (4 and 8 months), and middle-aged (16 months) rats. 

Our results demonstrate an age-related expression of 
the MT3 gene, thus, suggesting its implication in physio- 
logical changes associated with aging. 


2. Materials and methods 

2.3. Animal and experimental design 

Male Wistar rats (Charles River, Calco, Como, Italy) of 
2 months of âge (adolescent), with the same starting body 
weight (160 ± 10 g), were individually caged in a tempera- 
ture-controlled room (23 ± 1 °C) submitted to a 12-h light/ 
12-h dark cycle. Animais were housed in the Animal Care 
Facility at the Department of Biology, with ad libitum access 
to water and to a standard diet (Mucedola 4RF21; Settimo 
Milanese, Milan, Italy) up to 4 (social maturity, group 1), 8 
(adulthood, group 2), or 16 (aged, group 3) months. Eight 
animais of each group were anesthetized with chloral 
hydrate (40mg/100g body wt) and killed by décapitation. 
Group 4 was comprised 8 adolescent rats sacrificed four days 
after their arrivai in the animal facility. The brains were 
quickly removed and the cérébral cortex and hippocampus 
were dissected on ice. Samples of each brain région were snap 
frozen in liquid nitrogen immediately and stored at -80 C 
for subséquent RNA isolation. The livers of the rats belonging 
to group 4 were also collected, frozen and stored at -80 "C for 


RNA isolation. The protocols for animal care and use were 
approved by the Scientific Ethics Committee on Animal 
Expérimentation of the University Federico II of Naples. Ail 
experimental animal procedures were carried out in compli- 
ance with the national guidelines for the care and use of 
research animais (D.L.l 16/92, implémentation of EEC direc- 
tive 609/86). Ail efforts were made to minimize animal 
suffering and to reduce the number of animais used. 

2.2. RNA purification and cDNA synthesis 

Total RNA was extracted according to the TRI-Reagent 
(Sigma-Aldrich) protocol. The quality of each total RNA 
was checked by electrophoresis on 2% agarose gel stained 
with ethidium bromide and by measuring the optical 
density at 260/280 nm. A ratio of 1. 8-2.0 was accepted for 
further reverse transcription. The QuantiTect Reverse 
Transcription Kit (Qiagen) was used for the removal of 
genomic DNA contamination and for the subséquent cDNA 
synthesis. Approximately 1 p,g of total RNA was used, 
according to the kit's protocol. 

2.3. Quantitative real-time PCR analysis 

The real-time PCR reactions were carried out in 
quadruplicate in an Applied Biosystems 7500 Real-Time 
System by using the Power SYBR Green Master Mix PCR 
(Applied Biosystems) following procedures recommended 
by the manufacturer. Each SYBR Green reaction (total 
volume; 20 p,L) contained 12 p,L of real-time PCR Master 
Mix, 1 |xL of each of the forward and reverse primers 
(10 p-M), 2 p,L of 1:1 diluted cDNA, and 4 p,L of nuclease- 
free water. Nuclease-free water for the replacement of the 
cDNA template was used as a négative control. For the 
internai standard control, the expression of jd-actin gene 
was quantified. Primer sequences were designed using 
Primer Express software (Applied Biosystems). A single pair 
of spécifie primers for both MT1 and MT2 isoforms was 
designed on the nucléotide sequences of Rattus norvegicus 
MT1 (NM 138826.4) and MT2 (NM 001137564.1). Spécifie 
primers for MT3 were designed on the R. norvegicus 
template NM 053968.3, choosing the most divergent 
sequence tracts rather than MT1/2 isoforms. (3-actin 
primers were designed on the exon junction 75/76 (forward 
primer) on R. norvegicus template NM031 144.2. Ail the 
primers used in real-time PCR analysis are listed in Table 1. 
PCR was performed under the following conditions: holding 
stage of 95 °C per 10 min; cycling stage (45 cycles): 
95“Cxl0s - 60°Cxl0s - 72°Cxl0s; melting stage: 
95 °C x 5 s - 65 °C x 1 min - 95 °C x 30 s - 40 C x 30 s. The 
melting curve analysis of PCR products was performed in 


Table 1 

Spécifie primers used for semiquantitative real-time PCR analyses. 


Name 

Nucléotide sequence (5'— 3') 

Length 

Direction 

MTl /2 F 

ATGGACCCCAACTGCTCCTG 

20-mer 

Sense 

MT1/2_R 

CTTTGCAGACACAGCCCTGGG 

21-mer 

Antisense 

MT3_F 

CCCCTGTCCTACTGGTGGT 

19-mer 

Sense 

MT3_R 

CTGCATTTCTCGGCCTTG 

18-mer 

Antisense 

(3-act_F 

ACCCGCCACCAGTTCGCCAT 

20-mer 

Sense 

(3-act_R 

CGGCCCACGATGGAGGGGAA 

20-mer 

Antisense 
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order to ensure gene-specific amplification. Changes in the 
gene expression relative to the different samples were 
calculated according to the standard 2~ AAct method 
described by Livak and Schmittgen [28], 

2.4. Statistical analysis 

Data are presented as mean ± standard error of the mean 
(SEM) from four separate experiments in each sample. 
Statistical analyses were carried out by StatView software 
(Altura Software, Inc.). The différences between the mean 
values were analysed by one-way analysis of variance 
(ANOVA) followed by Fisher’s LSD test. The différences were 
considered significant when P< 0.05. 

3. Results 

3.1. Expression of MT isoforms in liver and brains of 
adolescent rats 

As a preliminary step, we firstly decided to compare the 
relative abundance of MT1/2 and MT3 transcripts in the 
liver, the cérébral cortex, and the hippocampus of 
adolescent rats. As expected, MTl and MT2 are the major 
forms in the liver, whereas MT3 is the most represented 
isoform (Fig. 1 ) in the two areas of the brain, with the level 
of the cortex higher than that of the hippocampus. 
However, a significant amount of MT3 transcripts was 
also détectable in the liver, as well as lower levels of MTl / 
MT2 isoforms were found both in the cortical and 
hippocampal areas (Fig. 1). 

3.2. Expression of MT isoforms in the aging brain 

Once the MTs gene expression profile in adolescent 
(2-month-old) rats established, we carried out the 
study evaluating MT1/2 and MT3 mRNAs levels in the 


a 
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Fig. 1. Real-time PCR analysis of metallothionein isoforms expression in 
liver, cérébral cortex and hippocampus of 2-month-old rats. MTl /2 and 
MT3 mRNA levels are shown relative to fi-actin mRNA. The data represent 
the mean ± SEM (eight rats per group). 


hippocampus and the cortex of adult (4- and 8-month-old), 
and middle-aged (16-month-old) rats. Real-time PCR 
revealed a huge increase in MT3 expression in both 
régions of the brain, in aged rats. In particular, the level of 
this transcript increased by 16 times compared to the 2- 
month-old rats in the cérébral cortex (Fig. 2a) and by about 
14 times in the hippocampus (Fig. 2b). Interestingly, we 
did not observe a linear increase in MT3 transcripts with 
the âge from 2 to 8 months in both the hippocampus and 
the cortex; on the contrary, a slight but significant decrease 
in MT3 transcripts was observed in adult rats when 
compared with young rats (Fig. 2). 

As regarding MTl /2 expression in brain areas, real-time 
PCR analysis demonstrated that the amount of these 
transcripts decreased in adults (4- and 8-month-old rats), 
and increased slightly only in aged rats (16 months), 



Fig. 2. Effect of aging on MT3 gene expression in rat cortical (a) and hippocampal (b) areas. The MT3 mRNA expression was normalized to that of fi-actin 
mRNA and converted infold change, compared with 2-month-old rats. The data represent the mean ± SEM (eight rats per group). Significance of différences is 
shown. *P< 0.05; "P< 0.001. “Significance vs. 2 months; ““Significance vs. 4 and 8 months. 



Fig. 3. Effect of aging on MT! and MT2 genes expression in rat cortical (a) and hippocampal (b) areas. The MT1 /2 mRNAs expression was normalized to that 
of (3-actin mRNA and converted in fold change, compared with the 2-month-old rats. The data represent the mean ± SEM (eight rats per group). The 
significance of différences is shown. ‘P < 0.05; a Significance vs. 2 months; “Significance vs. 4 and 8 months. 


reaching more or less the level measured in 2-months rats, 
both in cortical (Fig. 3a) and hippocampal (Fig. 3b) areas. 

4. Discussion 

Aging is a dominant risk factor for many CNS disorders, 
among which AD and PD [23,24,26,29], Aging contributes 
to the physiological décliné of cells and tissues; in the 
brain, it is associated with inflammation, increased 
oxidative damage, protein aggregate accumulation and 
demyelination, whereas it gives rise to a modest neuron 
and synapse loss [30,31 ]. Age-related changes in CNS gene 
expression levels hâve been reported and many genes, and 
relative proteins, hâve been suggested as prognostic 
biomarkers of neurodegenerative disorders [32-34], Re- 
cently, MT genes hâve been introduced in the list of early 
biomarkers of neuronal disorders [27,29]. It has been 
demonstrated that MTs overexpression in mouse brain 
provide neuroprotection against deleterious conséquences 
of oxidative and nitrosative stress. However, the authors 
did not identify the MT isoform(s) responsive to neuro- 
degeneration [27]. Therefore, a major aim of this study was 
to show whether changes in MT isoforms ratio occur with 
aging in rats’ hippocampus and cérébral cortex. 

MT3 is reported as the most abundant isoform 
expressed in the brain, and our data are consistent with 
this resuit. Here, we also show that the relative abundance 
of ail types of MT transcripts change during aging in both 
hippocampus and cortex; surprisingly, the first effect is a 
generalized decrease in the content of MTs transcripts 
from 2-month-old to 8-month-old rats. This pattern can be 
explained according to two hypothèses. In the first one, it is 
possible to assume that the rate of MT expression 
measured in adolescent rats is strictly linked to the 
physiological State of the brain, consisting of active, 
proliferating cells with a high metabolic rate. On the other 


hand, it cannot be excluded that the higher content of MTs 
detected in adolescent rats compared with adults is due to 
the stress that these animais undergo in their early stages 
of life. Indeed, once in the animal facility, the rats are 
maintained in individual cages under optimum conditions 
of humidity, température, photoperiod, with water ad 
libitum and right amount of food; these conditions would 
lead to a decrease in MTs expression, lacking the stressful 
factors that may cause their temporary increase. 

After passing middle âge, at 16 months, we observe a 
huge increase in MT3 transcripts in both cortical and 
hippocampal areas, while the MT1/2 mRNA content 
increases slightly, returning to the levels measured in 
young rats (2 months). These findings demonstrate an age- 
related expression of the MT3 gene, whereas the lower 
levels of MT1/2 genes might suggest a marginal rôle of 
these players, if any, in the aging process of the brain. 
Structural and functional studies on MT3 demonstrated 
the stronger Cu-thionein character of MT3 compared to 
that of MT1 and MT2 [35,36], corroborating the hypothe- 
sized rôle of this MT isoform in the brain as a key element 
in the homeostasis and metabolism of copper ions in the 
CNS [36], Recently, Fu et al. hâve ascertained an age- 
dependent increase in the Cu levels in rat brains [37], 

Therefore, it is conceivable that the increasing amount 
of MT3 expression in the aging brain might be linked to its 
different metal-binding behaviour with respect to the 
other MT isoforms and to its spécifie neuronal growth 
inhibitory bioactivity. Together, these features could better 
meet the needs of aged brain cells. 
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Preclinical researches on mesenchymal stem cells (MSCs) transplantation, which is used to 
treat hypoxic-ischemic (HI) brain damage, hâve received inspiring achievements. 
However, the insufficient migration of active cells to damaged tissues has limited their 
potential therapeutic effects. There are some évidences that hypoxia inducible factor-1 
alpha (HIF-la) promotes the viability and migration of the cells. Here, we aim to 
investigate whether overexpression of HIF-la in MSCs could improve the viability and 
migration capacity of cells, and its therapeutic efficiency on HI brain damage. In the study, 
MSCs with HIF-la overexpression was achieved by recombinant lentiviral vector and 
transplanted to the rats subséquent to HI. Our data indicated that overexpression of HIF- 
la promoted the viability and migration of MSCs, HIF-la overexpressed MSCs also had a 
stronger therapeutic efficiency on HI brain damaged treatment by mitigating the injury on 
behavioral and histological changes evoked by HI insults, accompanied with more MSCs 
migrating to cérébral damaged area. This study demonstrated that HIF-1 a overexpression 
could increase the MSCs’ therapeutic efficiency in HI and the promotion of the cells' 
directional migration to cérébral HI area by overexpression may be responsible for it, 
which showed that transplantation of MSCs with HIF-la overexpression is an attractive 
therapeutic option to treat Hl-induced brain injury in the future. 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 


1. Introduction 

Hypoxic-ischemic (HI) brain damage is a clinical 
syndrome with high mortality and morbidity [1], and 
the survivors hâve suffered widely mental and physical 
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disabilities [2], However, there are no effective treatments 
available for HI brain damage at présent. 

Mesenchymal stem cells (MSCs), which are easily 
obtained from many tissues, hâve aroused great interest 
in regenerative medicine because of their multi-directional 
différentiation potential and hypo-immunogenicity prop- 
erties [3], In recent years, MSCs hâve showed a significant 
neuroprotective effects in experimental rnodels of HI brain 
damage [4-6]. However, the main problem is that the low 
number of active cells reaching the target tissues after 
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transplantation blocks its potential therapeutic effects 
[7,8]. 

Previous studies suggested that delivering the desired 
genes into the cells may overcome some limitations of cells 
and facilitate therapeutic efficiency after cell transplanta- 
tion [9,10], Hypoxia inducible factor-1 alpha (HIF-la) is 
one of the most important regulators in response to oxygen 
levels, and the stabilization of HIF-1 a in cells may promote 
activation and migration [11-13], 

In the présent study, we intended to overexpress HIF- 
la in MSCs by recombinant lentivirus and investigate the 
issue whether HIF-la overexpression could promote the 
therapeutic efficiency of MSCs. 


2. Materials and methods 

2.1. Animais 

Ail procedures related to the use of animais complied 
with the guide published by NIH [14]. Adult male Sprague- 
Dawley rats (80-100 g) were purchased from the Animal 
Center of Zhejiang Chinese Medical University, Hangzhou, 
China (Laboratory Animal Certificate: scxk 2008-0016). 

2.2. Génération, culture and pretreatment of MSCs 

MSCs were obtained from rat fémoral and tibia bone 
marrow as previously described [15], Briefly, muscles, 
connective tissue and epiphyses were removed from fémur 
and tibia. Marrow was harvested by 1 .0 mL syringe with 
cell medium that consist of Dulbecco’s modified Eagle’s 
medium/Ham’s F-12 (DMEM/F-12; invitrogen, Grand 
Island, NY, USA), 10% fêtai bovine sérum (FBS; Life 
Technologies, Carlsbad, CA, USA) and 1% penicillin- 
streptomycin (HyClone Laboratories, Inc., South Logan, 
Utah, USA). Cells were then centrifuged for 10 min at 
1000 rpm, re-suspended in the cell medium and cultured 
at 37 °C, 5% C0 2 . After 24 h, ail non-adherent cells were 
removed by medium exchange; the fresh medium was 
subsequently replaced every 2 days. The adhèrent cells 
were trypsinized at 80% confluence, and passaged at 1:2. 
MSCs from passage 3 were used in the subséquent 
experiments. 

MSCs had received a hypoxia pretreatment during 24 h 
to mimic the microenvironment in the damaged area in 
order to activate the endogenous gene expression, which is 
the natural effect of cells in response to hypoxia. Briefly, 
the cells were placed in a hypoxia incubator and balanced 
with a gas mixture (1% 0 2 , 5% C0 2 , 94% N 2 ) for 10 min at a 
flow rate of 1 L/min. Then, the hypoxia incubator was 
sealed to maintain a hypoxie environment before being put 
into a cell incubator under 5% C0 2 at 37 °C. 

2.3. Lentivirus-mediated gene transfer and experimental 
groups of cells 

MSCs with HIF-la overexpression were constructed by 
a lentivirus System following the manufacturées instruc- 
tions (Invitrogen, Carlsbad, CA, USA). The cells were 
infected with a lentivirus encoding either GFP alone or 


both GFP and HIF-la. Then the MSCs were assigned into 
three experiment groups: HIF-la-MSCs (transfection with 
HIF-la and GFP), GFP-MSCs (transfection with GFP), NM- 
MSCs (undisturbed normal MSCs). 

The transfection efficiency of the cells was detected by 
fluorescence microscopy analysis after transfecting for 
48 h. The protein expression level of HIF-la in MSCs was 
detected by Western blotting. 

2.4. Western blotting 

The expression level of HIF-la was evaluated by 
western blotting after hypoxia pretreatment. Equal 
amounts of protein from MSCs were separated by 10% 
SDS-PAGE and transferred to polyvinylidene difluoride 
membranes (Millipore, Boston, MA). After being blocked 
with 5% nonfat dry milk for 2 h, the membranes were 
incubated with rabbit polyclonal antibodies against either 
HIF-la (1:1000 dilution, Abcam Technology, Cambridge, 
MA, USA) or GAPDH (1:2000 dilution, Abcam Technology, 
Cambridge, MA, USA) overnight at 4 °C. The membranes 
were then incubated with a second antibody (anti-rabbit 
IgG HRP, 1 :3000; Abcam Technology, Cambridge, MA, USA) 
and ECL détection Systems (Thermo Fisher Scientific, 
Waltham, MA, USA) were used for détection, GAPDH 
was used as a loading control. 

2.5. MTT assay 

An assay with 3-(4,5-dimethylthiazol-2-yl)-2,5-diphe- 
nyltetrazolium bromide (MTT) was performed to déter- 
mine the influence of HIF-la overexpression on MSC 
viability. Briefly, pretreated MSCs were seeded in 96-well 
plates at a density of 2000 cells per well, and cultured for 
48 h, and NM-MSCs were used as a control. 10 p,L (5 g/L in 
PBS) of MTT (Sigma Chemical Co., St. Louis, MO, USA) were 
added to each well and incubated at 37 °C for 4 h. The 
supernatants were carefully discarded and 150 p,L of 
dimethyl sulfoxide (Sigma Chemical Co., St. Louis, MO, 
USA) was added. Absorbance was measured at 490 nm by 
microplate reader (Bio-RAD Model 680, Cal, USA). 

2.6. CM-Dil labelling 

In order to make it easily observable, MSCs were stained 
with chloromethylbenzamido dialkylcarbocyanine (CM- 
Dil, Molecular Probes, USA) at the dose of 1 p,g CM-Dil per 
million cells, and washed twice after incubation at 37 °C for 
20 min. Labeled MSCs were used in the following 
experiments. 

2.7. Migration assay 

A cell migration assay was performed in Transwell 
chambers (pore size: 8 p,m, polycarbonate membrane; 
Corning Incorporated, NY, USA) to evaluate the migratory 
capacity of MSCs. Transwell chambers were put into six 
well plates. CM-Dil labeled cells were harvested in a 
serum-free medium; the density was adjusted to 
2 x 10 5 cells/ml. 200 p,L of cell suspension were plated 
in the upper chamber and 500 p,L of a fresh medium 


20 


D. Un et al./C. R. Biologies 340 (2017) 18-24 


containing 10% FBS were added to the lower chamber. 
After incubation at 37 °C under 5% C0 2 for 24 h, the cells on 
upper membrane layer were gently wiped away using a 
cotton swab and the cells migrating to the lower 
membrane were fixed with 4% paraformaldéhyde. The 
number of cells in the lower side of the membranes was 
counted microscopically. 

2.8. Hypoxic-ischemic mode! and experimental groups ofrats 

As HI model was established as in the previous method 
of Rice et al. [16], The rats were anesthetized with a 10% 
chloral hydrate solution (0.3 mL/1 00 g), and the left carotid 
artery was isolated and ligated doubly. After 2-3 h, the 
waking rats were put into a transparent container which 
was ventilated with a constant flow of mixed gas 
containing 8% oxygen and 92% nitrogen for 2 h. 

The animais were randomly assigned into four experi- 
mental groups after HI treatment: sham group, Hl-vehicle 
group, Hl-GFP-MSC group, HI-HIF-la-MSC group. For the 
rats in the sham group, the carotid artery was isolated but 
not ligated, and no hypoxia exposure was performed. The 
rats in the Hl-GFP-MSC group and the HI-HIF-la-MSC 
group had received 5 x 10 5 MSCs re-suspended in 0.5 mL 
of a physiological saline solution by intravenous injection 
in the tail 24 h after HI treatment, and the other rats 
received the same volume of physiological saline solution. 

2.9. Morris water maze test 

Rats (n = 6 per group) hâve been submitted to the 
Morris water maze test [17] on day 14 after HI to evaluate 
their spatial learning and memory abilities. The maze 
consisted of a black circular tank (160 cm diame- 
ter x 50 cm height) filled with water (température: 22- 
24 °C). An automatic tracking System (San Diego Instru- 
ments, San Diego, CA, USA) was used to record the 
swimming path and time. Each rat was tested for fîve 
training days and one test day. In the task, the rat was put 
in every quadrant of the pool sequentially, and swam for 
60 s to locate the submerged platform; if it failed, the 
expérimenter would guide it to the platform. Then the rat 


was allowed to remain on the platform for 10 s. On the last 
day, the platform was removed before the test, and the 
swimming time in the former platform quadrant was 
recorded. 

2.20. Tissue préparation and histological examination 

Rats brains were obtained on days 7, 14 and 21 after HI 
(each time point n = 6 per group), and each hippocampus 
was embedded in paraffin. The paraffin sections (4 p,m 
thick) were prepared for détection of the migration of 
MSCs and histological examination in the hippocampus. To 
measure the migrated cells’ number, the sections were 
observed with a fluorescence microscope to count the CM- 
Dil positive cells. The histological examination was 
determined by hematoxylin and eosin (HE) staining. 

2.2 2. Statistical analysis 

The data obtained were presented as mean ± standard 
déviation (SD). Différences between groups were evaluate by 
analysis of variance for multiple comparisons using SPSS 
19.0. A value of P <0.05 was considered to be statistically 
signifïcant. 


3. Results 

3.2. Overexpression of HIF- 2 a in MSCs 

After transfection for two days, GFP expression 
observed under the fluorescence microscope reached 
up to 90% of the cultured MSCs (Fig. IA). In the 
meantime, western blotting was performed to investi- 
gate the expression of HIF- la in MSCs after hypoxia 
pretreatment. In the présent study, the expression of 
HIF-la in HIF-la-MSCs was approximately two times 
that in NM-MSCs and GFP-MSCs, while the expression 
level in GFP-MSCs was in accordance with NM-MSCs 
(Fig. IB and C). The results indicated that HIF-la 
overexpressed MSCs were successfully constructed by a 
recombinant lentivirus System. 


NM-MSCs GFP-MSCs HIF-la-MSCs 



B 

C 


HIF-la 

GAPDH 



93KD 

36KD 



Fig. 1. Transfection effîciency of MSCs. A. The level of GFP expressed in MSCs was observed by fluorescence microscopy (magnification x 100). B. Protein 
level of HIF-la in different groups; GAPDH was used as an inner control. C. Quantification of HIF-la protein expression showed in B. **P < 0.01 versus the 
GFP-MSCs group. 
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Fig. 2. MTT assay used to quantify the MSCs viability in different group. 
"P< 0.01 versus the GFP-MSCs group. 


3.2. HIF-la overexpression improved the viability of MSCs 

MSC viability was assessed by an MTT assay to 
détermine the effect of HIF-lot overexpression. As shown 
by our results, HIF-la overexpression immensely im- 
proved the viability of MSCs, while there was no significant 
différence between NM-MSC and GFP-MSC groups 
(Fig. 2A). The resuit implied that MSCs with overexpression 
of HIF-la had better viability. 

3.3. HIF-la overexpression enhanced migration ofMSCs both 
in vitro and in vivo 

The migration ofMSCs was analyzed both in in vitro and 
in vivo experiments. To observe the migration of MSCs 


directly, the cells were labeled with CM-Dil, a fluorescence 
marker. The results showed that more than 95% of MSCs 
were CM-Dil positive after having been labeled (Fig. 3A). 

According to the Transwell assay, the number of cells 
migrating to the lower compartment in HlF-la-MSC group 
was more than that in NM-MSC and GFP-MSC groups after 
hypoxia pretreatment, while no significant différence was 
found between NM-MSC and GFP-MSC groups after 
culturing for 24 h (Fig. 3B). 

One day after HI, rats received an intravenous 
transplantation of MSCs, and the presence of migrated 
MSCs in the hippocampus after transplantation was 
investigated by observing CM-Dil positive cells. CM-Dil 
positive cells were clearly observed in the hippocampus on 
day 7 after HI. And from day 7 to day 21 subséquent to HI, 
the number of CM-Dil positive cells in the hippocampus 
increased gradually in a time-dependent manner. A 
remarkable increase in the recruitment of HIF-la trans- 
duced MSCs was observed in the hippocampus at every 
time point compared with GFP transduced MSCs (Fig. 3C). 
Statistical analysis showed that HIF-la overexpression 
improved the migration capacity of MSCs toward the sites 
of injury. 

The quantitation of migrated cells demonstrated 
significant increases in MSCs motility in HIF-la overex- 
pressed MSCs both in in vitro and in vivo experiments in 
comparison to GFP infected cells. 

3.4. HIF-la overexpression strengthen the therapeutic 
efficiency of MSCs on HI rats 

In order to investigate the therapeutic effect of HIF-la 
transduced MSCs, the Morris water maze test and HE 



B 




Fig. 3. MSCs migration assay in vitro and in vivo experiments. A. The phase-contrast and fluorescent image of CM-Dil labeled MSCs (magnification: x 100). 
B. Transwell migration assay, performed in vitro experiments. C. Transplanted MSCs homing to the hippocampus were counted at different time points. 
**P<0.01 versus the GFP-MSCs group. 
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Fig. 4. The therapeutic efficiency of HIF-la overexpressed MSCs in HI rats. A. Percentage of swimming time in the target quadrant in the Morris water maze. 
B. HE staining in different groups of the hippocampus. Représentative images in the sham group (1), the pathological changes in the Hl-Vehicle, HI-GFP- 
MSCs and HI-HIF-la-MSCs groups at different time points (2). ## P < 0.01 versus the sham group, AA P < 0.01 versus the Hl-Vehicle group, P < 0.05 versus 
the Hl-GFP-MSCs group. 


staining in the hippocampus were performed on HI rats. On 
day 14 after HI, the spatial performance of four groups was 
evaluated by the Morris water maze test. In the spatial 
probe test, spatial learning ability was significantly 
affected by HI, as Hl-treated rats spent less time in the 
target quadrant compared to the sham group. The results 
suggested that the rats suffered from cognitive dysfunction 
due to HI treatment. There was a significant increase of 
time in the target quadrant in both MSC-transplanted 
groups compared with Hl-vehicle group, but a more 
increase in the amount of time hâve been found in HI-HIF- 
la-MSC group (Fig. 4A). These results indicated that the 
Hl-induced rats achieved partial remission on memorial 
and cognitive functions after MSCs transplantation, while 
HIF-la overexpression could promote this effect. 

Some visualized pathological changes were found in HE 
staining. The data supported that the morphological 
structures of the hippocampus in the sham group were 
normal, with uniform dyeing (Fig. 4B1). In contrast with 
the sham group, the pathological changes in the hippo- 
campus in the Hl-vehicle group revealed serious injury, the 
cells were disordered, and cavitation was visible, suggest- 
ing that the HI model had been successfully established. 
Meanwhile, the rats in HI-GFP-MSC and HI-HIF-la-MSC 
groups showed similar but slight pathological changes at 
every time point after having been subjected to HI, and the 
changes in the MSCs with HIF-la transfection group were 
less marked than those of the GFP control group (Fig. 4B2). 
The results illustrated that the transplantation of MSCs 
ameliorated pathological changes after HI, while HIF-la 
overexpression could enhance therapeutic efficiency. 

4. Discussion 

In recent years, abundant researches show that MSCs 
transplantation has been considered a promising strategy 
in Hl-induced brain-damaged treatment [4, 18,19). Howev- 
er, the insufficient viability as well as the low number of 
cells transplanted to the target tissues limits the potential 
therapeutic effects [7,8]. Recent évidences hâve indicated 
that HIF-la is important in cell viability, prolifération and 


migration following HI injury [12,13]. In the présent study, 
we therefore investigated the therapeutic effect of HIF-la 
overexpressed MSCs on HI brain damage. Firstly, overex- 
pression was achieved using a recombinant lentiviral, and 
a robust overexpression of HIF-la at the protein level has 
been seen. Additionally, the functions of HIF-la overex- 
pression in MSCs showed that HIF-la overexpression 
resulted in increased viability and migration of MSCs. 
Finally, using a HI model, we found that HIF-la 
overexpressed MSCs had a stronger therapeutic efficiency 
on HI brain-damaged treatment by mitigating the injury 
on behavioral and histological changes evoked by HI 
insults, accompanied with more MSCs migrating to the 
damaged cérébral area. 

The hypoxia-inducible factor-1 (HIF-1) is recognized as 
the master regulator of hypoxia response. It is a 
heterodimeric transcription factor composed of an induc- 
ible HIF-la subunit and a HIF-1(3 subunit, and HIF-la is a 
crucial mediator of the cellular response to hypoxia in 
contrast to HIF-1 (3, implying that HIF-la is a key mediator 
of the bénéficiai effect in a low oxygen environment 
[20]. In this study, just as expected, we hâve seen the 
expression of HIF-la after hypoxia pretreatment, while 
overexpression increased the expression. 

The présent study showed that overexpression of HIF- 
la enhanced MSCs viability in vitro, which was consistent 
with a previous study [21]. Intriguingly, the results also 
showed in our présent study that HIF-la overexpression 
could further promote the viability of MSCs in vitro in 
comparison to the effect of HIF-la induced by hypoxia 
pretreatment, indicating that the expression level of 
endogenous HIF-la, which achieved during hypoxia, does 
not reach the maximum required for HIF-la transcription 
activities, and transfection with exogenous HIF-la impro- 
ves the transcription level. 

Cérébral ischemia and hypoxia will resuit in cognitive 
dysfunction [22,23] as well as pathological changes [24], 
especially the changes in the hippocampus [25], which 
showed in the Morris water maze test and HE staining in 
our study. Transplantation of MSCs after HI contributed 
to an improved function recovery of neurological déficits 
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and reduced the degree of pathological changes 
[18,19,26]. Meanwhile, it was observed that HIF-la 
overexpression led to a better efficiency of MSC transplan- 
tation. This resuit is not in keeping with studies by Lopez 
et al. [ 27 ], which showed that overexpression of HIF-1 a did 
not hâve protective actions upon hypoxia-mediated 
neuronal death. However, enhanced therapeutic efficiency 
of cells on brain damage through upregulation of HIF-la 
also been found in other groups [28,29]. These controver- 
sial results may be explained by the different damage 
degree (HI vs. mild hypoxia). 

We also investigated the recruitment of MSCs in the 
hippocampus in time course experiments; the présent 
study results together with that from previous study [30] 
demonstrated a number of MSCs in the hippocampus 
increasing continuously three weeks after transplantation. 
Consistently with in vitro experiments, HIF-la overex- 
pression also improved MSCs migration to the injured area. 
The upregulation of HIF-la increased the migration of 
MSCs, as well as it enhanced the therapeutic efficiency of 
MSCs, which was also found by other groups [31 ,32 ]. Thus, 
enhancing the number of MSCs to the site of injury may 
improve the therapeutic efficiency of MSCs transplanta- 
tion. In other words, the improvement of the therapeutic 
efficiency through overexpression of the HIF-la of MSCs 
could be due to increased cell migration to the damaged 
area. 

The limitation of this study is that the underlying 
associated mechanisms about HIF-la promoting the 
directional migration of MSCs to the cérébral HI area need 
to be confirmed, which we will investigate in further 
studies. 

In conclusion, this study demonstrated that HIF-la 
overexpression could increase the therapeutic efficiency 
MSCs in HI, which will correlate well with the promotion 
cells’ directional migration to the cérébral HI area. We 
propose that transplantation of MSCs with HIF-la 
overexpression is an attractive therapeutic option to treat 
Hl-induced brain injury in the future. 
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A high intake of sugars has been linked to diet-induced health problems. The aim of this 
study was to assess whether the long-term consumption of a high-carbohydrate diet (HCD) 
would cause the hepatic histopathological and metabolic abnormalities that characterize 
nonalcoholic steatohepatitis (NASH) in a desert gerbil, Gerbillus gerbillus. Compared to 
natural diet, HCD leads to several metabolic disorders including adiposity, dyslipidemia, 
insulin résistance, ectopic fat déposition in the liver, which were associated with higher 
levels of transcripts of genes involved with fat synthesis, endoplasmic réticulum (ER) 
stress, and fibrosis. In the same way, the experimented animais showed enhanced 
oxidative stress. Taken together, these results demonstrate that HCD consumption in 
gerbils induces metabolic disorders and damaged liver, which are key contributors to NASH 
development. These results suggest that this rodent represents a valuable natural model 
for human diet-induced metabolic disorders and nonalcoholic fatty liver disease (NAFLD). 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 


RÉSUMÉ 


Mots clés : 

Gerbillus gerbillus 
Insulino-résistance 
Lipogenèse 
Stress R.E. 

Stress oxydatif 


Les effets d’un régime riche en hydrates de carbone (RHC) sur l’induction de 
dysfonctionnements métaboliques et d’altérations histopathologiques au niveau du foie 
caractérisant la stéatose hépatique non alcoolique (NASH) ont été étudiés chez la gerbille 
Gerbillus gerbillus. Comparativement au régime naturel, l’alimentation hyperglucidique 
induit des perturbations métaboliques marquées par une adiposité, une dyslipidémie, une 
insulinorésistance et un dépôt ectopique de graisses dans le foie, associé à l’augmentation 
des ARNm de gènes impliqués dans la synthèse des lipides, le stress du réticulum 
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endoplasmique et la fibrose. Parallèlement, les animaux expérimentaux développent un 
état de stress oxydatif. Le RHC induit chez G. gerbillus des troubles métaboliques et des 
lésions hépatiques caractéristiques de la NASH. Nos résultats suggèrent que ce rongeur 
représente un modèle naturel de choix pour l'étude des dysfonctionnements 
métaboliques et hépatiques induits par une consommation excessive de glucides chez 
l’humain. 

© 2016 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés. 


1. Introduction 

Nonalcoholic fatty liver disease (NAFLD) is the main 
cause of hepatic dysfunction in developed countries and is 
closely related to components of metabolic syndrome such 
as obesity, dyslipidemia and type-2 diabètes [1], Excessive 
accumulation of triglycéride (TG) in hépatocytes is the 
hallmark of NAFLD. The spectrum of NAFLD ranges from 
hepatic steatosis or fatty liver to nonalcoholic steatohe- 
patitis (NASH), liver fibrosis, liver cirrhosis, and eventually 
hepatocellular carcinoma (HCC) [2], The précisé mecha- 
nisms of NAFLD remain poorly understood. The “multiple- 
hit hypothesis” is currently the most recognized theory to 
explain disease development and progression, the dysré- 
gulation in lipid metabolism being involved in the first hit 
[3], It was estimated that 30 % of the TG content in NAFLD 
livers came from de novo lipogenesis, underlying the 
importance of this pathway in the etiology of NAFLD 
[4,5], De novo lipogenesis (DNL) can be triggered by 
multiple mechanisms, including increased expression of 
lipogenic enzymes by several spécifie transcription factors; 
this is particularly true for members of the SREBP family. 
One of them, SREBP-lc, Controls hepatic DNL primarily by 
régulation of expression of genes involved in DNL, lipid 
homeostasis and glucose metabolism [5,6]. Accordingly, 
hepatic expression of SREBP-lc and its target genes are 
increased in human fatty liver, compared to healthy 
individuals [7,8], although DNL is an important détermi- 
nant of steatosis [9], The second of the two hits could be 
due to (1 ) oxidative stress, (2) proinflammatory cytokines, 
(3) mitochondrial dysfunction, or/and (4) endoplasmic 
réticulum stress. Recently, accumulating data hâve impli- 
cated the disruption of endoplasmic réticulum (ER) 
homeostasis, or ER stress, in both the development of 
steatosis and progression to NASH [10,11]. ER stress may 
lead to the activation of various intracellular stress 
pathways that can initiate or exacerbate insulin résistance 
(IR) and inflammation and, in some cases, culminate in 
hépatocyte cell death and liver damage, ail of which are 
important in the pathogenesis of NAFLD. In spite of 
growing knowledge, several aspects of NAFLD pathogene- 
sis are still unknown. 

Considering the difficulty in developing human studies 
to evaluate the influence of nutrition in the development of 
NAFLD and associated metabolic abnormalities, animal 
models constitute a reliable alternative way. Different 
animal models of NAFLD/NASH hâve been developed, but 
few of them replicate the entire human phenotype 
[12,13], These models may be classified into three basic 
categories: those caused by either spontaneous or induced 
genetic mutation; those produced by either dietary or 
pharmacological manipulation; and those involving genetic 


mutation and dietary or Chemical challenges. The dietary 
manipulations used in these last two types of models 
usually do not resemble the human dietary pattern. 

The aim of this study was to détermine the long-term 
impact of high-carbohydrate diet on liver morphology and 
function in a desert rodent ( Gerbillus gerbillus ). We 
therefore evaluated the markers involved in metabolic 
functions, i.e. lipogenesis, fibrosis, ER stress, histopatho- 
logical changes, and oxidative stress. In the présent study, 
we developed a model of obesity and obesity-related 
NAFLD in a desert rodent (G. gerbillus') using a simple 
carbohydrate-rich diet. 

2. Materials and methods 

2.3. Animais and diets 

G. gerbillus individuals were collected from the semi- 
desert Algerian région of Beni-Abbes (30°7’ latitude north 
and 2 = 10’ longitude west). The authorization to capture the 
animais in desert région was given by the Ministry of 
Higher Education (Algeria). G. gerbillus specimens were 
maintained under controlled température (22 ± 1 °C), 
humidity (50%) with a fixed 12-h light/dark cycle. After 
two weeks of acclimatation, the adult gerbils of both sexes 
were randomly divided into two groups. The control group 
(n = 6) received a natural diet composed of halophile fresh 
plants, seeds, dry plants. . . whereas the experimental group 
(n = 6) received an HC diet (25 % of barley and 75 % of dried 
dates) corresponding to a daily calorie intake of 22.5 calories/ 
animal. Details of the composition of the high-carbohydrate 
diet (HCD) are presented in Table 1. After six months of diet 

Table 1 

Diet compositions according to Nicole Tonelli and François Gallouin, Fruit 
and seeds comestibles worldwide, 7th édition. MedPharm Scientific 
Publishers/Taylor & Francis, 2008. 


Nutritional composition 

Barley 

(g/100 g) 

Dried dates 

(g/100 g) 

Carbohydrates 

63.30 

65 

Starch 

61.59 

- 

Simple sugars 

1.71 


Glucose 

- 

25 

Fructose 

- 

25 

Saccharose 

- 

14 

Sorbitol 

- 

1 

Protein 

11.2 

2 

Lipids 

2.1 

0.5 

Water 

12.2 

20 

Fibers 

9.8 

9 

Vitamins 

0.015 

0.012 

Minerais & Trace Eléments 

0.956 

0.862 

Calories 

314 kcal 

275 kcal 
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administration, the animais were killed by décapitation. The 
livers and fat pads (mesenteric. génital and perirenal) were 
totally removed, rinsed with cold 0.9 % NaCl, weighed, and 
eut. One part of the liver was fixed in 10% neutral buffered 
formalin for histopathological investigations. The other part 
was first frozen in liquid nitrogen and then stored at-80 C for 
later use in biochemical studies. 

2.2. Biochemical analyses 

Blood samples were collected after 24 weeks of feeding 
by puncture from the retro-orbital venous plexus into 
heparinized or EDTA-containing tubes. The plasma was 
immediately separated by centrifugation at 3000 g for 
10 min and rapidly stored at-20°C until analysis. Plasma 
glucose, total cholestérol, triacylglycerol (TG), HDL choles- 
térol and reliable markers of hepatic function (aspartate 
aminotransferase, AST ; alanine aminotransferase, ALT) were 
measured by using the corresponding kits for analysis 
(Spinreact, Spain). Plasma LDL cholestérol concentration 
was calculated by using the équation of Friedewald. Plasma 
insulin was measured by radioimmunoassay (1NSUL1N-CT 
kit, Cisbio assay, France). Homeostasis model assessment of 
insulin résistance (HOMA-1R) was calculated according to 
the following formula: HOMA-IR= fasting insulin (p,U/ 
mL) x fasting glucose (mmol/L)/22.5. 

2.3. Histopathology 

A portion of liver (100 mg) was fixed with 10% formalin 
for 24 h, dehydrated in graded éthanol, and embedded in 
paraffin. The samples were eut into 5-|xm-thick sections 
and stained with Masson’s trichrome for évaluation by 
light microscopy for fibrotic déposition. 

2.4. Hepatic lipid analyses 

Liver total lipids were extracted by the method of Folch 
et al. [14]. The concentrations of triglycérides and choles- 
térol in the hepatic lipid extracts were measured using the 
same enzymatic kits as those used for plasma analysis. 

3. Estimation of hepatic oxidative stress 

3.1. Liver homogenate préparation 

Liver samples from each gerbil were homogenized in 
Tris/EDTA buffer (50 mM Tris. 0.1 mM EDTA. pH 7.6) to 
obtain a 10% (w/v) liver homogenate. An aliquot of tissue 
homogenate was centrifuged at 1 ,000 rpm for 5 min at 4 °C 
for the estimation of MDA (malondialdehyde). The remain- 
ing volume of homogenate was centrifuged at 9,600 rpm for 
20 min at 4 °C. The supernatant was collected and used for 
the estimation of GSH (reduced glutathione), and for 
measuring the activity of SOD (superoxide dismutase), 
CAT (catalase), and GST (glutathione-S-transferase). 

3.2. Markers of lipid peroxidation 

Lipid peroxidation was assessed indirectly by the 
measurement of secondary products malondialdehyde 


(MDA) using a thiobarbituric acid (TBA) method according 
to the protocol. MDA formed a colored complex in the 
presence of thiobarbituric acid, which was detected by 
spectrophotometrical measurement of the absorbance at 
535 nm [15j. 

3.3. Evaluation of antioxidant enzyme activity 

3.3.1. Activity of SOD 

Superoxide dismutase was estimated using the method 
described in [16], Superoxide anions are generated during 
the oxidation of pyrogallol under test conditions and SOD 
enzyme inhibits the oxidation of pyrogallol. The enzyme 
activity is measured by monitoring the rate of decrease in 
optical density at 420 nm. One unit of enzyme was defined 
as the amount of enzyme required to produce 50% 
inhibition of pyrogallol autoxidation under the assay 
conditions and expressed as U/mg protein. 

3.3.2. Activity of catalase 

Catalase activity was assayed by the method of Aebi 
[17]. The rate of H 2 0 2 décomposition was followed by 
monitoring absorption at 240 nm. The enzyme activity was 
expressed as millimoles of H 2 0 2 consumed per minute per 
milligram of protein. 

3.3.3. Activity of GST 

GST activity was measured using the method of Habig 
et al. [ 1 8 ]. Briefly, 1 mM CDNB was added to a buffer solution 
containing 1 mM GSH and an aliquot of the sample to be 
tested. Upon addition of CDNB, the change in absorbance at 
340 nm was measured as a function of time. The extinction 
coefficient for this reaction is 9.6 mM^’ cra -1 . 

3.3.4. Concentration of GSH 

GSH level was estimated using a colorimétrie technique 
as mentioned by Ellman [19], based on the development of 
a yellow color when 5,5-dithio-bis-(2-nitrobenzoic acid) 
(DTNB) is added to compounds containing sulfhydryl 
groups. Absorbance was recorded at 41 2 nm. The total GSH 
content was expressed as nanomoles of GSH per milligram 
of protein. 

3.4. Protein estimation 

The protein content was determined using Bradford’s 
method [20], bovine sérum albumin being taken as the 
standard. 

3.5. RNA extraction and RT-qPCR 

Total RNA was isolated according to [21]. Real-time 
quantitative RT-PCR analyses were performed with 50 ng 
of cDNA/3 mM MgC12/250 nM sense and antisense primers 
(Proligo, Boulder, CO, USA) in a final reaction volume of 
25 |xl by using the qPCRTM Core Kit (Eurogentec, Brussels) 
and the MyiQ real-time PCR détection System (Bio-Rad). 
Spécifie primers were designed with Primer Express 
software (Table 2). The relative quantitation of each gene 
was calculated after normalization to 18S ribosomal RNA 
by using the comparative CT method. 
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Table 2 

Sequences of the primers used for quantitative real-time PCR. 


Genes 

primer sequenses 

FAS 

Sense 

5'-tgc tcc cag ctg cag gc -3' 

AntiSense 

5'-gcc cgg tag etc tgg gtg ta-3' 

SREBPlc 

Sense 

5'-gga gcc atg gat tgc aca tt-3' 

AntiSense 

5'-ggc ccg gga agt cac tgt-3' 

ACC 

Sense 

5'-tgg gea cag acc gtg gta g-3' 

AntiSense 

5'-ggt ctt aaa tgc aga gtc tgg gaa-3' 

SREBP2 

Sense 

5'-ccc ttg act tcc ttg ctg ca-3' 

AntiSense 

5'-gcg tga gtg tgg gcg aat c-3' 

HMG CoA synthase 

Sense 

5'-gcc gtg aac tgg gtc gaa-3' 

AntiSense 

5'-gca tat ata gea atg tet cct gea a-3' 

HMG CoA reductase 

Sense 

5'-gat tet ggc agt cag tgg gaa-3' 

AntiSense 

5'-gtt gta gcc gcc tat gct cc-3' 

aSma 

Sense 

5'-ttg gaa aag atc tgg cac cac-3' 

AntiSense 

5'-gca gta gtc aeg aag gaa tag-3' 

CollAl 

Sense 

5'-ccc cgg gac tcc tgg act t-3' 

AntiSense 

5 / -gctccgacacgccctctctc-3' 

TGF- fl 

Sense 

5'-gcctgagtggctgtcttttgac-3' 

AntiSense 

5'-catggatggtgcccaggt-3' 

GRP78 

Sense 

5'-gaa agg atg gtt aat gat gct gag-3' 

AntiSense 

5'-gtc ttc aat gtc ege atc ctg-3' 

ORP150 

Sense 

5'-tgt cct ctt ggc aga cct gtt g-3' 

AntiSense 

5'-ttt tcc tcc gag att cct tgt tc-3' 

18S 

Sense 

5'-ggg age ctg aga aac ggc-3' 

AntiSense 

5'-ggg teg gga gtg ggt aat tt-3' 


3.6. Préparation of nuclear extracts 

Nuclear extracts were prepared as described in 
[22]. Livers were removed and quickly rinsed in an ice- 
cold phosphate-buffered saline solution with pepstatin 
(5 pig/ml), leupeptin (5 p,g/ml), and aprotinin (2 ptg/ml) as 
protease inhibitors. The livers were then transferred into a 
beaker, on ice, containing the homogenization buffer (2 M 
sucrose, lOmM Hepes, pH 7.6, 25 mM KC1, 0.15 mM 
spermine, 0.5 mM spermidine, 1 mM EDTA, 10% glycerol, 
0.5 mM dithiothreitol) and minced with scissors. After 
Processing in a Dounce apparatus in fresh buffer on ice, the 
homogenates were layered on 2 M sucrose cushions and 
centrifuged at 80,000 g for 35 min at 0 °C. The nuclear 
pellets were then rinsed in a buffer containing lOmM of 
Tris, pH 7.4, 10 mM ofNaCl,3 mM ofMgCl 2 , resuspended in 
a 20 mM Hepes, pH 7.9, 0.4 M NaCl, 1 mM EDTA, 1 mM 
ECTA, 1 mM dithiothreitol buffer and incubated 15 min on 
ice. After spinning 5 min at 1 0,000 rpm, aliquots of the 
supernatants were stored at -80 °C. The protein content 
was determined as described by Bradford, using bovine 
sérum albumin as a standard. 

3.7. Statistical analysis 

Data were statistically analyzed by performing a 
nonparametric Mann-Whitney test using GraphPad Prism 


to compare experimental groups (v. 6.00 for Windows, San 
Diego, CA, USA). Data were expressed as mean < remove 
picture pageno 5> standard error of the mean (SEM) and 
considered significant at P< 0.05. 

4. Results 

4.1. Body and liver weights, bodyfat accumulation and 
biochemical parameters 

The administration of a high-carbohydrate diet for 
6 months resulted in significant increases in the final body, 
liver and abdominal fat weights compared with those in the 
gerbils from the control group (Table 3). The liver and fat 
indices increased 1.08 and 2.5 times, respectively, com- 
pared to the control diet, after a 24-week HC diet. The HCD 
group showed higher plasma concentrations of triglycér- 
ides (P < 0.001) and total cholestérol (P < 0.05). However, 
the mean plasma HDL cholestérol and LDL cholestérol levels 
were not statistically different between the two groups. 
Plasma glucose was not modified in experimental gerbils 
(P = 0.123), whereas plasma insulin level and HOMA-1R 
index were significantly increased (P < 0.05) with the HCD 
compared to the control diet (Table 3). The 24-week feeding 
of HCD to gerbils did not significantly affect ALT and AST 
activities (markers of liver function) (Table 3). 

4.2. Liver TG and cholestérol levels 

The liver TG and TC levels are shown in Table 3. HCD 
feeding significantly increased the hepatic TG levels 
(188.4 ± 26.20 mg/g) compared with the control group 
(82.97 ± 12.07 mg/g; P < 0.01). The hepatic TC content in 
the HCD group was increased two times compared to control 
animais. 


Table 3 

Physiological and plasma parameters in control and HCD gerbils after 
24 weeks of treatment. 



Control 
(n = 6) 

Experimental 
(n = 6) 

Body weight (g) 

27.68 ±0.88 

31.43 ± 0.70** 

Fat pads (mesenteric, génital 

0.49 ± 0.03 

1.35 ±0.15 ‘ 

and perirenal) (g) 

Fat to body weight ratio (%) 

1.76 ±0.10 

4.35 ±0.57 

Liver weight (g) 

0.98 ± 0.005 

1.19 ±0.005 

Liver to body weight ratio (%) 

3.51 ±0.21 

3.79 ±0.11 

Hepatic triglycérides (mg/g liver) 

82.97 ±12.07 

188.4 ±26.20 

Hepatic total cholestérol 

41.45 ±8.20 

86.86 ±3.85'' 

(mg/g liver) 

Plasma triglycérides (mg/dl) 

93.91 ±10.04 

208 .7 ±18.48 

Plasma total cholestérol (mg/dl) 

103.1 ±7.05 

136.8 ±9.03 

Plasma HDL-c (mg/dl) 

57.07 ±3.55 

51.52 ±5.20 

Plasma LDL-c (mg/dl) 

28.81 ±8.95 

42.33 ±8.14 

Plasma glucose (mg/dl) 

94.74 ±3.73 

107.8 ±6.78 

Plasma insulin (mU/mL) 

48.38 ±2.21 

183.9 ±47.26 

HOMA-1R index 

11.34 ±0.79 

48.91 ±12.78 

ASAT (U/l) 

32.08 ± 7.02 

56.84 ±25.71 

ALAT (U/l) 

18.95 ±2.68 

24.78 ±4.17 


Data are expressed as means ± SEM. 

* P <0.05. 

" P <0.01. 

P <0.001 versus the control group. 
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Fig. 1. Représentative Masson Trichrome stained liver sections from gerbils fed with normal diet (A) or high carbohydrate diet (B-D) during 24 weeks. 
Trichrome Masson staining was performed to detect fibrotic areas (in green, indicated by arrows). Sections were observed at 400-fold magnification (scale 
bars: 50 p,m). 


4.3. Histological analysis of the liver tissues 

Représentative photomicrographs of liver sections are 
described in Fig. 1. Liver sections of control animais 
showed a normal liver architecture formed of hepatic 
lobules. Briefly, each lobule is made up of radialing plates. 
Strands of cells form a network around a central vein 
(Fig. IA). The histological examination of liver tissue of 
gerbils submitted to the HC diet revealed variable changes. 
These changes were evidenced by a marked micro and 
macro steatosis characterizing the NAFLD status (Fig. IB). 
Trichrome staining showed a marked collagen déposition 
within the pericellular space at sites where fatty degener- 
ation of hépatocytes occurred (Fig. IC). Also, marked 
mononuclear cell infiltrâtes in the portai System space (red 
arrows) were observed (Fig. 1D). 

4.4. Lipid peroxidation of the liver 

As shown in Fig. 2A, the level of malondialdehyde 
(MDA), an index of lipid peroxidation in liver homogenates 
was significantly increased by 67% in the HCD group (P < 
0.05) compared to the control group. 


4.5. Activities of antioxidant enzymes in liver tissue 

The effect of a long-term administration of HCD on 
activities of liver antioxidant enzymes is depicted in 
Fig. 2(B-E). GSH (reduced glutathione) levels were 
statistically significantly lower in the HCD group than in 
the control group (Fig. 2B). SOD (superoxide dismutase) 
activities did not yield statistically significant différences 
due to the large variability observed (HCD: 
(0.600 ±0.129 U/mg protein) vs. (0.914 ±0.331 U/mg 
protein) in control animais). As illustrated in Fig. 2A and B, 
the GST (glutathione-S-transferase) and CAT (catalase) 
activities were unchanged between these two groups. 

4.6. Hepatic gene expression in gerbils 

4.6. 1. mRNA expression of SREBPic. ACC and FAS 

Sterol regulatory element-binding protein-1 (SREBP-1) 
is a nuclear transcription factor that transcriptionally 
activâtes nearly ail the genes involved in de novo 
lipogenesis [23]. Acetyl-CoA carboxylase (ACC) and fatty 
acid synthase (FAS) are important downstream target 
genes of SREBP-1 and are the rate-limiting enzymes of 
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Fig. 2. Oxidative stress in liver tissue. The intrahepatic contents in MDA(A), GSH (B), SOD(C), CST (D)and CAT(E) were respectively measured in frozen liver 
homogenate of control or treated gerbils. * P < 0.05. """ P < 0.0001 versus control group. 


palmitate formation. Compared with Controls, HCD feeding 
alone significantly up-regulated the hepatic mRNA expres- 
sion of SREBPl-c (P < 0.01), FAS (P < 0.001) and ACC 
(P < 0.01) (Fig. 3A-C). 


4.6.2. mRNA expression of SREBP2, HMG CoA synthase and 
HMG CoA reductase 

SREBP-2 preferentially régulâtes target genes involved 
in cholestérol biosynthesis and uptake in the liver, such as 
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Fig. 3. mRNA expression for genes involved in fatty acid synthesis in the liver of gerbils fed normal diet or high carbohydrate diet for 24 weeks. ’* P < 0.01 . 
"" P < 0.001 vs. the corresponding values in control gerbils. 


3-hydroxy-3-methylglutaryl-CoA synthase (HMGCS) and 
3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) 
[24], The mRNA expressions of SREBP-2 and HMG CoA 
synthase were signifïcantly increased in HCD gerbils as 
compared to Controls (Fig. 4A-B). However, HMGCoA 
reductase mRNA expression was not signifïcantly different 
between the groups, due to the large variability observed, 
but tended to increase as compared to Controls (Fig. 4C). 

4.6.3. mRNA expression of a-SMA, CollAI and TGF-/3 

Our results show that mRNA gene expression for a-SMA 
(a-smooth muscle actin, a gene related to stellate cell 
activation) was approximatively two-fold higher in HCD 
group compared to control group (P < 0.05). Similarly, the 
hepatic mRNA expression of CollAI (collagen 1A1) was 
up-regulated (1.5-fold) in HCD group compared to control 
group (Fig. 5B). However, liver expression of TGF-p 
(transforming growth factor) remained unchanged be- 
tween the two groups (Fig. 5C). 

4.6.4. mRNA expression of GRP7 8 and ORP150 

Abnormal endoplasmic réticulum (ER) stress contrib- 

utes substantially to hépatocyte cell death during altered 
lipid metabolism [25], To investigate whether ER stress 


was associated with the pathogenesis of hepatic steatosis, 
we studied GRP78 (glucose-regulated protein 78) by RT- 
PCR. GRP78 acts as an apoptotic regulator by protecting or 
promoting cell death during ER stress [26]. The RT-PCR 
results indicated that the GRP78 mRNA expression level 
was elevated in the experiment group, compared with the 
corresponding control groups (+ 44 %) (Fig. 6A). The 150- 
kDa oxygen-regulated protein (ORP150) is an anti- 
apoptotic ER résident chaperone that exerts a protective 
effect against ER stress-dependent apoptosis [27], As 
shown in Fig. 6B, ORP150 mRNA expression was signifi- 
cantly increased in the HCD group compared with the 
control group (P < 0.05). 

5. Discussion 

In this study, after 24 weeks of feeding with HCD (dates 
and barley), gerbils exhibited increased body weight, 
viscéral fat mass, viscéral fat/body weight ratio, and liver/ 
body weight ratio when compared with control animais. 
HCD treatment also caused significant incréments in the 
plasma TC and TG levels, elevated hepatic TG and TC 
contents, and marked lipid droplet déposition in the liver 
sections, which are consistent with signs of fatty liver. At 


Relative HMGCoA reductase gene expression Relative HMGCoA synthase gene expression Relative SREBP2 gene expression 
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this time, the liver enzymes (ALT and AST) levels were 
normal, suggesting that the presence of NASH does not 
necessarily correlate with higher levels of these liver 
transaminases. Indeed, ALT levels are persistently normal 
in more than one half of the patients with NAFLD and 
biopsy proven NASH [28], It is interesting to note that after 
24 weeks of dietary treatment, the HCD group exhibited no 
significant changes (but a 14% increased trend) in the 
plasma glucose level. In addition, HCD induced hyper- 
insulinemia and insulin résistance (confïrmed by HOMA-IR 
calculation). 

Body fat distribution appears to be even more impor- 
tant than the total amount of adipose tissue, and the 
viscéral fat mass is strongly linked to insulin résistance and 
NAFLD [29]. Viscéral fat-released FFAs are transported to 
the liver by the portai vein and may contribute to hepatic 
steatosis, production of triglycéride rich VLDL, and 
elevated (3-oxidation [30], 

A large body of evidence supports a complex interaction 
between NAFLD and insulin résistance [31], Some studies 
hâve suggested that the abdominal adipose tissue plays an 
important rôle in the development of insulin résistance 

[32] , Furthermore, viscéral adipose tissue (VAT), a harmful 
fat deposit, has been considered to induce liver insulin 
résistance and further induce systemic insulin résistance 

[33] , 

This study demonstrates that a high-carbohydrate diet 
was able to induce obesity-related NAFLD in our experi- 
mental model. Animal models in which NAFLD was 
induced by simple carbohydrate-rich diets (usually fruc- 
tose) are less numerous, and in most of them only hepatic 
steatosis was observed [34,35]. Hyperinsulinemia led to 
increased hepatic synthesis of fatty acids, triglycéride 
accumulation in the hépatocytes, with subséquent stea- 
tosis. Several in vivo and in vitro studies hâve indicated that 
hepatic TG déposition may not be harmful; rather, it may 
represent a protective mechanism against FFA-induced 
lipotoxic liver injury by storing FFAs in the form of TGs 
[36], 

Increased calorie intake, especially refined sugar and 
fructose, correlates with increases in dyslipidemia, IR, 
and nonalcoholic fatty liver disease (NAFLD) [37]. NAFLD 
is the first step in the hepatic diseases that can evolve 
into steatohepatitis (NASH), with inflammatory infiltra- 
tion resulting in cirrhosis, and even hepatocarcinomas 
[38], 

Epidemiologie studies suggested that a diet high in 
refined carbohydrates (HCD) (50-65% of energy from 
carbohydrates, high glycémie index, and low in fiber) could 
promote NAFLD progression [39]. Our work demonstrates 
that gerbils fed a high glycémie index diet (GI for dates is 
70) for 24 weeks had greater hepatic TGs and adiposity 
than gerbils fed a low glycémie index diet. Our results 
fit with those of Scribner et al. [40], In addition, 
several experiments showed that the high-glycemic index 


Fig. 4. mRNA expression of SREBP2 and SREBP2 (HMGCoA synthase, 
HMGCoA reductase) target gene involved in cholestérol metabolism in 
the liver of gerbils fed with normal diet or high-carbohydrate diet for 
24 weeks. * P < 0.05. *’ P < 0.01 vs. the corresponding values in control 
gerbils. 
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carbohydrate has a close relationship with obesity, IR, and 
increased plasma and hepatic TGs [40], The grade of 
hepatic steatosis seems to be associated with the glycémie 
index of the diet, regardless of total energy or carbohydrate 
intake [39]. The underlying mechanism is not fully 
elucidated, but some hypothèses hâve been proposed. 
Foods with high glycémie indices enhance the hepatic 
influx of glucose. The excess of hepatic glucose exceeds the 
ability of glycogen production; therefore, this carbohy- 
drate will be used for the synthesis of new TG through DNL 
within the hépatocytes. An elevated-glycemic index food 
might augment oxidative stress, which can contribute to 
NASH development [41]. 

Similarly, Maersk et al. [42] found that a daily intake of 
1 L of a sucrose-sweetened soft drink for six months 
increased liver lipids in overweight participants compared 
with the intake of the same amount of milk, diet cola, or 


Fig. 5. mRNA expression for genes related to hepatic fibrosis in the liver of 
gerbils fed with normal diet or high-carbohydrate diet for 24 weeks. 
* P < 0.05 vs. the corresponding value in control gerbils. 
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water. Carbohydrate-induced accumulation of fat in the 
liver can lead to nonalcoholic fatty liver disease (NAFLD) 
[43], 

The composition of diet and especially the amount of 
glucose/starch may influence the health effects of fructose. 
As fructose is présent with glucose in most food products, 
it is more practical and relevant to look at the effects of 
fructose and glucose together than at the effects of fructose 
alone. A larger increase in DNL after eating fructose and 
glucose together (50:50 glucose:fructose) rather than the 
same amount of pure glucose has been shown [44]. Eating 
glucose with fructose is likely to affect fructose’s health 
effects by stimulating the flow of fructose to DNL [45], This 
effect could be due to both increased absorption capacity 
for fructose when co-ingested with glucose and therefore 
greater availability of fructose carbon atoms going towards 
DNL and increased blood insulin levels when glucose is 
présent in the diet. Insulin stimulâtes DNL directly and 
indirectly by inhibiting other important metabolic path- 
ways for fructose, such as gluconeogenesis. 

In the présent study, HCD-fed gerbil had severe hepatic 
steatosis. This HCD-induced elevated steatosis was associ- 
ated with the significant upregulation of DNL genes, 
including FAS (an enzyme that catalyzes long-chain fatty 
acid biosynthesis through the condensation of acetyl-CoA 
and malonyl-CoA and ACC, as well as SREBP (a key 
regulator of fatty acid synthesis in the liver gene 
expression). This evidence is consistent with previous in 
vivo studies and clinical observations that HCD-induced 
NAFLD was associated with elevated hepatic DNL 
[39]. Donnelly et al. [4] observed that patients with NAFLD 
hâve elevated fasting DNL compared with healthy 
individuals. Altering dietary macronutrient contents can 
induce differential metabolic conséquences in the liver 
[46]. Excessive dietary carbohydrates induced glucose 
conversion to FA by DNL [46]. It is important to note that 
diets rich in complex carbohydrates and fibers did not 
induce clinical DNL in humans, suggesting a potential 
interaction between digestibility/structure of carbohydra- 
tes and DNL [47]. The deleterious effect of our diet is 
attributed to dry dates that contain essentially simple 
sugars. 

The liver is the major organ involved in cholestérol 
metabolism. Sterol regulatory element-binding protein-2 
(SREBP-2) promûtes the expression of target genes 
involved in cholestérol biosynthesis and uptake in the 
liver, such as 3-hydroxy-3-methylglutaryl-CoA synthase 
(HMGCS) and 3-hydroxy-3-methylglutaryl-CoA reduc- 
tase (HMGCR) [24], The data from the présent investiga- 
tion (Fig. 4) showed that the SREBP-2 precursor and 
nuclear active forms as well as the HMGCR and HMGCS 
target genes were significantly up-regulated by HCD in 
gerbil. 

Previous research has shown that a fatty liver induced 
by a high saturated fat diet shows signs of endoplasmic 
réticulum (ER) stress in a rodent model [48j. The 
endoplasmic réticulum is an important organelle respon- 
sible for protein synthesis, folding, assembly and transport 
in eukaryotic cells. GRP78 is a member of the HSP70 family 
abundant in the lumen of ER, and plays important rôles as a 
molecular chaperone that removes misfolded proteins in 


ER lumen. Overexpression of GRP78 has been shown to 
increase ER stress résistance and to hâve bénéficiai effects 
in several cell types [49]. In addition, the survival response 
activâtes genes that encode ER-residing chaperones such 
as GRP78/Bip, which uses energy derived from ATP 
hydrolysis to prevent the aggregation of ER proteins, and 
is considered the classical marker of UPR activation. 

To confirm the involvement of ER in fatty liver 
developed in our experimented gerbils, we studied 
GPR78 and ORP150 by quantitative RT-PCR, and the data 
in Fig. 5 show that GRP78 moderately increased in 
G. gerbillus liver. It was demonstrated that obesity is 
associated with the induction of ER stress predominantly 
in the liver and adipose tissues [50], Previous research has 
shown that fatty liver induced by a high sucrose diet and 
high saturated fat diet shows clear ER stress events in a 
rodent model [48], 

In the current study, we observed a coordinated 
transcriptional upregulation of lipogenesis genes together 
with ER stress. Lipid accumulation and lipotoxicity are 
known to positively regulate the ER stress response 
through a cycle that leads to further toxic insults upon 
the liver [51], 

Abnormal ER stress contributes substantially to hépa- 
tocyte cell death during altered lipid metabolism 
[24], Thus, modulation of hepatic ER stress becomes 
crucial in the understanding of the extent of its contribu- 
tion to the pathogenesis of various hepatic disorders. 
Induction of ORP150 exhibits protective rôles in the 
prévention of ER stress [27] and dysrégulation of calcium 
homeostasis and apoptosis. Our work showed that the 
levels of ORP150 in the liver of control gerbils were 
significantly greater compared with that of experimented 
gerbils (Fig. 5B), which supports the spéculation that 
ORP150 expression could be bénéficiai under the local 
challenging conditions of ER stress. 

The induction of ORP150 médiates the amelioration of 
hepatic ER stress and apoptosis by AMPK in vitro and in 
vivo. [52] Collectively, these results indicate that chronic 
HC diet altered, at least partially, endoplasmic réticulum 
functions and induced stress conditions. 

Collagen 1A1 (CollAl) is the major hepatic collagen 
subtype associated with NASH in rodents and humans 
[53]. HC-feeding for 24 weeks greatly elevated biochem- 
ical and histological markers of fibrosis. The fibrotic 
status of the liver is governed by a balance between 
extracellular matrix (ECM) déposition and turnover. The 
biochemical (expression of aSMA) and histological 
evidence of fibrosis (Fig. 1) in HC fed gerbils provide 
clear evidence that the rate of collagen production 
exceeds that of its dégradation. 

Although it has been widely proposed that stellate cells 
are the main contributor to ECM deposits in the liver, it is 
possible that other cell types are involved [54], MMPs are a 
family of zinc-dependent protéinases that are responsible 
for the dégradation of ECM proteins. MMP dysrégulation 
has been implicated in the pathophysiology of obesity and 
diabètes. 

Overwhelming evidence indicates that the perisinusoi- 
dal hepatic stellate cells (HSCs) are the major source of 
liver fibrosis of any etiology [55]. During liver injury, HSCs 
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transdifferentiate from their quiescent physiologie to the 
fibrogenic phenotype. Such a transformation is induced by 
inflammatory mediators (e.g., TNFa) and reactive oxygen 
species and apoptotic bodies arising from dying hépato- 
cytes, and activated HSCs are characterized by expression 
of aSMA, a marker not found in quiescent HSCs [56,57]. 

Activation of HSCs is further characterized by the 
overexpression of a-SMA, MMP-2, and MMP-9 [58], Al- 
though the rôle of HSC activation in nonalcoholic fatty liver 
disease (NAFLD) has not been completely clarified, several 
studies hâve reported increased HSC activation in nonal- 
coholic steatohepatitis (NASH) [59], 

Several studies hâve demonstrated that insulin résis- 
tance is associated with advanced stages of fibrosis in 
NAFLD [60]. Because insulin promûtes HSC activation and 
insulin sensitizers can attenuate hepatic fibrosis in NASH, 
it has been suggested that insulin résistance plays an 
important rôle in NASH-related fibrogenesis [61], It is 
understood that oxidative stress induces the activation of 
HSCs in NASH [62], 

In addition to liver fibrosis, it is well-established that 
oxidative stress also is very important in the development 
and progression of NASH. It has been shown that chronic 
oxidative stress, generated through the oxidation of 
cytotoxic free fatty acids, may lead to cytokine upregulation 
and déplétion of hepatic antioxidant levels [63]. In addition, 
enhanced lipid peroxidation leads to the génération of 
reactive by-products, such as MDA, which hâve been shown 
to further stimulate cytokine production [64]. In the current 
study, there was a significant decrease in GSH content and 
SOD activity in HC-fed gerbils, and a significant increase in 
MDA levels. These findings are in agreement with those of 
Leclercq [65], who reported that diets with high amounts of 
simple carbohydrates induce hypertriglyceridemia, result- 
ing in a réduction of the antioxidant reserves [65], 

Most previous studies reported alterations in the 
expression and activity of antioxidant enzymes in the 
liver of fructose-fed rodents, although the course of these 
alterations varied largely. Some authors observed that a 
fructose-rich diet reduced the antioxidant capacity and 
caused oxidative damages in the liver [66], whereas others 
reported the absence of oxidative stress [67] or adaptation 
to modest oxidative stress [68] in fructose-fed rats. 
Alternatively, Girard et al. [69] noticed that a fructose- 
rich diet enhanced the total antioxidant capacity of the 
liver. 

6. Conclusion 

In summary, the présent study demonstrates that a 
long-term (24-week) consumption of a high-carbohydrate 
diet induces adiposity, insulin résistance, up-regulation of 
de novo lipogenesis, ER stress, oxidative stress, and finally 
fatty liver. Altogether, these results suggest that G. gerbillus 
represents a valuable natural model for human diet- 
induced metabolic disorders and NAFLD. 
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Géographie populations of the two main sandflies généra présent in Thailand were studied 
for species and population identification. Size and shape of Phlebotomus stantoni and 
Sergentomyia hodgsoni from different island and mainland locations were examined by 
landmark-based géométrie morphometrics. Intraspecific and interspecific wing compari- 
son was carried out based on 1 2 anatomical landmarks. The wing centroid size of P. stantoni 
was generally larger than that of S. hodgsoni. Within both species, wings from the continent 
were significantly larger than those from island populations. Size variation could be 
significant between géographie locations, but could also overlap between généra. The wing 
venation geometry showed non-overlapping différences between two species. The within- 
species variation of géométrie shape between different geographical locations was highly 
significant, but it could not interfère with the interspecies différence. The lack of species 
overlapping in shape, and the high discrimination between géographie populations, make 
géométrie shape a promising character for future taxonomie and epidemiological studies. 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 


RÉSUMÉ 


Mots clés : 
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Quelques populations géographiques des deux genres principaux de phlébotomes présents 
en Thaïlande, Phlebotomus et Sergentomyia, ont été examinées sur la base de la géométrie 
alaire. En utilisant 12 repères anatomiques de l’aile, la taille et la forme des phlébotomes 
ont été comparées entre des populations insulaires et continentales. Deux espèces ont été 
examinées : Phlebotomus stantoni et Sergentomyia hodgsoni. La taille des ailes de P. stantoni 
était généralement plus importante que celle de S. hodgsoni, mais dans les deux espèces, les 
ailes des populations du continent étaient significativement plus grandes que celles des 
populations insulaires. La variation de la taille, même si elle pouvait être significative entre 
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certaines localités, se chevauchait entre elles, et même quelque peu entre les deux genres. 
La géométrie de la nervation alaire n'a montré aucun chevauchement entre ces derniers. Au 
sein de chaque espèce, la géométrie de l'aile a montré des divergences très significatives 
entre zones géographiques. Cette distinction nette entre espèces et l’excellente 
discrimination géographique dans chacune d’elles font de la forme géométrique de l'aile 
un complément utile au diagnostic morphologique des phlébotomes, et un caractère 
prometteur pour les futures études épidémiologiques et biologiques de ces insectes. 

© 2016 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés. 


1. Introduction 

Phlebotomus stantoni (Newstead 1914) and 
Sergentomyia hodgsoni (Sinton 1933) are sandflies that 
belong to subfamily Phlebotominae (Diptera: Psychodi- 
dae). Sandflies are vectors of various pathogens, among 
which Leishmania sp.\ therefore, their correct taxonomie 
identification is important for epidemiological investiga- 
tions. Phlebotomus Rondani & Berte and Sergentomyia 
Franca & Parrot are Leishmania protozoa vectors in the Old 
World [1], while Lutzomyia spp. transmit New World 
leishmaniasis [1,2], 

In Thailand, P. stantoni and S. hodgsoni are commonly 
found in cave dwellers [3-7], These two species are 
sympatric in some places, especially on isolated islands in 
Chumphon province in the Southern part of Thailand [8], In 
particular, Phlebotomus and Sergentomyia specimens were 
found in swiftlet caves on Lang Ga Jiew Island of 
Chumphon province, and S. hodgsoni was the prédominant 
sandfly species on isolated islands. The identification of 
sandfly species belonging to different généra is not difficult 
using morphological examination. Our objective was not to 
test for the morphometric récognition of different généra, 
but to assess the quality of species récognition within each 
genus in spite of possible géographie changes. 

To set up the control strategies against leishmaniasis 
transmission, essentially an entomological surveillance 
program, it may be important to improve our knowledge 
about the distribution of vectors, or potential vectors. To 
reach such objective and to ensure regular updates require 
the use of low-cost, fast and accurate tools to identify 
sandflies. There hâve been very few studies on species 
variation of sandflies from different geographical locations, 
especially none in Thailand, using géométrie morphome- 
trics (GM). In this study, we investigated intraspecific and 
interspecific wing size and shape variation of P. stantoni 
and S. hodgsoni in different geographical locations, 
including between the mainland and isolated islands. 

The information derived from this study could help to 
evaluate the morphometric approach as a routine charac- 
terizing tool for future epidemiological studies of sandflies 
in Thailand. 

2. Materials and methods 

2.3. Sandfly collection 

Using Center for Disease Control (CDC) light traps, 
sandflies were collected in 2011 in seven different 


geographical locations across Thailand (Fig. 1). Isolated 
islands and mainland location sites were selected in this 
study. Three isolated islands, including Lang Ga Jiew (LGJ), 
Ka (KA), and Mapraw (MAP) Islands, were selected for 
sandfly collection; these three islands are located in the 
Gulf of Thailand, in Chumphon province, Southern 
Thailand. LGJ Island is a small offshore island that covers 
an area of 0.063 km 1 2 and is located approximately 8 km 
away from the shore. MAP and KA islands are located 
approximately 1 km away from the shore. Four other 
mainland location sites were selected: a small sample from 
an urban location in Nonthaburi province (NTB), which is 
located in central Thailand, and other from limestone caves 
of Kanchanaburi province (KCB, western Thailand), Rat- 
chaburi province (RCB, central Thailand) and mountainous 
areas of Kamphaeng Phet province (KPP, northern 
Thailand). 

2.2. Sandfly préparation 

Collected female sandflies were stored in 95% ethyl 
alcohol and mounted in Hoyer medium on glass micro- 
scopie slides. Morphological identification was based 
mainly on the keys of Lewis [9,10], Subsequently, 
P. stantoni and S. hodgsoni were photographed using a 
Nikon DS-Ril SIGHT digital caméra connected to stereo- 
microscope Nikon AZ 100 M (Nikon Corp., Tokyo, Japan) 
with a size scale provided on the picture. A total of 
157 right wing pictures from P. stantoni and S. hodgsoni 
across different locations were processed for landmark- 
based GM analysis (Table 1 ). 

2.3. Landmark-based GM analysis 

Twelve landmarks were digitized (Fig. 2). When 
connected after a Procrustes superimposition onto the 
consensus configuration, they appear as polygons allowing 
an easier Visual comparison of mean wing shape between 
sandfly species and populations. The Procrustes superpo- 
sition was conducted according to the Generalized 
Procrustes Superimposition (GPA) procedure [11,12], 
generating both size (centroid size [CS]) and shape (partial 
warps [PW]) variables. 

CS is defined as the square root of the sum of the 
squared distances between the center of the configuration 
of landmarks and each individual landmark [13]. Size 
comparisons were performed using a non-parametric 
permutation test (1000 cycles), and comparisons were 
illustrated by quantile boxes. 
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Fig. 1 . Map of sandfly collection sites in this study. 


The shape variables (partial warps scores [PW]) were 
processed by standard muitivariate analysis to statistically 
compare groups. The principal components (PCA) of the 
PW (i.e„ the relative warps [RW]) were computed to 
visualize the morphospace for sandfly species from 


different locations. The discriminant analysis (DA) used 
as input variables the RW, or a few first of them according 
to the sample sizes, and was illustrated by the correspond- 
ing factor maps. For these DA, the two localities with very 
low sample sizes were removed: NTB (6 individuals) and 
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Table 1 

Number of female Phlebotomus stantoni and Sergentomyia h odgsoni wing pictures used for géométrie morphometric analysis. 


Species 

Locations 

Symbols 

Area description 

Number of 
wing pictures 

P. stantoni 

Lang Ga Jiew Island, Chumphon province 

LGJ 

Limestone swiftlet cave on an isolated island 

21 

P. stantoni 

Nonthaburi province 

NTB 

Urban area in central Thailand 

6 

P. stantoni 

Kamphaeng Phet province 

KPP 

Mountainous areas in northern Thailand 

17 

S. hodgsoni 

Ka Island, Chumphon province 

KA 

Limestone swiftlet cave on an isolated island 

7 

S. hodgsoni 

Lang Ga Jiew Island, Chumphon province 

LGJ 

Limestone swiftlet cave on an isolated island 

57 

S. hodgsoni 

Mapraw Island, Chumphon province 

MAP 

Limestone swiftlet cave on an isolated island 

13 

S. hodgsoni 

Kanchanaburi province 

KCB 

Limestone cave in western Thailand 

14 

S. hodgsoni 

Ratchaburi province 

RCB 

Limestone cave in central Thailand 

22 

Total 




157 



Fig. 2. Position of 1 2 landmarks digitized for female Phlebotomus stantoni and Sergentomyia hodgsoni. 


KA (7 individuals). Validated reclassification scores were 
estimated after a DA on the two taxa, after the DA on the 
géographie populations within each species. In this 
validated reclassification, each individual was allocated 
to its closest group (Mahalanobis distance) without being 
used to help détermine a group center [14], The statistical 
significance (P < 0.05) of shape variation was analyzed 
using non-parametric methods (1000 runs) applied to the 
pairwise Mahalanobis distances. A neighbor-joining tree 
was generated based on Procrustes distances between 
sandflies from different locations. 

The allometric effect on the discrimination between 
collecting sites was explored by regressing the discrimi- 
nant factors on centroid size and computing the coefficient 
of détermination. 

2.4. Software 

CLIC version 97 [15,16], freely available at http : / / 
mome-clic.com, was used to conduct the digitization of the 


images and subséquent multivariate analyses on the 
coordinate data. The following modules were used: COO 
for landmarks collection, MOG to generate centroid size 
and shape variables, as well as to compute Procrustes 
distances and VAR to compare the size between groups. 
The Procrustes distances were used to compute the 
neighbor-joining tree using the PHYL1P neighbor module 
[17], and the resulting tree was visualized using the 
NJPLOT software [18j. The module PAD of CLIC was used to 
perform discriminant analyzes between groups, as well as 
to compute Mahalanobis distances and the corresponding 
statistical significance. 

3. Results 

3.3. Wing size variation 

Mean wing CS comparison of P. stantoni and S. hodgsoni 
from different geographical locations (Table 1) showed 
size variation (Fig. 3, Table 2). Mean wing CS was 
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Fig. 3. Variation of the wing centroid size (in mm) of female Phlebotomus stantoni and Sergentomyia hodgsoni front different localities. Vertical bars 
represent individuals. P 25%, percentiles 25%, P 75%, percentiles 75%. Boxes are quantile boxes between P 25% and P 75%, with indication of the médian 
position. 


• P. stantoni 



Fig. 4. Polygons as connected rnean landmark positions to visualize wing shape différences between female Phlebotomus stantoni (black) and 
Sergentomyia hodgsoni (gray) after Procrustes superposition. No amplification applied. 


S. hodgsoni RCB 
S. hodgsoni KCB 
S. hodgsoni MAP 
S. hodgsoni LGJ 
S. hodgsoni KA 
P. stantoni KPP 
P. stantoni NTB 
P. stantoni LGJ 

mm. 


significantly (P <0.05) larger in P. stantoni than in 
S. hodgsoni (1.22 mm and 0.91mm, respectively). 
P. stantoni from mainland in Kamphaeng Phet province 
had the largest wings (mean, 1.28 mm), significantly 
different from P. stantoni from Lang Ga Jiew Island and 
an urban location in Nonthaburi province (P <0.05). 
However, the mean wing CS of P. stantoni from caves on 
Lang Ga Jiew Island and Nonthaburi province were not 
significantly different (P > 0.05). The wing CS of S. hodgsoni 
from a mainland cave in Ratchaburi province was the 
largest (mean, 1.02 mm), whereas S. hodgsoni from Lang Ga 
Jiew Island had the smallest wing CS (0.89 mm). Almost ali 


populations of S. hodgsoni were significantly different 
(P <0.05) from each other. However, there was no 
significant (P >0.05) différence of wing size between 
S. hodgsoni from Ka Island and Lang Ga Jiew Island, Ka 
Island and Kanchanaburi province, or Mapraw Island and 
Kanchanaburi province. 

3.2. Wing shape variation 

After superimposition of the mean landmark configu- 
rations onto the consensus one, P. stantoni and S. hodgsoni 
wing shapes were clearly different (Fig. 4). Polygons of 
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Fig. 5. Polygons as connected mean landmark positions to visualize the wing shape différences in female Phlebotomus stantoni from various geographical 
locations, as obtained after Procrustes superposition. No amplification applied. 


0 S. hodgsoni KA 
9 S. hodgsoni LGJ 



Fig. 6. Polygons as connected mean landmark positions to visualize wing shape différences of female Sergentomyia hodgsoni from various geographical 
locations, as obtained after Procrustes superposition. No amplification applied. 


mean wing shape of within species samples are shown in 
Fig. 5 (P. stantoni ) and Fig. 6 (S. hodgsoni ). Consistent 
intraspecific variation was visible in both sandfly species. 

In agreement with the polygon superimpositions, the 
RW-based morphospace showed no overlapping of wing 
shapes between P. stantoni and S. hodgsoni, as well as 
between some géographie localities within S. hodgsoni 
(Fig. 7). Within each species, significant géographie 
différences (P <0.05) could be found (Table 3). 
P. stantoni from Lang Gajiew Island and Kamphaeng Phet 
province were significantly different (P <0.05), as were 
P. stantoni from Nonthaburi province and Kamphaeng Phet 
province (P< 0.05). 

The factor map of the DA clearly showed the large shape 
divergence between taxa, as well as some non-overlapping 


groups within both species: P. stantoni was almost 
completely subdivided into island and mainland samples 
(Fig. 8), while in S. hodgsoni the trend was similar, showing 
divergence between mainland (KCB, RCB) and island 
samples (LGJ, MAP), and, within these areas, between 
sampling sites (KCB versus RCB, and LGJ versus MAP) 
(Fig. 9). Statistical significance was frequently found, 
except with localities NTB and KA, represented by only a 
few individuals (Table 3). For S. hodgsoni the contribution 
of size to shape-based discrimination ranged from 0.099 
(with DF1) to 0.165 (with DF2), and for P. stantoni (a single 
DF since two populations were compared), it was 0.117. 

The validated reclassification scores (Table 4) con- 
firmed the complété séparation of the two taxa, providing 
an average correct assignment of 79% between two 
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Fig. 7. Wing shape morphospace based on RW variables of female Phlebotomus stantoni (1, 2, 3) and Sergentomyia hodgsoni (4, 5, 6, 7, 8) samples, 
represented by convex hulls. The horizontal axis is the first relative warp (RW1), and the vertical axis is the second relative warp (RW2). 1,P. stantoni LGJ; 2, 
P. stantoni NTB; 3, P. stantoni KPP; 4, S. hodgsoni KA; 5, S. hodgsoni LGJ; 6, S. hodgsoni MAP; 7, S. hodgsoni KCB; 8, S. hodgsoni RCB. 



Black = Croup 1 21; 736410000000 
White = Group 2 17: 000101452310 

Fig. 8. Factor map of the discriminant factor separating LGJ island (black 
bars) and KPP northern mainland (grey bars) samples of female 
Phlebotomus stantoni, as derived from shape variables (16 first RW). 
The NTB sample has been removed from this analysis due to low sample 
size (6 individuals). 

samples in P. stantoni (scores ranging from 71 to 88%), and 
of 65% between four samples in S. hodgsoni (scores ranging 
from 63 to 71%). 

Based on Procrustes distances, the neighbor-joining 
tree first separated the two généra and then, within 
species, continental and islands locations (Fig. 10). 


4. Discussion 

This is the first GM study on sandflies from Thailand, 
here represented by various mainland and island géo- 
graphie locations of two taxa: P. stantoni and 

5. hodgsoni. We selected the simplest and fastest 
géométrie approach, i.e. the landmark-based approach. 
The wings of sandflies are hairy, which could represent a 
problem in accurately identifying some landmarks 

[19] . However, in our samples, manipulating the samples 
for dissection and mounting was generally enough to clear 
most of the hairy cover. Our results suggest that with no 
further changes in the entomological practice of dissect- 
ing and mounting sandflies, relevant information was 
provided with only 1 2 landmarks. We used the right wing 
of the specimens disregarding possible asymmetric 
effects. However, this source of variation has been 
estimated at a 1 or 2% of the interindividual variation 

[20] , which should not interfère with our comparisons 
based on one side of the insect. 

Size and shape provided taxonomie information of 
unequal values. The geometry of the wing venation 
distinguished both taxa without any possible confusion. 
This was in agreement with the likely high evolutionary 
divergence between two généra. The same géométrie 
shape could also allow one to distinguish the géographie 
location within each species. P. stantoni from two 
geographical locations were almost completely distinct. 
S. hodgsoni from mainland and island, as well as within 
these areas, showed on average diverging géométrie 
shapes, although with some overlapping. Discriminating 
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Fig. 9. Factor map of the two first discriminant factors (DF) derived from shape variables of Sergentomyia hodgsoni originating from various géographie 
locations: 1, LGJ island; 2, MAP island; 3, KCB western mainland; 4, RCB central mainland. Each point represents an individuai. The horizontal axis is DF1; 
the vertical axis is DF2. The 1<A sample iras been removed from this analysis due to low sample size (seven individuals). 


between géographie localities was weakly under the 
influence of size variation, suggesting true shape variation. 

The shape divergence between intraspecific sites has 
two possible origins: the environmental effects, or the 
genetic drift due to géographie isolation. There is a growing 
consensus claiming that, even if the environmental effect is 
probably acting in géométrie shape différentiation, genetic 
drift is likely to be the primary force affecting it [21 -26]. An 
argument supporting that hypothesis in our sample is that, 
in spite of significant and many times non-overlapping 
différence between géographie locations, there was no 
interférence of this shape variation with the interspecific 
shape variation. Similar observation was already eviden- 
ced in a heterogeneous species sample of neotropical 
sandflies from various altitudes [27], and hâve been 
published for various other arthropods [23-26,28-30]. 

The size of the wings was computed as a global 
estimation, the centroid size, allowing one to compare it 
among samples based on a global, single value. Because 
wing size indicates specimen size [24], we can infer that 
P. stantoni specimens were generally larger ones than 


S. hodgsoni. More interestingly, in both taxa, specimens 
from the mainland were significantly larger than those 
from island populations. Alternatively, it can be said that in 
both species, northern populations showed larger sizes 
than Southern ones. Such observation has been already 
made in sandflies: the New World Lutzomyia ayrozcii and 
L geniculata showed larger size in higher altitudes [27], the 
Old World P. papatasi also showed size différence in 
relation to different altitudes in Turkey [31 ]; and the same 
species from the northern Atlas Mountains of Morocco was 
larger than that from the Southern Atlas Mountains 
[32]. This size variation may be related to environmental 
factors including température, host population, and other 
ecological or climatic effects [32], and appears to be in 
agreement with the Bergmann’s rule [33], This “rule” 
suffers some exceptions: for instance in our P. stantoni 
sample from an urban area, the specimens were smaller 
than those from the northern areas and island areas. 

As such, the GM approach provides a powerful 
characterizing feature, the géométrie shape, and repre- 
sents a low-cost and useful tool to identify species and 
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Fig. 10. Neighbor-joining tree based on Procrustes distances, with indication of the patristic distances values, for female Phlebotomus stantoni and 
Sergentomyia hodgsoni from different geographical locations. 


Table 2 

Mean wing centroid size of female Phlebotomus stantoni and Sergentomyia hodgsoni. 


Species 

Locations 

N 

Mean (Min-Max) (mm.) 

S.D. 

S.E. 

P. stantoni 

LGJ 

21 

1.19 (0.99-1.32) 

0.09 

0.02 

P. stantoni 

NTB 

6 

1.12 (0.94-1.23) 

0.10 

0.04 

P. stantoni 

KPP 

17 

1.28 (1.15-1.38) 

0.07 

0.02 

S. hodgsoni 

KA 

7 

0.91 (0.84-0.98) 

0.06 

0.02 

S. hodgsoni 

LGJ 

57 

0.89 (0.76-0.98) 

0.05 

0.007 

S. hodgsoni 

MAP 

13 

0.96 (0.90-1.04) 

0.04 

0.01 

S. hodgsoni 

KCB 

14 

0.94 (0.83-0.93) 

0.05 

0.01 

S. hodgsoni 

RCB 

22 

1.02 (0.92-1.16) 

0.07 

0.01 

Min: minimum; Max: 
correspond to Table 1. 

maximum; mm: millimeter; 

Mean: average centroid 

size; S.D.: standard déviation; S.E.: 

standard error. 

The location symbols 


Table 3 

Mahalanobis distances between wing shapes of female Phlebotomus stantoni and Sergentomyia hodgsoni. 

P. stantoni 




S. hodgsoni 







LGJ 

NTB 

KPP 


KA 

LGJ 

MAP 

KCB 

RCB 

LGJ 

0.00 



KA 

0.00 





NTB 

3.78 

0.00 


LGJ 

2.37 

0.00 




KPP 

3.73 a 

5.65 a 

0.00 

MAP 

3.03 

2.29 a 

0.00 







KCB 

4.46 a 

00 

4.91 a 

0.00 






RCB 

4.72 a 

3.67 a 

4.32 a 

3.04 a 

0.00 


a Indicates statistical signifkance at P < 0.05 after a Bonferroni test. Distances and significance of populations with either NTB (six individuals) or KA (seven 
individuals) are indicative, and should be confirmed with larger samples. 
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Table 4 

Validated reclassification scores between two species, 
Phlebotomus stantoni and Sergentomyia hodgsoni, between LGJ and KPP 
within P. stantoni, and between LGJ, MAP, KCB and RCB within S. hodgsoni. 


Species 

Assigned 

Observed 

Percent (%) 

P. stantoni 

38 

38 

100 

S. hodgsoni 

106 

106 

100 

P. stantoni 

LGJ 

15 

21 

71 

KPP 

15 

17 

88 

Total 

30 

38 

79 

S. hodgsoni 

LGJ 

36 

57 

63 

MAP 

9 

13 

69 

KCB 

10 

14 

71 

RCB 

14 

22 

64 

Total 

69 

106 

65 


géographie populations. It could become a routine 
complément to morphological studies for future epidemi- 
ological investigation of leishmaniasis vectors in Thailand. 
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PAPUA NIUGINI biodiversity expédition 


The unusual clingfish Protogobiesox asymmetricus n. gen, n. sp. is described on the basis of 
four specimens collected in deep water off the north coast of Papua New Guinea in 
2012. The species is characterized by its 9-10 dorsal rays, 8 anal rays, 17-24 pectoral-fin 
rays, 15 principal caudal-fin rays, 3 gills, third arch with 3 gill rakers, 34-35 total vertebrae, 
with asymmetrical latéral bending starting behind the skull, bent at an angle of 85°-92°; 
skull asymmetrical in frontal view; skin naked, surface of head and body without striae; 
dise without adhesive papillae. A new subfamily Protogobiesocinae is described for this 
species and Lepadicyathus mendeleevi Prokofiev, 2005, which is redescribed. The new 
subfamily is compared within the family; keys to the subfamilies of Gobiesocidae and the 
species within the new subfamily are presented; its phylogenetic relationship to other 
gobiesocids is inferred based on a multi-locus DNA dataset. 

© 2016 Académie des sciences. Published by Elsevier Masson SAS. Ail rights reserved. 
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Le porte-écuelle insolite Protogobiesox asymmetricus n. gen, n. sp. est décrit à partir de 
quatre spécimens collectés en 2012 dans les eaux profondes au large de la côte nord de la 
Papouasie-Nouvelle-Guinée. L’espèce se caractérise par sa nageoire dorsale avec 9- 
10 rayons, sa nageoire anale avec 8 rayons, sa nageoire pectorale avec 17-24 rayons, 
15 rayons principaux sur la nageoire caudale, 3 branchies, un troisième arc branchial avec 
3 branchiospines, une colonne vertébrale avec un total de 34-35 vertèbres, le corps arrière 
du crâne présentant une pré-asymétrie latérale et plié à un angle de 85°-92, un crâne 
asymétrique en vue de face, une peau nue, une surface de tête et un corps sans stries, un 
disque sans papille adhésive. Une nouvelle sous-famille, les Protogobiesocinae, est décrite 
pour cette espèce, et Lepadicyathus mendeleevi Prokofiev 2005 est décrit à nouveau. La 
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nouvelle sous-famille est comparée aux autres membres de la famille; des clés pour les 
sous-familles de Gobiesocidae et les espèces incluses dans cette nouvelle sous-famille 
sont présentées ; sa relation de parenté au sein des gobiesocidés est inférée à partir d'un 
jeu de données ADN multi-locus. 

© 2016 Académie des sciences. Publié par Elsevier Masson SAS. Tous droits réservés. 


1. Introduction 

The clingfishes of the family Gobiesocidae are predom- 
inantly marine fishes distributed worldwide in tropical 
and temperate seas, with only eight species of the genus 
Gobiesox Lacepède 1800 living in freshwater streams of the 
tropics. They occur on hard substrata, usually on rocky 
bottom or in coral reefs, mostly in shallow waters, but 
some live on the continental or insular shelf down to 337 m 
( Kopua nuimata Hardy 1984 from New Zealand; see Hardy 
[1]; Hutchins [2]). Clingfishes are small in body size 
(usually less than 50 mm in total length), characterized by 
an adhesive dise formed by the pelvic fins, the head 
depressed, the skin naked, one dorsal and anal fin each, and 
several specialized osteological characters. The family was 
revised by Briggs [3], who distinguished 30 généra and 
85 species, with some divided in several subspecies. 
Springer & Fraser [4] examined the osteology of the species 
in the enigmatic perciform family Cheilobranchidae, and 
assigned its members (arranged in a single genus Alabes 
Cloquet 1816) to the Gobiesocidae. The following généra 
and species were described subsequently to the révisions 
of Briggs [3] and Springer & Fraser [4]: Lecanogaster, 
L. chrysea and Opeatogenys cadenati by Briggs [5], 
Aspasmodes andÆ briggsi by Smith [6], Gobiesox marijeanae 
by Briggs [7], Gobiesox fluviatilis and G. mexicanus by Briggs 
& Miller [8], Lepadichthys bolini by Briggs [9], Gobiesox 
lucayanus by Briggs [10], Lepadichthys ctenion and 

L. erythraeus by Briggs & Link [11], Diplecogaster bimaculata 
euxinica by Murgoci [12], Apletodon knysnaensis by Smith 
[13] (synonym of A. pellegrini according to Briggs [14]), 
Lepadichthys lineatus by Briggs [15], Lissonanchus and 
L lusheri by Smith [16], Tomicodon prodomus by Briggs [17], 
Derilissus and D. nanus by Briggs [18], Lepadichthys caritus 
by Briggs [19], Arcos decoris, Rimicola brevis and Tomicodon 
bidens by Briggs [20], Tomicodon rhabdotus by Smith-Vaniz 
[21 ], Gymnoscyphus and G. ascitus by Bôhlke & Robins [22], 
Derilissus kremnobates and D. vittiger by Fraser [23], 
Derilissus altifrons by Smith-Vaniz [24], Discotrema and 
D. crinophilum by Briggs [25] as D. crinophila ), Modicus, 

M. minimus and M. tangaroa by Hardy [26], Cochleoceps 
bassensis by Hutchins [27], Propherallodus, P. briggsi, 
PheraUodichthys and P. meshimaensis by Shiogaki & Dotsu 
[28], Kopua and K. nuimata by Hardy [29], Aspasmogaster 
occidentalis by Hutchins [30], Tomicodon abuelorum by 
Szelistowski [31], Kopua kuiteri by Hutchins [2], Cochleo- 
ceps viridis, C. bicolor and C. orientalis by Hutchins [32], 
Posidonichthys and P. hutchinsi by Briggs [33], Gobiesox 
juniperoserrai by Espinosa et al. [34], Apletodon incognitus 
by Hofrichter & Patzner [35], Lepadichthys springeri by 
Briggs [36], Tomicodon reitzae by Briggs [37], Rimicola 
cabrilloi by Briggs [38], Tomicodon briggsi, T. clarkei, 


T. cryptus, T. lavettsmithi and T. leurodiscus by Williams & 
Tyler [39], Alabes elongata, A. gibbosa, A. obtusirostris, 
A. occidentalis and A. scotti by Hutchins & Morrison [40], 
Acyrtus pauciradiatus by Sampaio et al. [41 ], Lepadicyathus 
and L mendeleevi by Prokofiev [42], Alabes bathys and 
A. springeri by Hutchins [43], Apletodon wirtzi by Fricke 
[44], Discotrema monogrammum and D. zonatum by Craig & 
Randall [45], Briggsia and B. hastingsi by Craig & Randall 
[46], Apletodon barbatus by Fricke et al. [47 ], Aspasmichthys 
alorensis and Lepadichthys akiko by Allen & Erdmann [48], 
Kopua japonica by Moore et al. [49], Derilissus lombardii by 
Sparks & Gruber [50], Acyrtus lanthanum by Conway et al. 
[51 ], Unguitrema and U. nigrum by Fricke [52], Diplecogaster 
tonstricula by Fricke et al. [53], Kopua vermiculata by 
Shinohara & Katayama [54], and Lepadichthys bilineatus by 
Craig, Bogorodsky & Randall in Craig et al. [55], 

At présent, 166 species belonging to 48 généra of 
gobiesocids are recognized; they are currently classified 
into nine subfamilies (Briggs [3]; Conway et al. [56]): 
Aspasminae Briggs 1955, Cheilobranchinae Günther 1870, 
Chorisochisminae Briggs 1955, Diademichtyinae Whitley 
1950, Diplocrepinae Briggs 1955, Gobiesocinae Bleeker 
1859, Haplocyclicinae Briggs 1955, Lepadogastrinae 
Canestrini 1871-1872, Trachelochisminae Briggs 1955. A 
list of a generic classification within these subfamilies was 
provided by Conway et al. [56], 

During the PAPUA NIUGIN1 Biodiversity Expédition in 
2012, specimens of an unusual, laterally asymmetrical 
clingfish were collected in deep water off the north-eastern 
coast of Papua New Guinea. They were compared with the 
other known gobiesocid species and other taxa in 
gobiesocid-related fish families, and were found to 
represent an undescribed subfamily, genus and species 
of gobiesocid fishes, which is described in the présent 
paper. In addition, the genus Lepadicyathus Prokofiev 2005, 
which is currently classified in the Aspasminae, is revised; 
the only known species, Lepadicyathus mendeleevi Pro- 
kofiev 2005, is redescribed. 

2. Methods and materials 

Our descriptive methods follow Briggs [3] and Fricke 
et al. [57]. Sex détermination based on the presence or the 
absence of urogénital papillae is confirmed by internai 
examination (gonads). The abbreviation “SL” refers to the 
standard length (measured from the tip of the snout to the 
middle of the caudal fin base), “TL” to the total length 
(measured from the tip of the snout to the end of the caudal 
fin). The adhesive dise is divided into three different areas: 
région A is the anterior portion, région B is the posterior 
one, and région C is the centre of the dise (as illustrated by 
Briggs [3]). In the description of new taxa, data of the 
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holotype are given first, followed by data of the paratypes 
in parenthèses. Fin rays are counted using the method of 
Fricke (1983), where spines are expressed as Roman 
numerals, unbranched soft rays are expressed as lower- 
case Roman numerals, and branched rays as Arabie 
numerals. Authors and dates of family-group taxa follow 
Laan et al. [58]. The updated key to subfamilies of 
Gobiesocidae is based on Briggs [3 ]. To elucidate the 
relationships of the unusual clingfish discovered in this 
study and to assign it to its proper higher hiérarchie 
taxonomie rank, its multi-locus DNA sequences were 
generated and compared to those obtained from the 
représentative taxa in the Gobiesocidae and other gobie- 
socid-related families determined by previous studies 
investigating acanthomorph phylogeny (Chen et al. 
[59,60]; Dettaï & Lecointre [61]; Near et al. [62]). The list 
of taxa used for the molecular analysis in this study was 


provided in Table 1. The collection of new molecular data 
was carried out by the following procedures. Small pièces 
of muscle or fin tissue were excised from the specimens, 
preserved in 95% éthanol, and stored at -20 °C in the 
Marine Biodiversity and Phylogenomics Laboratory at the 
Institute of Oceanography, National Taiwan University 
(NTU), Taipei. Genomic DNA was extracted using an 
automated DNA-extractor 199 (LabTurbo 48 Compact 
System with LGD 480-220 kits: Taigene 200 Bioscience 
Corporation, Taipei) following the manufacturées 201 pro- 
tocol. The polymerase chain reaction (PCR) was used to 
amplify the following targeting gene fragments: recombi- 
nation activating gene 1 (RAG1), rhodopsin (RH), early 
growth response 2B (EGR2B), cytochrome oxidase 
subunit I (COI), 12S and 16S ribosomal RNAs (12S and 
1 6S). The first three are encoded in the nuclear genome and 
the remaining three are part of the mitochondrial genome. 


Table 1 

Taxa, genes, and Genbank accession numbers for the gene sequences of représentative species. 


Taxon 

Tissue ID 

Voucher 

GenBank accession no. 

RAG1 Rhodopsin 

EGR2B 

COI 

12S 

16S 

Pleuronectidae 1 

Hippoglossus stenolepis Schmidt 1 904 

WJC148 

NA 

KF311999 

KF312141 

KF312079 

AM749128 

AM749128 

AM749128 

Carangidae 1 

Parastromateus niger (Bloch 1795) 

WJC11 

NA 

EF095654 

EF095616 

KC442131 

KC442079 

EF095562 

EF095590 

Embiotocidae 

Ditrema temminckii Bleeker 1853 

WJC171 

NA 

KY1 26055 

KY1 26046 

KY1 26033 

AP009129 

AP009129 

AP009129 

Mugilidae 

Liza aurata (Risso 1810) 

WJC18 

NA 

KF017112 

KF017144 

KF017049 

JQ060457 

EU7 15447 

GQ252691 

Pomacentridae 

Dascyllus aruanus (Linnaeus 1758) 

WJC377 

NA 

EF095674 

EF095632 

KY1 26034 

KJ968020 

AF081228 

AF1 19402 

Cichlidae 

Astronotus ocellatus (Agassiz 1831) 

WJC168 

NA 

EF095671 

EF095629 

JN231060 

KC442077 

EF095575 

EF095603 

Ambassidae 

Parambassis wolffii (Bleeker 1850) 

WJC49 

NA 

EF095647 

EF095612 

KY1 26035 

EF095558 

EF095586 


Bedotiidae 

Bedotia geayi Pellegrin 1907 

WJC399 

NA 

EF095640 

AY141267 

KC442117 

AY290799 

AY141339 

AY141409 

Adrianichthyidae 

Oryzias latipes (Temminck & Schlegel 1846) 

NA 

NA 

EF095641 

AB001606 


AP00442 1 

EF095555 

EF095583 

Clinidae 

Heterostichus rostratus Girard 1854 

WJC118 

SIO 01-179 

HQ1 68866 

HQ1 68988 

KY1 26036 

HQ1 68631 

HQ1 58801 

AY822076 

Labrisomidae 

Labrisomus nuchipinnis 

WJC86 

NA 

KY1 26056 

KY1 26047 

KY1 26037 

KF930015 

AY098808 

AY098848 

(Quoy & Gaimard 1824) 

Gobiesocidae 

Aspasminae 

Aspasma minima (Dôderlein 1887) 

Cheilobranchinae 

Alabes hoesei Springer & Fraser 1976 

WJC5173 

1.45630-050 

KY1 26057 

EU637943 

KY1 26038 

NC_008130 

KY1 26065 

NC_008130 

KY1 26072 

NC_008130 

KY126080 

Alabes scotti Hutchins & Morrison 2004 

WJC5174 

1.46305-001 

KY1 26058 

KY1 26048 

KY1 26039 

KY1 26066 

KY1 26073 

KY1 26081 

Diademichthyinae 

Lepadichthys lineatus Briggs 1966 

WJC5008 

14-1309 Dili 

KY1 26059 

KY1 26049 

KY1 26040 

KY1 26067 

KY1 26074 

KY126082 

Gobiesocinae 

Acyrtus lanthanum Conway, 

WJC5006 

01-033-Arsp 

KY1 26060 

KY1 26050 

KY1 26041 

KY1 26068 

KY1 26075 

KY126083 

Baldwin & White 2014 

Gobiesox strumosus Cope 1870 

WJC65 

NA 

KY1 26061 

KY126051 

KY1 26042 

KY1 26069 

KY1 26076 

KY126084 

Acyrtops beryllinus 

WJC82 

NA 

KY1 26062 

KY1 26052 

KY1 26043 

KY1 26070 

KY1 26077 

KY126085 

(Hildebrand & Ginsburg 1927) 
Tomicodon humeralis (Gilbert 1890) 

Lepadogastrinae 

Apletodon dentatus (Facciolà 1887) 

WJC305 

SIO 02-1-1 

NA 

HQ1 68876 

KY1 26063 

HQ1 69001 

KY1 26053 

KY1 26044 

HQ1 68639 

KJ616456 

KY1 26078 

KY126086 

Lepadogaster lepadogaster 
(Bonnaterre 1788) 

Protogobiesocinae 

Protogobiesox asymmetricus n. sp. 

PNG1169 

NTUM 10633 

KY1 26064 

AY141273 

KY1 26054 

KY1 26045 

KF369136 

KY1 26071 

AY141347 

KY1 26079 

AY141417 

KY126087 


a Outgroups used for rooting the inferred phylogenetic tree. 
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Table 2 

Primers used in this study. 


Locus/Primer 1 


Primer sequences (5'-3') 

Source 

RGA1 

RI 2533F 

CTGAGCTGCAGTCAGTACCATAAGATGT 

Lôpez et al., 2004 


RI 4078R 

TGAGCCTCCATGAACTTCTGAAGRTAYTT 

Lôpez et ai, 2004 


RI 4061 R 

AATACTTGGAGGTGTAGAGCCAGT 

Chen et al., 2007 


RI 4090R 

CTGAGTCCTTGTGAGCTTCCATRAAYTT 

Lôpez et ai, 2004 

EGR2B 

E2B 252F 

CGCAACCAGACTTTCACCTAY 

Chen et ai, 2013 


E2B 261F 

TTCACCTAYATGGGNAAGTTCTCMAT 

Chen et ai, 2013 


E2B 270F 

ATGGGRAAGTTCTCCATCGAC 

Chen et ai, 2013 


E2B 278F 

AGTTTTCCATCGACTCSCAGTA 

Chen et ai, 2008 


E2B 287F 

TTGACTCSCAGTATCCAGGTAAC 

Chen et ai, 2008 


E2B 10 78 R 

AATTTGCGNCCGCAGSAGTC 

Chen et ai, 2013 


E2B 1 078R-bis 

GAACTTACGNCCGCAGAARTC 

Chen et ai, 2013 


E2B 1108R 

TmGTGTGTCTCTTrCTYTCGTC 

Chen et ai, 2008 


E2B 1 1 12R 

ATTTTNGTGTGTCGYTTYCTC 

Chen et ai, 2013 


E2B 1117R 

AGGTGGATTTTGGTGTGTCTYTr 

Chen et ai, 2008 


E2B 1121R 

CCTCAGGTGGATTTTAGTGTGTC 

Chen et ai, 2013 

Rhosopsin 

RH 1F 

ATGAACGGCACAGARGGAC 

Chen et ai, 2003 


RH 193F 

CNTATGAATAYCCTCAGTACTACC 

Chen et ai 


RH 1039R 

TGCTTGTTCATGCAGATGTAGA 

Chen et ai 

COI 

CoxLFishFl 

TCAACCAACCACAAAGACATTGGCAC 

Ward et ai, 2005 


CoxI_FishF2 

TCGACTAATCATAAAGATATCGGCAC 

Ward et ai, 2005 


CoxI_FishRl 

TAGACTTCTGGGTGGCCAAAGAATCA 

Ward et ai, 2005 


CoxI_FishR2 

ACTTCAGGGTGACCGAAGAATCAGAA 

Ward et ai, 2005 

12S 

L1090 

AAAGCACGGCACTGAAGATGC 

Palumbi et ai, 1991 


H1478 

TTTCATGTTTCCTTGCGGTAC 

Palumbi et ai, 1991 

16S 

arL 

CGCCTGTTTAACAAAAACAT 

Palumbi et ai, 1991 


brH 

CCGGTCTGAACTCAGATCAGT 

Palumbi et ai, 1991 


RAG1, recombination activation gene 1; EGR2B: early growth response protein gene 2B; COI: mitochondrial cytochrome oxidase subunit 1; 12S, 16S, 
mitochondrial ribosomal RNA; R = A or G, Y = C or T, S = C or G. 
a Reverse primers in italics 


These genes were selected for their ability to provide 
Sound phylogenetic information for inter- and intra- 
familial' relationships in several previous phylogenetic 
studies focusing on the percomorph fishes (e.g., Chen et al. 
[60]; Lo et al. [63]). Primer sequences in this study are 
listed in Table 2. Température cycling profiles for 
amplification consisted of an initial dénaturation stage 
(95 °C, 5 min) followed by 35 cycles, each with a 
dénaturation step (95 °C, 40 s), an annealing step (53 °C 
for RAG1 and EGR2B; 55 °C for the other genes, 40 s) and an 
élongation step (72 °C, 90 s for RAG1 and 60 s for the other 
genes), before a final extension stage (72 °C, 7 min). PCR 
Products were purified using the AMPure magnetic bead 
clean-up protocol (Agencourt Bioscience Corp., USA). 
Purified PCR products were sequenced by Sanger sequenc- 
ing using dye-labelled terminators. Sequence détermina- 
tions from Sanger reaction products were generated on 
ABI 3730 analyzers (Applied Biosystems, USA) at Genomics 
BioSci & Tech (Taipei) and at the Center of Biotechnology 
(National Taiwan University). For each locus and specimen, 
the chromatograms of the sequences were checked and 
edited manually using the program CodonCode Aligner 
3. 7. 2. 2 (Codoncode Corporation, Dedham, MA, USA). 
Chromatograms were examined for evidence of sequenc- 
ing artefacts or problematic base calls. Verified comple- 
mentary chromatograms were then assembled to build 
consensus sequences. Ail sequences from the orthologous 
gene were aligned using the automatic multiple-alignment 
program Muscle (Edgar [64]). The resulting multiple 
sequence alignments were adjusted by eye using the 
inferred arnino acid translation as a guide for inferred gap 


placement in protein coding genes and the loop/stem 
structure of 12S and 16S ribosomal RNAs to identify loop 
régions where site homology was difficult to ascertain with 
confidence. Alignments were edited using Se-Al v2.0all 
(Rambaut [65]). Short régions of the alignments from 
EGR2B, 12S and 16S in which homology assessment is 
uncertain were excluded from the further steps of 
phylogenetic analysis. Phylogenetic analysis was conduc- 
ted based on the combined data matrix complied from the 
aligned DNA sequences from the six genes with a 
partitioned maximum likelihood (ML) method as imple- 
mented in the parallel version of RAxML 7.4.2 (Stamatakis 
[66]). A mixed-model analysis was used that allows the 
independent estimation of individual models of nucléotide 
substitution for each partition. Three partitions (by codon) 
were assigned for cytochrome b and rhodopsin. Fourteen 
partitions (by gene and by codon position of protein coding 
gene) were assigned. Because RAxML only provides GTR- 
related (Yang [67]) models of rate heterogeneity for 
nucléotide data (Stamatakis [66]), the nucléotide substi- 
tution model GTR + r+I was employed for the analysis. 
Optimal ML tree search was conducted with 100 separate 
runs using the default algorithm of the program for each 
run and the final tree with the best ML score was selected 
among the 100 ML trees of these runs. Finally, nodal 
support was assessed with bootstrapping (BS) (Felsenstein 
[68]) with the maximum likelihood (ML) criterion, based 
on 500 pseudo-replicates. The inferred phylogenetic tree 
was rooted by two non-“Ovalentaria” taxa: Hippoglossus 
stenolepis (Pleuronectidae) and Parastromateus niger 
(Carangidae). The specimens cited in the présent paper 
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are deposited in the following collections: AMS, The 
Australian Muséum, Sydney, Australia; MNHN, “Muséum 
national d’histoire naturelle”, Paris, France; NTUM, Na- 
tional Taiwan University Muséums, Taipei, Taiwan; Z1N, 
Laboratory of Ichthyology, Zoological Institute, Russian 
Academy of Sciences, St. Petersburg, Russia. The abbre- 
viations of repositories of additional materials follow 
Fricke & Eschmeyer [69], 

2.1. Comparative material for morphological examination 

Acyrtops amplicirrus : CAS-SU 18101 (2 paratypes), St. 
Thomas, Virgin Islands. Acyrtus artius: CAS-SU 23254 
(holotype), Curaçao, Netherlands West Indies; CAS-SU 
47942 (1 paratype), Curaçao, Netherlands West Indies. 
Alabes dorsalis : SMNS 1694 (4), Murray River estuary, 
South Australia; SMNS 21395 (1), Cape Dromedary, New 
South Wales, Australia. Apletodon barbatus : SMNS 26427 
(holotype), Sâo Tiago Island, Cape Verde Islands; MNHN 
2009-1592 (1 paratype), same data as the holotype; SMNS 
24604 (1), Sal Island, Cape Verde Islands; SMNS 24605 (1 ), 
Sal Island, Cape Verde Islands; SMNS 26428 (19 paratypes), 
same data as the holotype; USNM 396967 (1 paratype), 
same data as the holotype. A. dentatus: CCML uncat. (2), 
Alegranza Island, Canary Islands; CCML uncat. (2), 
Lanzarote Island, Canary Islands; SMNS 12664 (1), Genoa, 
Italy. A. incognitus : NMW 93029 (holotype), Banyuls-sur- 
Mer, France; CCML uncat. (2), Gran Canaria, Canary Islands; 
CCML uncat. ( 1 ), La Graciosa, Canary Islands ; SMNS21814 
(1), Faial, Azores Islands. A. peUegrinv. SAIAB 10255 (1), 
Knysna, South Africa; SAIAB 10852 (11), Knysna, South 
Africa; SAIAB 14599 (7), South Africa; SAIAB 43756 (5), 
False Bay, South Africa; USNM 198169 (1 paratype), 
Knysna, South Africa; USNM 270272 (1), Algoa Bay, South 
Africa. A. wirtzi: SMNS 24130 (holotype), Bombom Island, 
Sâo Tomé and Principe; MNHN 2005-0170 (1 paratype), 
same data as the holotype; SMNS 24446 (5 paratypes), 
same data as the holotype; SMNS 24132 (2 paratypes), 
same data as the holotype; SMNS 25472 (5), Limbe, 
Cameroon; USNM 381374 (1 paratype), same data as the 
holotype. Arcos nudus: SMNS 11301 (4), Rio San Juan, 
Dominican Republic; SMNS 11401 (12), Rio Baonico, 
Dominican Republic; SMNS 11408 (3), Rio Nazaito, 
Dominican Republic. Aspasma minima: NMW 76995 
(2 syntypes), Sagami Bay, Japan. Aspasmichthys ciconiae: 
CAS-SU 7136 (1 syntype), Inland Sea, Japan; USNM 50760 
(1 syntype), Inland Sea, Japan. Aspasmogaster costata : 
SMNS 14754 (9), Red Head, New South Wales, Australia. 
Chorisochismus dentex : USNM 119223 (1 ), Algoa Bay, South 
Africa. Cochleoceps spatula: SMNS 12519 (1), Cervantes 
Island, Western Australia. Conidens laticephalus : SMNS 
24701 (10), Northeast Cape, Taiwan. C. samoensis: SMNS 
21982 (1), Province Sud, Grande Terre, New Caledonia; 
SMNS 22836 (3), Province Sud, Grande Terre, New 
Caledonia. Creocele cardinalis : SMNS 14806 (1), Tasmania, 
Australia. Dellichthys morelandi : SMNS 13754 (5), North 
Island, New Zealand; SMNS 13762 (1), North Island, New 
Zealand; SMNS 13780 (4), North Island, New Zealand; 
SMNS 13827 (1), North Island, New Zealand; SMNS 14142 
(4), North Island, New Zealand. Derüissus altifrons: ANSP 
112690 (holotype), Dominica. D. vittiger: ANSP 109626 


(holotype), Venezuela. Diademichthys lineatus : SMNS 
18156 (2), Fiji, Viti Levu; SMNS 21975 (3), New Caledonia, 
Récif Goéland; SMNS 21977 (1), New Caledonia, Nouméa; 
SMNS 21979 (2), New Caledonia, Baie Maâ; SMNS 21980 
(2), New Caledonia, Bancs Nord; SMNS 22995 (1), New 
Caledonia, Baie de Pritzbue; SMNS 23522 (1), New 
Caledonia, Nouméa; SMNS 23997 (1), New Caledonia, 
Nouméa; SMNS 25414 (1), New Caledonia, Nouméa; SMNS 
26545 (1), New Caledonia, Baie de Goro. Diplecogaster 
bimaculata: HUJ 20577 (7), Balearic Islands, S Mallorca; 
HUJ 20593 (1), Balearic Islands, NE Mallorca; HUJ 20602 
(2), Balearic Islands, NW Mallorca; HUJ 20623 (1), Balearic 
Islands, SSE Mallorca; SMNS 12541 (1), Pyrénées Orien- 
tales, France; SMNS 13177 (1), Giglio Island, Italy; SMNS 
14049 (2), Giglio Island, Italy; SMNS 19061 (2), Karavas 
Alsavcak Bay, Northern Cyprus; SMNS 19204 (2), Giglio 
Island, Italy; SMNS 20163 (8), Caniço de Baixo, Madeira; 
SMNS 20347 (1), Tabarca, Tunisia; SMNS 21202 (2), Porto 
Novo, Madeira. D. ctenocrypta: ZMUC P9037 (holotype), 
Gran Canaria, Canary Islands. D. pectoralis: SMNS 11916 
(4), Faial Island, Azores Islands. D. tonstricula : ZSM 40089 
(holotype), Dakar, Sénégal; CCML uncat. (2 paratypes), 
Fuerteventura, Canary Islands; ZSM uncat. (5 paratypes), 
Dakar, Sénégal. Diplocrepis puniceus : SMNS 10012 (2), 
North Island, New Zealand; SMNS 13755 (2), North Island, 
New Zealand; SMNS 13769 (1), North Island, New Zealand; 
SMNS 13782 (1), North Island, New Zealand; SMNS 13828 
(8), North Island, New Zealand; SMNS 13832 (4), North 
Island, New Zealand; SMNS 13922 (42), South Island, New 
Zealand; SMNS 13936 (32), South Island, New Zealand; 
SMNS 13945 (24), South Island, New Zealand; SMNS 14003 
(1), South Island, New Zealand; SMNS 14043 (1), South 
Island, New Zealand; SMNS 14069 (1), South Island, New 
Zealand. Discotrema crinophilum : SMNS 16677 (1), 
Philippines, Balicasag Island; SMNS 16678 (1), Philippines, 
Balicasag Island; SMNS 17896 (1), Society Islands, Moorea 
(new record); SMNS 21981 (1), New Caledonia, Nouméa. 
D. monogrammum : BPBM 39040 (holotype), Papua New 
Guinea, New Britain; BPBM 36504 (2), Indonesia, Flores. D. 
zonatum: BPBM 38972 (holotype), Fiji, Charybdis Reef. 
Gastrocyathus gracilis: SMNS 10011 (1), North Island, New 
Zealand; SMNS 14124 (1), North Island, New Zealand. 
Gastroscyphus hectoris : SMNS 13763 (1), North Island, New 
Zealand; SMNS 13770 (1), North Island, New Zealand; 
SMNS 13829 (3), North Island, New Zealand; SMNS 13834 
(27), North Island, New Zealand; SMNS 13854 (1), North 
Island, New Zealand; SMNS 13862 (1), North Island, New 
Zealand; SMNS 13923 (26), South Island, New Zealand; 
SMNS 13924 (76), South Island, New Zealand; SMNS 1 3938 
(6), South Island, New Zealand; SMNS 13947(9), South 
Island, New Zealand; SMNS 13948 (17), South Island, New 
Zealand; SMNS 13949 (30), South Island, New Zealand; 
SMNS 14008 (5), South Island, New Zealand; SMNS 14018 
(1), South Island, New Zealand. Gobiesox strumosus : SMNS 
13612 (1), Rio de Janeiro, Brazil. Gouania willdenowi : SMNS 
421 (4), Hvar Island, Croatia; SMNS 576 (2), Hvar Island, 
Croatia; SMNS 13175 (1), Giglio Island, Italy; SMNS 13573 
(1), Costa Brava, Spain; SMNS 13574 (1), Crete Island, 
Greece; SMNS 16808 (6), Sicily, Italy; SMNS 22064 (6), Cres 
Island, Croatia. Gymnoscyphus ascitus : ANSP 113587 
(holotype), St. Vincent. Haplocylix littoreus: CAS-SU 
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47678 (2), North Island, New Zealand. Kopua nuimata : 
SMNS 10013 (1), North Island, New Zealand; SMNS 13737 
(6), North Island, New Zealand; SMNS 13781 (1), North 
Island, New Zealand; SMNS 13826 (22), North Island, New 
Zealand; SMNS 13833 (2), North Island, New Zealand; 
SMNS 13943 (4), South Island, New Zealand; SMNS 13996 
(20), South Island, New Zealand; SMNS 14030 (1), South 
Island, New Zealand; SMNS 14042 (2), Stewart Island, New 
Zealand; SMNS 14070 (1), South Island, New Zealand; 
SMNS 14127 (1), South Island, New Zealand; SMNS 14151 
(2), South Island, New Zealand. Lecanogaster chrysea: 
BMNH 1958.7.3.1 (holotype), Lingo, Ghana. Lepadichthys 
caritus : BPBM 34140 (1 ), Indonesia, Flores. L. coccinotaenia: 
BPBM 21711 (1), South Africa, KwaZulu-Natal. L. ery- 
thraeus : SMNS 22550 (1), Red Sea, Gulf of Aqaba, Egypt, 
14 km north of Nuweiba. Lfrenatus : SMNS 21978 (1), New 
Caledonia, off Nouméa. L. lineatus : SMNS 22573 (4), Red 
Sea, Gulf of Aqaba, Egypt, 141<m north of Nuweiba. L. 
minor: SMNS 15130 (1), Fiji, Viti Levu; SMNS 17822 (2), 
Cook Islands, Aitutaki; SMNS 20909 (2), La Réunion, Les 
Filaos; SMNS 21020 (2), La Réunion, Les Filaos; SMNS 
21161 (1), La Réunion, Les Filaos; SMNS 21178 (2), La 
Réunion, Saint-Leu; SMNS 22980 (1 ), Loyalty Islands, Lifou; 
SMNS 26966 (3), New Caledonia, Nouméa. Lepadichthys 
sp.: SMNS 21976 (1), New Caledonia, Baie Maâ; SMNS 
23962 (15), New Caledonia, Ile Nou; SMNS 26967 (4), New 
Caledonia, Nouméa; SMNS 26968 (1), New Caledonia, 
Nouméa. Lepadogaster candolii: SMNS 579 (1), Hvar Island, 
Croatia; SMNS 993 (1), Trieste, Italy; SMNS 2662 (6), 
Naples, Italy; SMNS 8996 (1 ), Senj, Croatia; SMNS 9416 (1 ), 
Cres Island, Croatia; SMNS 13176 (1), Giglio Island, Italy; 
SMNS 13575 (5), Saronic Gulf, Greece; SMNS 13576 (1), 
Elba Island, Italy; SMNS 13677 (1), Peloponnes, Greece; 
SMNS 13578 (3), Crete Island, Greece; SMNS 13579 (3), 
Muros, Galicia, Spain; SMNS 13580 (1), Saronic Gulf, 
Greece; SMNS 13581 (1), Costa Brava, Spain; SMNS 13582 
(2), Elba Island, Italy; SMNS 13613 (1), Biograd, Croatia; 
SMNS 13614 (1), Sicily, Italy; SMNS 13615 (1), Sicily, Italy; 
SMNS 13616 (2), Corfu Island, Greece; SMNS 13620 (2), 
Corfu Island, Greece; SMNS 14700 (6), Reis Magos, 
Madeira; SMNS 14721 (2), Cres Island, Croatia; SMNS 
1 6020 ( 1 ), Caniçal, Madeira; SMNS 1 6069 ( 1 ), Siciliy, Italy; 
SMNS 25479 (2), Caniço de Baixo, Madeira. L. lepadogaster: 
SMNS 1079 (1), Azores; SMNS 13556 (1), Biograd, Croatia; 
SMNS 13557 (1), Golfo di Trieste, Italy; SMNS 13558 (1), 
Bari, Italy; SMNS 13559 ( 1 ), Karlobag, Croatia; SMNS 13560 
(2), Corfu Island, Greece; SMNS 13561 (5), Elba Island, 
Italy; SMNS 13562 (1), Collioure, Pyrenées-Orientales, 
France; SMNS 13563 (1), Crete Island, Greece; SMNS 13564 
(1), Cap Negro, N Tétouan, Morocco; SMNS 13565 (2), Elba 
Island, Italy; SMNS 13566 (2), Sea of Marmara, Turkey; 
SMNS 13567 (1), istanbul Province, Turkey, Black Sea; 
SMNS 13568 (1), Granada Province, Spain; SMNS 13569 

(1) , Banyuls-sur-Mer, France; SMNS 13570 (2), Cran 
Canaria, Canary Islands; SMNS 13617 (5), Costa Brava, 
Spain; SMNS 13618 (1), Corfu Island, Greece; SMNS 14699 
(4), Reis Magos, Madeira; SMNS 16726 (6), Fuerteventura, 
Canary Islands; SMNS 19070 (1), Karavaç Alsavçak Bay, 
Northern Cyprus; SMNS 22609 (8), La Palma, Canary 
Islands; SMNS 25316 (1), Bodrum, Turkey; SMNS 25481 

(2) , Caniço de Baixo, Madeira. L. purpurea: SMNS 10295 (1), 


Lanzarote, Canary Islands; SMNS 12549 (1), Collioure, 
France; SMNS 13571 (1), Biograd, Croatia; SMNS 13572 (2), 
Saronic Gulf, Greece; SMNS 15286 (2), Cap Béar, Pyrénées 
Orientales, France; SMNS 15315 (2), La Palma, Canary 
Islands; SMNS 15328 (1), La Palma, Canary Islands; SMNS 
16727 (4), Fuerteventura, Canary Islands; SMNS 16766 (2), 
Fuerteventura, Canary Islands; SMNS 22610 (4), La Palma, 
Canary Islands; SMNS 22623 (1 ), La Palma, Canary Islands; 
SMNS 23501 (5), El Hierro, Canary Islands. Liobranchia 
stria : USNM 149910 (holotype), Saipan, Marianas Islands. 
Modicus minimus: BMNH 1982.6.17.55 (1 paratype), North 
Island, New Zealand. Opeatogenys cadenati: MNHN 1959- 
0059 (holotype), Sénégal; MNHN 1959-0060(2 paratypes), 
Sénégal. Parvicrepis parvipinnis: AMS 1.7708 (5 syntypes), 
Sydney, New South Wales, Australia. Pherallodichthys 
meshimaensis : NSMT-P 46752 (holotype), Kyushu, Japan. 
Pherallodiscus funebris: CAS-SU 20 (9 paralectotypes), Gulf 
of California, Mexico. Pherallodus indicus: CAS-SU 47683 

(1) , Raroia, Tuamotu Archipelago. P. smithi : CAS-SU 31349 
(holotype), KwaZulu-Natal, South Africa. Posidonichthys 
hutchinsi: AMS 1.17615-002 (holotype), Fiddler’s Bay, 
South Australia. Propherallodus briggsi : NSMT-P 46749 
(holotype), Kyushu, Japan. Rimicola muscarum : CAS-SU 
3030 (holotype), Monterey Bay, California, USA. Sicyases 
sanguineus: SMNS 390 (2), Chile. Tomicodon humeralis : 
CAS-SU 70 (5 syntypes), Gulf of California, Mexico. 
Trachelochismus melobesia: SMNS 13926 (1), South Island, 
New Zealand; SMNS 13940(1), South Island, New Zealand. 
T. pirmulatus: SMNS 13925(1), South Island, New Zealand; 
SMNS 13937 (1), South Island, New Zealand; SMNS 13939 

(2) , South Island, New Zealand; SMNS 13946 (4), South 
Island, New Zealand; SMNS 13950 (6), South Island, New 
Zealand; SMNS 13951 (34), South Island, New Zealand. 
Unguitrema nigrum: NTUM 10603 (holotype), Madang, 
Papua New Guinea; MNHN uncat. (1 paratype), same data 
as the holotype. 


3. Results 

3.3. Gobiesocid phylogeny 

The combined dataset for inferring gobiesocid phylog- 
eny contains the sequences of RAG1 gene (1476 bp), 
Rhodopsin gene (819 bp), EGR2B gene (855 bp), COI gene 
(633 bp), 32S gene (344 bp) and 16S gene (334 bp) from a 
total of 22 examined taxa that include 11 gobiesocid 
représentatives from five out of their nine currently 
recognized subfamilies (Table 1 ). The total of 4461 nuclé- 
otides comprise 2184 variable sites (49%) and 736 parsi- 
mony-informative sites. The resulting ML tree is shown in 
Fig. 1. 

From the inferred tree, the monophyly of the Gobie- 
socidae as well as its subfamilies Gobiesocinae and 
Cheilobranchinae is confirmed with 100% bootstrap nodal 
support. The Lepadogastrinae is monophyletic, but it 
received only week support. Among the outgroup taxa, 
the blennioids are the most closely related to the 
Gobiesocidae and their sister-group relationship is sup- 
ported by a strong bootstrap value of 99%. Within the 
family, our resuit divides the gobiesocid taxa sampled in 
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Fig. 1. Phylogenetic tree of the Gobiesocidae inferred by the partitioned maximum-likelihood method with the GTR + r+I nucléotide substitution model 
based on the combined dataset. Branch lengths are proportional to the inferred nucléotide substitutions. Numbers at nodes represent bootstrap values in 
percentage. Values below 50% are not shown. An undetermined specimen of Arcos sp. was further determined to represent a recently described species, 
Acyrtus lanthanum Conway, Baldwin & White 2014 by COI sequence similarity. 


this study into two sub-groups (I and II), one containing the 
taxa from the Gobiesocinae and the other with the 
remaining gobiesocids in four subfamilies plus the newly 
discovered clingfïsh from deep water off the north coast of 
Papua New Guinea (see the description below) (Fig. 1). 
According to our phylogenetic resuit, we can induce that 
the new taxon is a gobiesocid of sub-group II, and it 
appears to be an independent gobiesocid lineage or 
subfamily as it is not tightly clustered with any known 
subfamilies of the Gobiesocidae, at least sampled in this 
study (Fig. 1). 

3.2. Protogobiesocinae new subfamily (Figs. 2-6) 

3.2 A. Type genus 

Protogobiesox n. gen. 

3.2.2. Diagnosis 

A subfamily of gobiesocid fishes with a partial latéral 
asymmetry; head slightly asymmetrical when seen from 
the front; vertébral column laterally bent, at an angle of 
40-92°; anus situated on the right-hand side of the body in 
males, and on the left-hand side in females; position of 
guts slightly asymmetrical to strongly asymmetrical; 
gills 3; pelvic fin dise rudimentary, in two parts, at most 


with a few rudimentary papillae basally in région B. Pelvic 
bones présent, pelvic fin forming a double dise; scapula, 
pectoral radial and rays présent; total vertebrae 33-35. 

3.2.3. Description 

Dorsal-fin rays ix-xii; anal-fin rays viii-x; pectoral-fin 
rays xvii-xxiv; caudal-fin rays (iv-v),x-xv,(iv-v). Gills 3; gill 
rakers on 3rd arch 0-3. 

Head latéral line System with 2 pores in the nasal canal, 
1 -2 pores in the postorbital canal, 0-2 pores in the lacrimal 
canal, 0-2 pores in the preopercular canal, and 0-2 pores in 
the mandibular canal. Vertebrae 33-35 (15 + 18-20). 
Pleural ribs 12-13 pairs. 

Head broad, depressed. Head and skull at least slightly 
asymmetrical in frontal view. Gill membranes attached to 
isthmus. Pelvic bones présent. Body bent to the left-hand 
side (males) or right-hand side (females); vertébral 
column bent at angle of 40-92°, bending starts between 
skull and level of anus. Anus situated on right-hand side of 
the body (males) or the left-hand side of the body 
(females); male with a short urogénital papilla. Skin 
naked; surface of head and body may bear striae. 

Dorsal fin origin on the left-hand side of the body 
(males), and on the right-hand side (females); anal fin 
origin on the right-hand side of the body (males) or the 
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left-hand side of the body (females). Dorsal and anal-fin 
bases attached to caudal-fin base by membranes. Pelvic fin 
forming a double dise; scapula, pectoral radial and rays 
présent. Disc membrane inserting at the base of the 12th- 
16th pectoral-fin ray. Disc rudimentary, incomplète, 
completely or mostly without adhesive papillae. 

3.2.4. Etymology 

The name of the new subfamily is based on the type 
genus, Protogobiesox; its meaning is “ancestral gobiesocid," 
referring to the pelvic fins resembling the putative 
ancestral form of gobiesocids. 

3.2.5. Distribution 

The two described species are distributed off the north 
coast of the main island of Papua New Guinea, western 
Pacific Océan (Fig. 7). 

3.2.6. Comparions 

Within the family Gobiesocidae, this subfamily is 
unique in its latéral asymmetry. It further differs from 
the subfamily Cheilobranchinae in the presence of pelvic 
bones and a dise (pelvic bones and dise greatly reduced or 
absent in Cheilobranchinae), the presence of a scapula, 
pectoral radiais and pectoral-fin rays (ail these absent in 
the Cheilobranchinae, and the total vertebrae number 33- 
35 (60-78 in Cheilobranchinae), and from the other 
subfamilies in the presence of three gills (three and one 
half in the Chorisochisminae, Haplocyclicinae, Lepadogas- 
trinae, Trachelochisminae), the gill membranes which are 
attached to the isthmus (free from isthmus in Diplocre- 
pinae, Gobiesocinae, Haplocyclicinae, Trachelochisminae), 
the double dise (dise single in Chorisochisminae, Diade- 
michthyinae, Gobiesocinae, Haplocyclicinae), and the dise 
région B lacking papillae or with very few rudimental 
papillae in 1-2 rows (with 6-10 rows of flattened papillae 
in the Aspasminae). A key to the subfamilies of the 
Gobiesocinae is presented below, as well as a key to the 
species of the Protogobiesocinae n. subfam. 


9a. Body laterally symmetrical, vertébral Aspasminae 

column not bent; dise région B with 6- 
10 rows of flattened papillae (indo-Pacific) 

9b. Body laterally asymmetrical, vertébral Protogobiesocinae 

column bent at an angle of 40-92 ; dise n. subfam 
région B lacking papillae or with very few 
rudimental papillae in 1-2 rows (Western 
Pacific) 


3.2.8. Key to the species of the subfamily Protogobiesocinae 
n. subfam 


la. Dorsal fin with 9-10 rays; anal fin with 8 rays; 
caudal fin with 15 principal rays; vertébral 
column strongly bent at an angle of 85-92 ; 
bending starts behind skull; skin naked, without 
striae 

lb. Dorsal fin with 11-12 rays; anal fin with 

10 rays; caudal fin with 10 principal rays; 
vertébral column bent at an angle of 40-45 ; 
bending starts just before level of anus; skin 
naked, with numerous, closely set, oblique striae 
consisting of minute tubercles 


Protogobiesox 
asymmetricus 
n. gen., n. sp. 


Lepadicyathus 

mendeleevi 


3.3. Protogobiesox new genus (Figs. 2-5) 

3.3.1. Type species 

Protogobiesox asymmetricus n. gen., n. sp. 

3.3.2. Diagnosis 

A genus of protogobiesocine fishes with a partial latéral 
asymmetry; head slightly asymmetrical when seen from 
the front; vertébral column strongly bent towards the left 
in males, towards the right in females, at an angle of 85- 
92 ; bending starts behind the skull; anus situated on the 
right-hand side of the body in males, on the left-hand side 
in females; position of guts strongly asymmetrical; gills 3; 
pelvic fin dise rudimentary, in two parts, without adhesive 
papillae. Dorsal fin with 9-10 rays; anal fin with 8 rays; 
caudal fin with 15 principal rays. Skin naked, without 
striae. 


3.2.7. Key to the subfamilies of the family Gobiesocidae 


la. Pelvic bones présent, pelvic fin forming a 2 
single or double dise; scapula, pectoral radial 
and rays présent; total vertebrae 25-54 
1 b. Pelvic bones and pelvic fin greatly reduced Cheilobranchinae 
or absent; scapula, pectoral radial and rays 
absent; total vertebrae 60-78 
2a. Gills three and one-half 
2b. Gills three 

3a. Gill membranes free from isthmus 
3b. Gill membranes attached to isthmus 
4a. Disc double (Indo-Pacific to Japan) 

4b. Disc single (New Zealand) 

5a. Disc double (Eastern Atlantic and South 
Africa) 

5b. Disc single (South Africa) 

6a. Gill membranes free from isthmus 
6b. Gill membranes attached to isthmus 
7a. Disc double (South Africa to Indo-Pacific; 

New Zealand) 

7b. Disc single (America plus South Africa) 

8a. Disc double (Indo-Pacific) 

8b. Disc single (South Africa to Indo-Pacific) 


3.3.3. Description 

Dorsal-fin rays x (ix-x); anal-fin rays viii (viii); 
pectoral-fin raysxxiv(xvii-xxiv); caudal-fin rays (v),xv,(iv) 
[(iv-v),xv,(iv)]. Gills 3 (3); gill rakers on 3rd arch 3 (3). 

Detailed description see below under Protogobiesox 
asymmetricus n. gen. n. sp. 

3.3.4. Etymology 

The name of the new genus means “ancestral 
gobiesocid;" it refers to the pelvic fins resembling the 
putative ancestral form of gobiesocids. The gender is 
masculine. 

3.3.5. Remarks 

This genus is monotypic, with Protogobiesox asymme- 
tricus n. gen., n. sp. as the only known species. 

3.4. Protogobiesox asymmetricus n. sp. (Figs. 2-5) 

Asymmetrical deep-water clingfish. 


3 
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5 
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Fig. 2. Protogobiesox asymmetricus gen. nov., sp. nov., NTUM 10628, 
holotype, male, 60.0 mm SL, Papua New Guinea, Madang Province, 
northeast of Taviltae. A (upper) Latéral view; B (centre) dorsal view; C 
(lower) ventral view. Bar: 5 mm. 

3.4.1. Type material 

Holotype. NTUM 10628, male, 60.0 mm SL, Papua New 
Guinea, Madang Province, northeast of Taviltae, 4°31’S 


145°31’E, 380-382 m depth, St. CP4035, sunken logs, some 
covered with sponges, W.-J. Chen, R/V Alis, 17 Dec. 2012. 

Paratypes. Four specimens: MNHN 2016-03333, 1 fe- 
male, 59.0 mm SL, same data as the holotype; NTUM 
10633, 1 male, 53.6 mm SL, same data as holotype; NTUM 
10629, 1 male, 56.0 mm SL, Papua New Guinea, Sandaun 
Province, northwest of Aitape, 3°01’S 142°16’E, 400- 
560 m in depth, St. CP4056, sand bottom with sunken 
logs, W.-J. Chen, R/V Alis, 20 Dec. 2012; NTUM 10674, 
1 male, 51.8 mm SL, Papua New Guinea, West Sepik 
Province, 13.51<m northwest of Aitape, 3°03’S 142°18’E, 
370-374 m depth, St. CP 4055, Wei-Jen Chen, R/V Alis, 
20 Dec. 2012. 

3.4.2. Diagnosis 

Dorsal fin with 9-10 rays, anal fin with 8 rays, pectoral 
fin with 17-24 rays, principal caudal-fin rays 15, three 
gills, third arch with 3 gill rakers. Vertébral column with a 
total of 34-35 vertebrae, with asymmetrical latéral 
bending starting behind the skull, bent at an angle of 
85°-92°. Skull asymmetrical in frontal view. Skin naked, 
surface of head and body without striae. Disc without 
adhesive papillae; posterior margin may be partially 
fringed; région B small, left- and right-hand sides separate, 
larger on the right-hand side than on left-hand side. 




Fig. 3. Protogobiesox asymmetricus gen. nov., sp. nov., NTUM 1 0628, holotype, male, 60.0 mm SL, Papua New Guinea, Madang Province, northeast of Taviltae. 
A (upper left) Head, dorsal view, indicating latéral line canals and pores. B (upper right) head, frontal view, showing skull, eye and snout asymmetry. 
C (lower) pelvic fin dise. Bars: 3 mm. 
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Fig. 4. Protogobiesox asymmetricus gen. nov., sp. nov., NTUM 10628, 
holotype, male, 60.0 mm SL, Papua New Guinea, Madang Province, 
northeast of Taviltae. Photograph immediately after collection. Vertébral 
column bent towards the left-hand side. 



Fig. 5. Protogobiesox asymmetricus gen. nov., sp. nov., MNHN XXXX- 
XXXX, paratype, female, 59.0 mm SL, Papua New Guinea, Madang 
Province, northeast of Taviltae. Photograph immediately after 
collection. Vertébral column bent towards the right-hand side. 

3.4.3. Description 

Dorsal-fin rays x (ix-x); anal-fin rays viii (viii); 
pectoral-fin rays xxiv (xvii-xxiv); caudal-fin rays (v),xv,(iv) 
[(iv-v),xv,(iv)]. Gills 3 (3); gill rakers on 3rd arch 3 (3). 
Measurements of the holotype and paratypes, see Table 3. 


Both upper and lower jaws with several sériés of 
undifferentiated conical teeth, no canines or incisors. 

Head latéral line System with 2 pores in the nasal canal, 
1 pore in the postorbital canal, lacrimal and preopercular 
canals missing, and 2 pores in the mandibular canal 
(Figs. 2 and 3). 

Vertebrae 35 (34-35) [15 + 20 (15 + 19-20)]; asymmet- 
rical bending of the vertébral column starting from the lst 
(lst) vertebra. Pleural ribs 13 (13) pairs. 

Head broad, depressed. Head and skull asymmetrical in 
frontal view (Fig. 3B). Head length 34.5 (32.0-33.6)% SL 
(2.9-3. 1 in SL). Maximum body depth 10.7 (8.5-11.8)% SL 
(8.5-11.8 in SL). Maximum body width 19.2 (18.3-23.5)% 
SL (4.2-5.4 in SL). Maximum (horizontal) orbit diameter 
8.8 (9.1-12.2)% SL (2. 6-3.9 in head length). Snout relatively 
long, pointed (Fig. 3A). Preorbital length 7.3 (4.0-7.3)% SL 
(4.5-7. 9 in head length). Interorbital distance narrow, 2.3 
(1.1 -2.7)% SL (12.3-30.0 in head length). Upper jaw length 
13.5 (9.3-13.0)% SL. Lower jaw length 10.7 (7.9-12.4)% SL. 
Gill membranes attached to isthmus. Body bent to the left- 
hand side (males) or right-hand side (females); vertébral 
column bent at angle of 87° (85-92°), bending starts 
behind skull. Anus situated on the right-hand side of the 
body (male) or the left-hand side of the body (female); 
doser to anal-fin origin than the dise; male with a short 
urogénital papilla, its length 0.5 (0.4-0.5)% of SL (62-90 in 


Table 3 

Counts and measurements of the holotype and paratypes of Protogrammus asymmetricus n. gen. n. sp. 



Holotype 

NTUM 

10628 

Male 

Paratype 

NTUM 

10674 

Male 

Paratype 

NTUM 

10629 

Male 

Paratype 

NTUM 

10633 

Male 

Paratype 

MNHN 

2016-03333 

Female 

Standard length 

60.0 

51.8 

56.0 

53.6 

59.0 

Head length 

20.7 

16.6 

18.5 

18.0 

19.3 

Body depth 

6.4 

4.4 

5.4 

6.3 

6.0 

Body width 

11.5 

9.5 

10.8 

12.6 

11.4 

Orbit diameter 

5.3 

6.3 

5.1 

4.9 

5.8 

Preorbital length 

4.4 

2.1 

4.1 

3.9 

4.1 

Interorbital distance 

1.4 

0.7 

1.5 

0.6 

0.9 

Upper jaw length 

8.1 

4.8 

7.3 

6.7 

7.7 

Lower jaw length 

6.4 

4.1 

6.0 

5.5 

7.3 

Preanus length 

36.7 

32.9 

33.3 

31.5 

35.9 

Urogénital papilla 

0.3 

0.2 

0.3 

0.2 

- 

PC length 

5.3 

5.5 

4.6 

4.7 

6.1 

PC depth 

6.5 

4.6 

5.2 

4.9 

6.0 

PC width 

1.4 

1.6 

1.1 

2.0 

1.4 

Predorsal length 

43.9 

35.9 

42.1 

38.4 

42.4 

Preanal length 

45.6 

40.0 

41.4 

40.4 

46.2 

Distance dise to anus 

14.6 

? 

10.4 

10.9 

9.9 

Distance anus to anal-fin origin 

9.2 

8.1 

7.1 

7.7 

6.9 

Prepectoral length 

20.6 

17.0 

20.7 

20.7 

22.0 

Prédise length 

13.2 

11.0 

12.4 

11.9 

14.2 

1 st dorsal ray 

5.3 

3.2 

5.1 

4.8 

4.1 

La st dorsal ray 

6.0 

5.5 

4.0 

4.6 

5.3 

Dorsal-fin base 

14.6 

7.9 

11.0 

8.9 

12.2 

lst anal ray 

4.5 

2.5 

3.9 

2.7 

2.9 

Last but one anal ray 

6.5 

3.1 

5.8 

4.3 

4.6 

Last anal ray 

5.1 

2.4 

4.7 

3.9 

4.1 

Anal-fin base length 

9.5 

6.8 

9.5 

8.3 

9.5 

Pectoral-fin length 

10.4 

9.7 

10.7 

9.8 

10.6 

Pectoral-fin base 

3.6 

3.6 

3.9 

3.9 

4.3 

Disc length 

16.0 

? 

11.4 

11.0 

13.9 

Caudal-fin length 

12.8 

10.6 

11.4 

10.2 

12.4 


mm: measurements; PC: Caudal peduncle. 
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head length); distance between dise and anus 24.3 (16.8- 
20.3)% SL, distance between anus and anal-fin origin 15.3 
(11.7-15.6)% SL. Pre-anus length 61.2 (58.8-63.5)% SL 
(1.6-1. 7 in SL). Sides of body with an elongate dermal fold 
on left-hand side below the anal-fin base. Caudal-peduncle 
length 8.8 (8.2-10.6)% SL (9.4-1 2.2 in SL). Caudal-peduncle 
depth 10.8 (8.9-10.2)% SL (9.2-11.3 in SL). Caudal- 
peduncle width 2.3 (2.0-3. 7)% SL (26.8-50.9 in SL). Skin 
naked; surface of head and body without striae. 

Predorsal-fin length 73.2 (69.3-75.2)% SL (1.3-1.4 in 
SL). Pre-anal-fin length 76.0 (77.2-78.3)% SL (1.3-1.4 in 
SL). Dorsal fin origin on the left-hand side of the body 
(male), on the right-hand side of the body (female); anal fin 
origin on the right-hand side of the body (male), on the 
left-hand side of the body (female). Dorsal and anal-fin 
bases attached to caudal-fin base by membranes. Pre- 
pectoral-fin length 34.3 (32.8-38.6)% SL (2.6-3.0 in SL). 
Prédise length 22.0 (21.2-24.1)% SL (4.2-4.7 in SL). Disc 
length 26.7 (20.4-23.6)% SL (3. 8-4.9 in SL). Disc membrane 
inserting at the base of the 13th (13th) pectoral-fin ray. 
Disc rudimentary, incomplète, with left- and right-hand 
parts of régions A and B separate, without adhesive 
papillae; fin-ray branches at distal margins with short 
filaments; région B small (Fig. 3C). Caudal-fin length 21.4 
(19.0-21.0) % SL (4.7-5.2 in SL). 

Colour in life (Figs. 4 and 5). Head and body brown, with 
an irregular, slightly darker banding, anteriorly reddish, 
posteriorly yellowish; eyes dark grey. Pectoral-fin rays 
reddish, caudal fin rose; dorsal and anal-fin rays rose, 
membranes blackish. 

Colour in alcohol. Head and body uniformly yellowish 
brown, lower sides of head and body pale; eyes dark grey; 
occiput transparent, brain visible, yellowish white; peri- 
toneum black. Dorsal and anal fins and distal margin of 
caudal fin dark grey. 

3.4.4. Distribution 

North coast of Papua New Guinea (Madang Province; 
Sandaun Province) (Fig. 7). This new species was collected 
at depths of 381-480 m. 

3.4.5. Etymology 

The name of the new species, asymmetricus, refers to the 
unusual asymmetry of head and body. 

3.4.6. Comparisons 

This new species is distinguished from the related 
species Lepadicyathus mendeleevi by its 9-10 dorsal-fin rays 
(11-12 in L mendeleevi ), the 8 anal-fin rays (10 in L. 
mendeleevi ), 15 principal caudal-fin rays (10 in L. mende- 
leevi), the vertébral column strongly bent at an angle of 85- 
92° (40-45° in L. mendeleevi), the bending that starts behind 
the skull (just before the level of the anus in L. mendeleevi), 
and the naked skin without striae (with numerous, closely- 
set, oblique striae consisting of minute tubercles in 
L. mendeleevi ). From other gobiesocids, it is distinguished 
by the characters of the subfamily (see above). 

3.4.7. Remarks 

This species is highly asymmetrical. While the head 
is normally orientated, the skull is already slightly 


asymmetrical when seen from the front (Fig. 3B), with an 
oblique mouth slit and the left eye situated lower than 
the right eye in the male holotype; the vertébral column 
inserts slightly on the left-hand side of the skull (Fig. 2B), 
and then strongly turns left, meeting the left-hand side 
of the body on the level of the dorsal fin insertion; it 
continues running along the left-hand side towards the 
hypural plate. The anal fin, consequently, inserts on the 
right-hand side of the body; the caudal fin lies fiat on the 
ground. Asymmetries are also found in the position of 
the internai organs, which are restricted to the right half 
of the body (Fig. 2C); the anus is situated on the right- 
hand side. Asymmetries are observed also in the pectoral 
and pelvic fin skeletons; the fins and supporting 
éléments on right-hand side and are much smaller than 
the corresponding éléments on the left-hand side (dise 
asymmetry, see Fig. 3C). In the female paratype, the 
asymmetries are the other way around, with the 
vertébral column turning towards the right-hand side 
of the fish, and the anus and the internai organs shifted 
towards the left-hand side. These asymmetries are not 
artificiai, as ail four specimens of the type sériés show 
exactly the same pattern, which was already apparent in 
the freshly collected fish. 

The right pectoral fins, and right lobes of the pelvic dise, 
were partially removed for molecular studies prior to the 
morphological examination. 

This new species is found in relatively cool water in an 
area with temporary upwelling along the northern coast of 
Papua New Guinea (Kuroda [70], Lee et al. [71], Hasegawa 
et al. [72]). It was found in association with sunken logs; its 
specializations and especially the latéral asymmetry may 
be useful to be attached to these logs. 

3.5. Lepadicyathus Prokofiev, 2005 (Fig. 6) 

3.5.1. Synonymy 

Lepadicyathus Prokofiev 2005: Prokofiev [42]: 559 [546] 
(type species: Lepadicyathus mendeleevi Prokofiev 2005 by 
original désignation (also monotypic). 

3.5.2. Diagnosis 

A genus of protogobiesocine fishes with a partial 
latéral asymmetry; head slightly asymmetrical when 
seen from the front; vertébral column bent towards the 
left in males, towards the right in females, at an angle of 
40-45°; bending starts slightly before level of anus; anus 
situated towards the right-hand side of the body in 
males, towards the left-hand side in females; position 
of guts slightly asymmetrical; gills 3; pelvic fin dise 
rudimentary, in two parts, with a few rudimentary 
papillae basally in région B. Dorsal fin with 11-12 rays; 
anal fin with 10 rays; caudal fin with 10 principal rays. 
Skin naked, surface of head and body covered with a 
dense pattern of obliquely arranged striae consisting of 
minute tubercles. 

3.5.3. Description 

Dorsal-fin rays xi-xii; anal-fin rays x; pectoral-fin rays 
xxiv; caudal-fin rays (v),x,(v). No gill rakers on 3rd arch; 
1-2 gill rakers on the lst and 2nd arches. 
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Fig. 6. Lepadicyathus mendeleevi Prokofiev, 2005, ZIN 53413, holotype, female, 13.6 mm SL, Papua New Guinea, Madang Province, near Bongu. A (upper) 
Latéral view. B (centre) dorsal view. C (lower left) detail of postorbital région showing striae. D (lower right) pelvic dise. Bars: 1 mm. 


Premaxillary with a single row of incisor-like teeth, 
some with a hook-like appendage. Dentary with a single 
row of smaller incisor-like teeth, fused at the bases. 

Head latéral line System with 2 pores in the nasal canal, 
2 pores in the postorbital canal, 2 pores in the lacrimal 
canal, and 2 pores in the preopercular canal, pores in 
mandibular canal lacking. 

Pelvic dise rudimentary, in two parts; anterior margin 
with small fringes; dise only with a few rudimentary 
papillae basally in région B. 

Skin naked; surface of head and body covered with a 
dense pattern of obliquely arranged striae consisting of 
minute tubercles (Fig. 6C). 

3.5.4. Remarks 

The dermal striae found in Lepadicyathus are an 
unusual, specialised character not found in other gobie- 
socids. Prokofiev [42] compared it with Liobranchia stria of 
Briggs [3], which lias deep grooves in the skin, but the State 
of Liobranchia is quite different from that of Lepadicyathus. 
Liobranchia has just grooves surrounded by folds, while the 
skin of Lepadicyathus is even, but covered with striae 
consisting of minute tubercles. This seems to be an 
independent spécialisation, though we are unaware of 
the function of these striae. The tubercles may be based on 
reduced scales, which would support the primitive State of 
Lepadicyathus. 

3.6. Lepadicyathus mendeleevi Prokofiev, 2005 (Fig. 6) 

Mendeleev’s clingfish. 


3.6.3. Synonymy 

Lepadicyathus mendeleevi Prokofiev 2005: Prokofiev 
[42]: 560 (547) (Papua New Guinea, near village of Bongu, 
Madang Province, 1 m depth. Holotype: ZIN 53413). 
Prokofiev [73]: 181 (discussion on systematic position). 

3.6.2. Material 

ZIN 53413 (holotype, female, 13.6 mm SL), Papua New 
Guinea, Madang Province, near Bongu, 1 m depth, in 
Millepora coral, R/V Dmitrii Mendeleev, Cruise 18,14 Feb. 
1977. ZIN 53413a (1 paratype, male, 13.9 mm SL), same 
data as the holotype. 

3.6.3. Diagnosis 

Dorsal fin with 11-12 rays, anal fin with 10 rays, 
pectoral fin with 24 rays, principal caudal-fin rays 10, 
three gills, third arch without gill rakers. Vertébral column 
with a total of 33 vertebrae, with asymmetrical bending 
staring just before the level of the anus, bent at an angle of 
40-45°. Skull not asymmetrical in frontal view. Skin 
naked, surface of head and body, with numerous, closely 
set, oblique striae consisting of minute tubercles. Disc 
without adhesive papillae, except a few basally in région 
B, arranged in two rows; posterior margin not fringed; 
région B small, left and right-hand sides connected by a 
membrane. 

3.6.4. Description 

Dorsal-fin rays xii (xi); anal-fin rays x (x); pectoral-fin 
rays xxiv (xxiv); caudal-fin rays (v),x,(v) [(v),x,(v)[. Gills 3 
(3); no gill rakers on 3rd arch; 1-2 gill rakers on lst and 2nd 


R. Fricke et al./C. R. Biologies 340 (2017) 47-62 


59 


Table 4 

Counts and measurements of the holotype and paratype of Lepadicyathus 
mendeleevi Prokofiev, 2005. 



Holotype 

ZIN 

53413 

Female 

Paratype 

ZIN 

53413a 

Male 

Standard length 

13.6 

13.9 

Head length 

4.9 

4.4 

Body depth 

1.5 

1.5 

Body width 

2.9 

2.8 

Orbit diameter 

1.4 

1.7 

Preorbital length 

1.2 

1.5 

Interorbital distance 

1.0 

1.0 

Upper jaw length 

0.6 

0.6 

Lower jaw length 

0.4 

0.4 

Preanus length 

8.4 

10.3 

Urogénital papilla 

0.05 

- 

PC length 

0.5 

0.6 

PC depth 

1.1 

1.2 

PC width 

0.7 

0.8 

Predorsal length 

10.1 

10.4 

Preanal length 

9.7 

11.3 

Distance dise to anus 

1.2 

2.4 

Distance anus to anal-fin origin 

1.2 

2.45 

Prepectoral length 

5.5 

5.8 

Prédise length 

4.4 

4.1 

lst dorsal ray 

1.1 

1.2 

Last dorsal ray 

ca. 0.8 

1.0 

Dorsal-fin base 

3.0 

4.4 

lst anal ray 

1.0 

1.1 

Last but one anal ray 

1.4 

1.5 

Last anal ray 

0.9 

1.0 

Anal-fin base length 

2.4 

2.9 

Pectoral-fin length 

1.5 

1.3 

Pectoral-fin base 

1.0 

1.4 

Disc length 

2.5 

2.6 

Caudal-fin length 

1.9 

ca. 1.8 


mm: measurements; PC: Caudal peduncle. 


arches. Measurements of the holotype and paratypes: see 
Table 4. 

Premaxillary with a single row of incisor-like teeth, 
some with a hook-like appendage. Dentary with a single 
row of smaller incisor-like teeth, fused at bases. 

Head latéral line System with 2 pores in the nasal canal, 
2 pores in the postorbital canal, 2 pores in the lacrimal 
canal, and 2 pores in the preopercular canal, pores in the 
mandibular canal lacking. 

Vertebrae 33 (15 + 18); asymmetrical bending of 
vertébral column starting from llth vertebra. Pleural ribs 
12 pairs. 

Head moderately broad, depressed. Head and skull 
slightly asymmetrical in frontal view. Head length 36.0 

(31.6)% SL (2.8-3. 2 in SL). Maximum body depth 11.9 

(10.8)% SL (9. 1-9.4 in SL). Maximum body width 21.3 

(20.1)% SL (4.7-5. 0 in SL). Maximum (horizontal) orbit 

diameter 10.3 (12.2)% SL (2.6-3. 5 in head length). Snout 
elongate, slightly pointed, tip rounded (Fig. 6B). Preorbital 
length 8.8 (10.8)% SL (2.9-4.1 in head length), in the male 
not longer than in the female. Interorbital distance 7.4 (7.2) 
%SL (4.4-4.9 in head length). Upper jaw length 4.4 (4.3)% 
SL. Lower jaw length 2.9 (2.9)% SL. Gill membranes 
attached to the isthmus. Body bent to the left-hand side 
(male) or the right-hand side (female); vertébral column 
bent at an angle of 40° (45 ), bending starts just before the 


level of the anus. Anus situated in the middle between the 
dise and anal-fin origin; both sexes with a short urogénital 
papilla; distance between dise and anus 8.8 (17.3)% SL, 
distance between anus and anal-fin origin 8.8 (17.6)% SL. 
Anus situated only slightly towards the right-hand side of 
the body (male) or the left-hand side of the body (female). 
Pre-anus length 61.8 (74.1)% SL (1.4-1. 6 in SL). Sides of 
body without an elongate dermal fold below the base of 
the anal fin. Caudal-peduncle length 3.7 (4.3)% SL (23.2- 
27.2 in SL). Caudal-peduncle depth 8.1 (8.6)% SL ( 1 1 .6-1 2.4 
in SL). Caudal-peduncle width 5.1 (5.8)% of SL (17.3-19.4 in 
SL). 

Predorsal-fin length 74.3 (74.8)% SL (1.3-1.4 in SL). Pre- 
anal-fin length 71.3 (81.3)% SL (1. 2-1.4 in SL). Dorsal fin 
origin towards the left-hand side of the body (male), the 
right-hand side of the body (female); anal fin origin near 
the right-hand side of the body (male), the left-hand side of 
the body (female). Prepectoral-fin length 40.4 (41.7)% SL 
(2.4-2. 5 in SL). Prédise length 32.3 (29.5)% SL (3.1 -3.4 in 
SL). Disc length 18.4 (18.7)% SL (5.3-5.4 in SL). Disc 
membrane inserting at the base of the 12th (16th) 
pectoral-fin ray. Pelvic dise rudimentary, in two parts; 
anterior margin of région A and posterior margin of 
région B with small fringes; dise only with a few 
rudimentary papillae basally in région B. Caudal-fin length 
14.0% SL (7.2 in SL). Skin naked; surface of head and body 
covered with a dense pattern of obliquely arranged striae 
consisting of minute tubercles (Fig. 6C). 

Colour in life.-Unknown. 

Colour in alcohol.- The holotype and paratype are 
reddish brown, the head and anterior dorsal surface 
carmesine red, with 5 white stripes; one stripe running 
from the tip of the snout along the dorsal midline to the 
dorsal-fin base (Fig. 6B), a latéral stripe on each side from 
posterior margin of eye to the dorsal part of the caudal 
peduncle, and another latéral stripe on each side from 
below the dentary to the lower part of the caudal peduncle. 
Eye dark grey. Occiput not transparent. Fins reddish; 
caudal fin with a basal brown spot, central part of fin 
carmesine red. 

3.6.5. Distribution 

North coast of Papua New Guinea (Madang Province) 
(Fig. 7). This new species is known only from the type 
locality. It was collected at a depth of 1 m, inhabiting a 
Millepora coral. 

3.6.6. Remarks 

Prokofiev [42] noticed that the paratype of Lepadicya- 
thus mendeleevi was “crumpled”, but thought this to be 
artificial, and was apparently not aware of the latéral 
asymmetry. During a re-examination, however, the 
moderate asymmetry of this species was discovered both 
in the holotype and the paratype. In the female holotype, 
an ovary full of moderately large eggs was visible in the X- 
ray image. 

While the head is normally orientated, the skull is 
already slightly asymmetrical when seen from the front, 
with a slightly oblique mouth slit and one cheek larger 
than the other; in the female holotype, the vertébral 
column is bent towards the right-hand side (Fig. 6B), with 
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Fig. 7. Geographical distribution of the species of the subfamily Protogobiesocinae n. subfam. in Papua New Guinea. A. Lepadicyathus medeleevi. 
B. Protogobiesox asymmetricus gen. nov., sp. nov. 


the anus placed slightly towards the left-hand side of the 
body; the vertébral column is bent to the left-hand side in 
the male paratype, with the anus doser to the right-hand 
side. A slight asymmetry is also found in the position of the 
internai organs, which are shifted towards the left-hand 
side of the body in the female holotype, but towards the 
right-hand side of the body in the male paratype. 

In a note published subsequently to the original 
description, Prokofiev [73] noted that the absence of 
papillae on the dise of Lepadicyathus mendeleevi is not 
unique in Gobiesocidae, but resembles young specimens of 
the genus Lepadichthys. However, the holotype of L. 
mendeleevi is certainly to be considered as an adult 
specimen, as the ovary is full of eggs. 

4. Discussion 

Chen et al. [60], in a molecular phylogenetic analysis of 
the Acanthomorpha, demonstrated that the order Per- 
ciformes is poiyphyletic, and that several of the discovered 
clades need formai récognition. From the phylogenetic 
results of Chen et al. [60], the Gobiesocidae formed a 
separate clade (Clade D) together with the blennioids; this 
clade was nested within the inferred clade Q (Chen et al. 
[60]; Dettaï & Lecointre [61]) or Ovalentaria (Wainwright 
et al. [74]) of the Percomorpha contained Atherinomorpha, 


Mugilidae, Pomacentridae (“Labroidei”), Cichlidae 
("Labroidei”), Embiotocidae (“Labroidei”), Ambassidae, 
and some other perciform fishes from Pholidichthyidae, 
Pseudochromidae, Polycentridae, Plesiopidae, Opistogna- 
thidae, and Grammatidae (Near et al. [62]). Our phyloge- 
netic resuit, by including the main taxa from the clade Q, 
confirms the previous hypothesis (Fig. 1), and suggest the 
clade including Blennioidei and Gobiesocidae should be 
formally named in the future. Very recently, Nelson et al. 
[75] classified the Gobiesociformes as a separate order for 
the single family Gobiesocidae, which they considered as 
closely related to the Blenniiformes (= Blennioidei) (see Lin 
& Hastings [76]). However, the authors these studies did 
not give the reason for their rearrangement of the 
taxonomie rank at the ordinal level; the classification 
should further be settled. 

Within the family Gobiesocidae, a total of 9 subfamilies 
were previously distinguished (Briggs [3], Conway et al. 
[56]); five of these were included for the présent molecular 
phylogenetic investigation. The new subfamily Protogo- 
biesocinae is morphologically most similar to the subfam- 
ily Aspasminae. Species of the Aspasminae are distributed 
from East Africa east to Mariana Islands, north to Japan, 
south to Southern Australia and New Zealand; they are 
known from very shallow to moderately deep water. The 
deepest occurring aspasmine species are Aspasma minima 
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(Dôderlein, 1887) from Japan (183-274 m), and Modicus 
tangaroa Hardy, 1983 from New Zealand (20-149 m). 
However, our molecular results revealed that the sister- 
group of the Aspasminae is the Diademichthyinae, a group 
of the Indo-West Pacific distributed clingfishes specialized 
for a commensal relationship with certain shallow-reef 
crinoids. Interestingly, the enigmatic clingfishes from the 
subfamily Cheilobranchinae turned out to be the sister- 
group of the Aspasminae/Diademichthyinae clade in the 
inferred phylogenetic tree (Fig. 1). The species in the new 
subfamily Protogobiesocinae also possess a reduced or 
“ancestral” form of their adhesive dise formed by the pelvic 
fins, similar to the species of Cheilobranchinae. According 
to our phylogeny, this feature is apparently a derived State 
and has been evolved twice within the Gobiesocidae 
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